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Chapter 1 

Introduction 



Quantum Chromodynamics (QCD) is widely accepted as the theory of the strong 
interaction. QCD is a gauge theory for the strong interaction of quarks and gluons 
based on the SU(3) color group. Quarks have been with us since the 60's when the 
first three flavors, up (u), down (d) and strange (s), were proposed as members 
of the fundamental representation for the approximate SU(3) flavor symmetry 
observed in the hadron spectra. In the simplest picture baryons are made up of 
three quarks, while most mesons can be explained as made up of a quark and an 
antiquark. Deep inelastic scattering (DIS) on the proton confirmed this picture of 
hadrons as bound states of more fundamental particles. The data supported the 
idea of charged spin-1/2 components or partons inside the nucleon. Those partons 
were identified with the quarks. 

Quarks carry a new quantum number that was named color. The color quan- 
tum number can take on three different values and was introduced to solve the 
statistics problem associated with ground state baryons with three equal quarks. 
Experimental evidence for the color quantum number comes from the analysis of 
the a{e^e~ hadrons) /a (e"'"e~ — > fJ.^fJ-~) cross sections ratio or the — > 77 
decay width. The absence of color will render the theoretical predictions a factor 
of 3 and 9 respectively below experimental data. Hadrons do not carry the color 
quantum number, and thus quarks have to combine in such a way that hadrons 
are "colorless" or color singlets. 

Quarks have never been observed as free particles and the consensus is that 
they are confined in the physical colorless hadrons. While confinement has never 
been proven in a rigorous way in QCD, there are hints, obtained in Monte Carlo 
simulations performed in a discretized space-time (see, for example Ref. [l|), that 
QCD leads to confinement. 

In the 70's 't Hooft, Politzer 3, Gross and Wilczek 0,11] found that in non- 
abelian gauge theories the effective coupling constant decreases at large momen- 
tum transfers or short distances while it increases at low momentum transfers or 
large distances. Both features, called asymptotic freedom and infrared slavery 
respectively, are necessary to describe quark dynamics. While DIS data, where 
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large momentum transfers are involved, shows the existence of almost free par- 
tons inside the nucleon, the absence of free quarks suggests confinement at large 
distances. Fritzsch and Gell-Mann proposed that quark dynamics was de- 

scribed by a nonabelian gauge theory based on the color quantum number, QCD. 
The appropriate gauge group was color SU(3). SU(3) has a fundamental repre- 
sentation of dimension 3 that can accommodate the three values for the quark 
color quantum number. This representation is complex and the antiquark color 
states belong to the complex conjugate of the fundamental representation. The 
combination of three quarks or of a quark and an antiquark can give rise to color 
states which are SU(3) color singlets. 

The QCD gauge boson is called gluon. The gluon appears in 8 different color 
states corresponding to the adjoint representation of color SU(3). Gluons are not 
observed as free particles and as in the quark case, we think they are bound in 
the physical hadrons. Evidence for gluons comes from DIS data both from the 
observation of three jets events, and from the fact that only half of the proton 
momentum is carried by charged partons. The rest must be attributed to non 
charged partons, the gluons. 

Another important feature of QCD with three quark flavors, and which is 
essential to understand its low energy spectrum, is chiral symmetry. Chiral sym- 
metry is an approximate symmetry (it is explicitly broken by the quark masses) 
of the QCD Lagrangian. This symmetry is not shared by the ground state and 
one talks of a spontaneously broken symmetry. This property manifest itself in 
the existence of an octet of almost-massless pseudoscalar Goldstone bosons. 

With only three quark flavors, the theory of electroweak interactions involving 
quarks leads to the existence of weak neutral currents with AS* = 1. Those pro- 
cesses were not observed experimentally. The easiest way to avoid this unwanted 
feature is through the GIM mechanism (from Glashow, Illiopoulos and Maiani 0] ) 
that required a new quark flavor, the charm {cj^. This new flavor implied the exis- 
tence of new types of hadrons containing the charm quark or antiquark. The J/^ 
meson discovered in 1974 0, was interpreted as a bound state of the charm 
quark and its antiquark. Two years later came the first observation of mesons 
with open charm [111 . Two more quark flavors, later named bottom (6) and 



top (t), were predicted in 1973 by Kobayashi and Maskawa [IJ] in order to get 
a realistic model of CP violation in weak theories. The discovery of the r lepton 
in 1974-1977 led further theoretical support to the existence of the bottom and 
top quarks as the cancellation of the triangle anomaly in the electroweak theory 
required as many quarks as there were leptons (also that the former appeared in 
three colors). The T meson discovered in 1977 ij] was interpreted as being a bb 
state. Top discovery had to await for another 17 years when it was finally seen in 
proton-antiproton collisions [isi ]. 

In the large energy regime QCD can be treated perturbatively and predictions 
of the theory have shown a remarkable agreement with experiment. On the other 



^References on earlier works on charm can be found in 
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hand at low energies QCD becomes nonperturbative and it is not known how 
to solve it exactly. In this low energy regime there are approaches which make 
approximations to the exact solutions like QCD sum rules and Lattice QCD. 

Another approach to QCD in the low energy regime is the use of phenomenolog- 
ical models which are not directly derived from QCD but include some of its basic 
properties. Constituent quark models (CQM) are among these. In CQM models 
hadrons are simply modeled, as bound states of three valence quarks (baryons) 
or a quark and an antiquark (mesons). Those quarks are quasiparticle degrees 
of freedom with the same quantum numbers as QCD quarks, differing from the 
latter in their masses and in the fact that they could have an internal structure. 

Among CQM, nonrelativistic quark models (NRQM) have shown to be phe- 
nomenologically very successful. In these models constituent quarks are treated 
nonrelativistically and they interact through potentials that mimic QCD asymp- 
totic freedom and confinement. In fact it was this simple nonrelativistic picture 
that led to the necessity for the color quantum number. The vast literature on 
NRQM calculations forbids to quote here but a few works. An early success of 
NRQM calculations was the satisfactory account of the magnetic moments of the 
octet baryons [16|]. Baryon spectra was studied within the NRQM by Isgur and 



Karl in a series of papers [17j, Il8|, Il9|, |20|, |21|, [22] where it was shown that 



model with nonrelativistic point-like quarks moving in a flavor independent con- 
fining potential, perturbed by color hyperfine interactions, was able to explain the 
main features of the spectra. Meson spectra was equally well described within the 
NRQM [23^. Color hyperfine interactions among quarks are generated through the 
one gluon exchange potential first introduced in The spin-spin part of this po- 
tential explains in a simple way the difference in mass between octet and decuplet 
baryons and between pseudoscalar and vector mesons. The NRQM was extended 
to the two nucleon system in the 80's being able to explain, in a more fundamental 
picture, the short range nucleon-nucleon repulsion in S-wave scattering as due to 
one gluon exchange between quarks with a simultaneous quark exchange between 
the two three-quark clusters 24, 2^, 26]. Chiral symmetry and its spontaneously 
breaking was also incorporated into the NRQM 27| and inter-qua rk p otentials 
coming from one-Goldstone boson exchanges naturally appeared 28l. l29l. l30l| . In a 
linear realization of chiral symmetry with two quark flavors not only the Golstone 
bosons (pions) but also their chiral partner (sigma) was considered. The inclu- 



sion of the latter improved the theoretical prediction of A^A^ phase shifts [31 
and deuteron properties [s^]. Some authors advocated that the one-gluon ex- 
change part of the potential be discarded altogether in favor of only Goldstone 
boson exchanges 3j]. This exclusive view had big problems and was heavily 
critizied [sp]. Most versions of the NRQM use a mixture of both one-gluon and 
one-Golstone boson exchanges. More recently the study of dibaryons [36|, l37|. 



tribaryons [38j, |39(] and tetraquarks [40|, |4l|, I44] systems have also been addressed 



in the NRQM. For a recent review on few baryon systems studies see 



tromagnetic [4^, M, M, M, M, M, M, |5l| and weak [52, l53|, [54 




factors have also been analyzed in the NRQM and exchange current contributions 
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to these observables have been studied with detail 

In this thesis our interest is in the analysis, within the framework of a NRQM, 
of static and dynamic properties, the latter related to weak decays, of hadrons 
with one and two heavy quarks c and/or 6. In systems with a heavy quark with 
mass much larger than the QCD scale (A nnn) a new symmetry known as heavy 



quark symmetry (HQS) arises [57j, |5a, |59|, |60|, l6l|) I64]. HQS is an approximate 
SU(A'^p) symmetry of QCD, being Np the number of heavy flavors (c, b,. . . ), that 
appears in systems containing heavy quarks with masses much larger than the 
typical quantities (Aqot, m^, m^, m^, . . . ) which set up the energy scale of the 
dynamics of the remaining degrees of freedom. In that limit the dynamics of the 
light quark degrees of freedom becomes independent of the heavy quark flavor 
and spin. This is similar to what happens in atomic physics where the electron 
properties are approximately independent of the spin and mass of the nucleus (for 
a fixed nuclear charge). HQS can be cast into the language of an effective theory 
(HQET)[63| that allows a systematic, order by order, evaluation of corrections 
to the infinite mass limit in inverse powers of the heavy quark masses. HQS and 
HQET have proved very useful tools to understand bottom and charm physics and 
they have been extensivel y u sed to describe the dynamics of systems containing 
a heavy c or 6 quark 6J, |65|. For instance, all lattice QCD simulations rely on 
HQS to describe bottom systems [66|, |67|, 168], |69|, In the case of systems with 
two heavy quarks one can not directly apply HQS. It is known that the kinetic 
energy term, needed in those systems to regulate infrared divergences, breaks 
heavy flavor symmetry ff^. Still, there is a symmetry that survives: heavy quark 
spin symmetry (HQSS) ^]. This symmetry amounts to the decoupling of the 
two heavy quark spins since the spin-spin interaction vanishes for infinite heavy 
quark masses. 

The constraints imposed by HQS and HQSS have not been systematically 
employed in the context of nonrelativistic quark models. We intend to study 
hadrons with one and two heavy c and/or b quarks in a NRQM making full use 
of those constraints. For instance, those constraints will simplify notably the 
resolution of the three-body problem in baryons allowing for the use of a simple 
variational ansatz. A NRQM treatment of those systems will comply with the 
model independent predictions of HQS and HQSS in the infinite heavy quark 
mass limit, as in that limit all spin-spin interactions terms involving heavy quarks 
vanish. Another question is to what extend the deviations from the infinite heavy 
quark mass limit in a NRQM calculation agree with the predictions of HQET. 

This thesis is organized as follows. In chapter [2] we briefly introduce the quark- 
quark potentials that have been used throughout this work. 

In chapter El we study leptonic decays of pseudoscalars B, D and vector B* , D* 
mesons, and B ^ D and B D* semileptonic decays. In the latter case full form 
factors are computed and used to get differential and total decay widths. When 
possible these form factors are improved using corrections from HQET. We make 
a prediction for the value of the \Vcb\ Cabbibo-Kobayashi-Maskawa (CKM) matrix 
element which compares well with experimental determinations. 
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Also in this chapter we will use partial conservation of the axial current (PCAC) 
to relate axial current matrix elements with the pion emission amplitude and in 
this way make a prediction for the strong coupling constants gBB*n and gnvn- 

Exclusive semileptonic and nonleptonic two-meson decays of the B~ meson 
have been analyzed in chapter HI where a fairly exhaustive study of these decays 
is done. In semileptonic decays we have excluded processes that involve a b ^ u 
transition at the quark level where the NRQM fails due to the large recoil momenta 
and the ignorance of resonance exchanges. Similarly in nonleptonic two-meson 
decays we shall only consider channels with at least a cc or B final meson. 

In chapter [5] we study baryons with two heavy quarks. A procedure to solve the 
three-body problem in an approximate way using a variational ansatz is presented. 
We obtain the spectrum and wave functions that are further used to compute dif- 
ferent static properties. The infinite heavy quark mass limit of the wave functions 
is studied and as expected they comply with HQSS model independent predictions. 

In chapter E] we use the wave functions computed in the previous chapter 
to study semileptonic decays of doubly heavy baryons. Full form factors, decay 
widths, and angular asymmetries are studied. 

Finally chapter \7\ is devoted to the study of the strong one pion decay of 
baryons with a heavy quark. As in chapter [3] PCAC is used to relate weak axial 
current matrix elements with the pion emission amplitude. 



Chapter 2 



Quark quark potentials used in 
this work 



2.1 Introduction 

In this chapter we give a brief account of the quark-quark potentials that we use in 
this work to get the hadron wave functions. We shall use five different two-body 
phenomenological potentials. One of these potentials (BHAD) has been taken 



from Ref. [73] and four others (ALl, AL2, API, AP2) from Refs. [TJ, l75|]. The 
use of different potentials will allow us to check the sensitivity of our results to 
the inter-quark interaction. 

All these potentials share a common structure that includes a confining term 
plus Coulomb and hyperfine terms coming from the one-gluon exchange potential. 
They differ in the form factor used in the hyperfine and 1/r terms and/or in the 
power chosen for the confining term. All of them neglect the tensor and spin-orbit 
pieces also present in the one-gluon exchange potential @] on the account that 
they are not essential to describe the global features of hadron spectroscopy. As 
a consequence both the total orbital angular momentum L and the total spin S 
are separately good quantum numbers. Chiral symmetry is not incorporated in 
these potentials and thus no terms generated by Goldstone boson exchange are 
included. Explicit expressions appear in what follows. 

2.2 BHAD, ALl, AL2, API and AP2 potentials 

For the qq case they can be given in the general form 



(^) = 1 + Arf - A 



+ ao^^ + 1 - e-"/- a,a, (2.1) 

rmmj rr^ SmiUij \ J t^^x^ J 
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where i,j = {l{u.d), s, c,b), a are the spin PauU matrices, rrii the quark masses 
and 

-B 



Xo{mi,mj) 



, , 2m,- mi 
A' ^ 



(2.2) 



^ m j + rrij 

Parameters for the different potentials are summarized in Table 12.11 Note that 





BHAD 


ALl 


AL2 


API 


AP2 




0.52 


0.5069 


0.5871 


0.4242 


0.5743 




0. 


0. 


0.1844 


0. 


0.3466 


P 


1 


1 


1 


2/3 


2/3 


A [GeV(i+P)] 


0.186 


0.1653 


0.1673 


0.3898 


0.3978 


A [GeV] 


0.9135 


0.8321 


0.8182 


1.1313 


1.1146 


Go 


1 














ro [GeV"i] 


2.305 










k' 


0. 


1.8609 


1.8475 


1.8025 


1.8993 


rriu = rud [GeV] 


0.337 


0.315 


0.320 


0.277 


0.280 


rus [GeV] 


0.600 


0.577 


0.587 


0.553 


0.569 


mc [GeV] 


1.870 


1.836 


1.851 


1.819 


1.840 


nib [GeV] 


5.259 


5.227 


5.231 


5.206 


5.213 


B 




0.2204 


0.2132 


0.3263 


0.3478 


A [GeV^-i] 




1.6553 


1.6560 


1.5296 


1.5321 



Table 2.1: Parameters for the BHAD, ALl, AL2, API and AP2 potentials taken 
from Refs. 3, 75] 

the oo and k' terms are mutually exclusive. The power p chosen for the confining 
term is p = 1 for the BHAD, ALl, AL2 potentials. This linear confinement is 
suggested by lattice results [To]- The API and AP2 potentials use instead P = § 
which is needed, in a nonrelativistic approach, to get the correct asymptotic Regge 
trajectories for large angular momentum mesons [t?]. The hyperfine term comes 
from the one-gluon exchange potential and the original 5(r) dependence has been 
smeared out in order that nonperturbative calculations are possible. The BHAD 
potential uses a Yukawa form factor, while ALl, AL2, API and AP2 potentials use 
a gaussian form factor. In these latter cases the depth and width of the gaussian 
depends on the constituent quark masses allowing for a better reproduction of the 
hyperfine splitting of heavy mesons. The phenomenological formula in Eq. (j2.2p . 
originally proposed in 78|], seems to be well suited for that purpose [79|. The 
AL2 and AP2 potentials include a further form factor in the 1/r and hyperfine 
terms with a shape given by (1 — exp(— r/rc)). This factor simulates asymptotic 
freedom and its inclusion seems to be more important for heavy mesons than for 
light ones. Finally note that the k and k' parameters in the ALl, AL2, API and 
AP2 potentials would be the same if one simply replaced the original 5{f) function 
in the hyperfine term by the gaussian form factor. The fact that they are different 
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allows for a great improvement in the results [79] . This is also true for the BHAD 
potential where the strength for the hyperfine term, where the 6{f) function has 
been replaced by a Yukawa form factor, is six times larger than the one used in 
the 1/r term. This, in effect, is allowing for sources of hyperfine interaction other 
than the one-gluon exchange. 

All free parameters in the potentials were adjusted in the original works to 
reproduce the light and heavy-light meson spectra (AL1,AL2,AP1,AP2), or to fit 
the low-lying levels of charmonium (BHAD). In this latter case nib ^-nd were 
adjusted to reproduce the ground state of the T and cj) mesons respectively while 
the light quark mass rrii was chosen from magnetic moments considerations. 

To obtain the quark-quark potential we shall use, as has been done in Refs. 0, 



7J, |75|, the prescription 



= ^ (2.3) 



This relation assumes the potential has a global color dependence given by , 
where A*^ stands for the Gell-Mann matrices in the color space. 



Chapter 3 



Leptonic and semileptonic 
decays of mesons with a heavy 
quark 



3.1 Introduction 

This chapter is devoted to the calculation of different weak observables of pseu- 
doscalar and vector mesons with a heavy quark. These weak observables are of 
great interest since they help to probe the quark structure of the involved hadrons, 
and also provide information to determinate the absolute value of elements of the 
CKM matrix. 

In Ref. [so'], in which a similar model was used to study the — > l~i>i and 

— > I'ui reactions, it was shown that a direct nonrelativistic calculation does 
not meet HQET constraints. This problem was solved imposing HQET relations 
among form factors. That led to a determination of the \ Vcb\ CKM matrix element 
in good agreement with experimental results. 

One would expect this problem will also show up in this calculation, so that 
constraints coming from HQET will be included at some point in the calculation 
in order to improve the results. 

In the case of mesons with a heavy quark, HQS leads to many model in- 
dependent predictions. For instance in the HQS limit the masses of the lowest 
lying (s-wave) pseudoscalar and vector mesons with a heavy quark are degener- 
ate. Nonrelativistic quark models satisfy this constraint: the spin-spin terms, that 
distinguishes vector from pseudoscalar, are zero if the mass of the heavy quark 
goes to infinity. HQS also predicts that the masses and leptonic decay constants 
of pseudoscalars and vector mesons are related via 



fpmp 



= fvmv 

HQS 



(3.1) 

HQS 



relation that is also satisfied in the quark model in the HQS limit. If one looks 
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now at the form factors for the semileptonic B ^ D and B D* decays HQS 
predicts relations among different form factors that are also met by the quark 
model in the HQS limit. The question is to what extent the deviations from the 
HQS limit evaluated in the quark model agree with the constraints deduced from 
HQET. In addition we will make use of these HQET constraints to improve the 
quark model results and thus come up with reliable predictions. 

Apart from lattice QCD and QCD sum rules (QCDSR) calculations with which 
we shall compare our results, and that will be quoted in the following, the different 
observables analyzed in this work have been studied in the quark model starting 
with the pioneering work of Ref. 



81 



83( 1 within a non relativistic version, 



to continue with different versions of the relativistic quark model applied to the 



determination of decay constants 
factors and differential d ecay widths [9 



form 



Isgur-Wise fu nctions |93.ll03 . 



constants fill. Eoi.El 



112 



104 



113 




114 ]. 



HLu ^jj, i^ij,. i^m. i-'^- • l_itiii > l -i-wvji , LLUJj, l -i-w^ J , 

, liod, IiOTI, llOa, llOa, lud] or strong coupling 



For cal culat ions in this chapter we have used experimental meson masses taken 
from Ref. |lli |. 

3.2 Meson states 



For a meson M we use the following expression for the wave function in the Fock- 
space representation. 



M,XP 



NR 



i27r)-2y^2EjMWfM) 



q, Qi pi 



q, a2 P2 



' P-p 

P + p 



(3.2) 



where P stands for the meson three momentum and A represents the spin projec- 
tion in the meson center of mass, ai and 02 represent the quantum numbers of 
spin (s), flavor (f) and color (c) 



a 



{s,f, c) 



(3.3) 



of the quark and the antiquark, while Ef-^(pi), pi and i^jjfe), P2 are their respec- 
tive energies and three-momenta, is the mass of the quark or antiquark with 
flavor /. The factor (— 1)2~'*2 jg included in order that the antiquark spin states 
have the correct relative phas^. The normalization of the quark and antiquark 



^Note that under charge conjugation (C) quark and antiquark creation operators are related 
via CcJ,(p)C^ = (— 1) 2 ). This means that the antiquark states with the correct spin 

relative phase are not |0) — \q, a p ) but are instead given by (— 1) rft (p) |0) = 

(-1)4-^ I g, ap). 
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states is 

{a' p'\ap) = (5,,, „ (2^)3 2Ef{p) 6{p'-p) (3.4) 

Furthermore, (t)2f,a2{p) is the momentum space wave function for the relative 
motion of the quark-antiquark system. Its normahzation is given by 

O-l 0-2 

and, thus, the normahzation of our meson states is 

( M, A' P' I M, A P ) = 5y , A (2vr)3 5(P' - P ) (3.6) 

For the particular case of ground state pseudoscalar (P) and vector iV) mesons 
we can assume the orbital angular momentum to be zero and then we will have 



Sl,S2,0 



si,S2, A (3.7) 



= 1^5,^,,A-^)^^^M\p\)ymCp) 

where (ji, J2) Jalw-i, ""^2, ""^s) is a Clebsch-Gordan coefficient, Iqo = l/\/47r is the 
Z = m = spherical harmonic, and (j) j^^j^ d I) Fourier transform of the radial 
coordinate space wave function. The phases are introduced for later convenience. 

To evaluate the coordinate space wave function we shall use the potentials 
presented in chapter [2j This will provide us with a spread of results that we will 
consider, and quote, as a theoretical error to the averaged value that we will quote 
as our central result. Another source of theoretical uncertainty, that we can not 
account for, is the use of nonrelativistic kinematics in the evaluation of the orbital 
wave functions and the construction of our states defined above. While this is a 
very good approximation for mesons with two heavy quarks (as the Be that will 
be studied in chapter H]), it is not that good for mesons with a light quark. That 
notwithstanding note that any nonrelativistic quark model has free parameters 
in the inter-quark interaction that are fitted to experimental data. In that sense 
we think that at least part of the ignored relativistic effects are included in an 
effective way in their fitted values. 



3.3 Leptonic decay of pseudoscalar and vector B and 
D mesons 

In this section we will consider the purely leptonic decay of pseudoscalars (P, D) 
and vector (P*, D*) mesons. The charged weak current operator for a specific 
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pair of quark flavors /i and /2 reads 



(ci,Si), (C2,S2) 



^ai(0)7^(l-75)^a2(0) (3.8) 



with ^ctx a quark field of a definite spin, flavor and color. The hadronic ma- 
trix elements involved in the processes can be parametrized in terms of a unique 
pseudoscalar fp or vector fy decay constant as 



(o| j/^/^o) P,p) = (o| -JiV^^(o) 

7/1/2 



p.p 



(0| j/i^^(0) y,AP\ = (0| 4y^(0) V, \P)=e(^)^{P)mvfv (3.9) 



where the meson states are normalized such that 

M, A' P' I M, A P \ = ,5v, A (2vr)3 2^m(P) 8{P' - P) 



(3.10) 



with Enj{P) the energy of the meson. 

In the first of Eqs. (j3.9p P^ is the four-momentum of the meson, while in the 
second my and e(^x^ ( P ) are the mass and the polarization vector of the vector 
meson. In both cases /i and /2 are the flavors of the quark and the antiquark 
that make up the meson. Expressions for the different polarization vectors used 
in this memory are presented in appendix lAl 

Concerns about the experimen tal d etermination of the pseudoscalars decay 
constants have been raised in Ref. 116l |. There the effect of radiative decays was 
analyzed concluding that for B mesons the decay constant determination could 
be greatly affected b y rad iative corrections. In the vector sector, and as rightly 
pointed out in Ref. 117l |. the vector decay constants are not relevant from a 
phenomenological point of view since B* and D* will decay through the electro- 
magnetic and/or strong interaction. They are nevertheless interesting as a mean 
to test HQS relations. 

For mesons at rest we will obtain 



fp = —{0\ J{^J%0) 
mp 

/y = — (0| 
my 



P,0 
V, 00 



(3.11) 



with mp the mass of the pseudoscalar meson. In our model, and due to the 
different normalization of our meson states, we shall evaluate the decay constants 
as 



fv = -\ — 

my 



p,o) 

/ NR 

V, 00) 

/ NR 



(3.12) 
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The corresponding matrix elements are given in appendix [Bl 

The results that we obtain for the different decay constants appear in Tables [3T] 
and 13.21 Starting with fu and Jd^ our results are larger that the ones obtained 
in the lattice by the UKQCD Co llaboratio n 118l | or the ones evaluated using 



QCD spectral sum rules (QSSR) [ui, Not only the independent values 

are larger but also the ratio fos/ fo is larger in our case. On the other hand 



our results are in better agreement with other lattice determinations 12ll . Il22l | 
They als o coinpar e very wel l with the experimental measurements of /d and /d 



in Refs. [lii, [lii, EH H [ll^, As for /s and /s,, we find a very good 



agreement in the case of between our results and the ones obtained in the 
lattice or with the use of QSSR. For fs our result is smaller and then also our 
ratio JbsI fs is larger. 

For the vector meson decay constants we obtain the values 



fD* = 2231?^ MeV fn- = 326+1] MeV 
fB* = l5lt\lMeV fBt 



(3.13) 



which are very much the same as the values obtained for the decay constants of 
their pseudoscalar counterparts. This almost equality of pseudoscalar and vector 
decay constants is expected in HQS in the limit where the heavy quark masses go 
to infinity where one would have j6l|, 



fv mv = fp mp 



my = mp 



(3.14) 



Our decay constants satisfy the above relation within 2%. On th e oth er hand 
UKQCD lattice data show deviations as large as 20% for D mesons lisl ]. 

In orde r to compare the values of the vector decay constants with lattice data 
from Ref . [118l | we give in Table 13.21 the quantity fv = my / fy . We find good 
agreement for fu* and Jb* but not so much for the other two. Also our ratios 
fo*/ foi and fB*/fB* sue larger than the ones favored by lattice calculations. 

On the other hand the ratio 



1.138 



rO.Oll 

-0.008 



(3.15) 



is in very good agreement with the expectation in Ref. 133l | where they would get 
1.05 ~ 1.20 for that ratio. 



3.4 Semileptonic decay of B into Dli/ and 

In this case the hadronic matrix elements are parametrized as 

(^D, P' \r^'{0)\ B, P) (^D,P' 

^mBmo ^JmBmo 

= {v + v')^h+{w) + {v-v')^h^{w) (3.16) 



B, P 



22 CHAPTER 3. DECAYS OF B,D,B* AND D* 





fD [MeV] 


/o. [MeV] 


foJfD 


This work 


2431^^ 


34lt^ 




Experimental data 








CLEO 
ALEPH fl25] 
OPAL [126] 
BEATRICE \127] 
E653 fl29] 
BES fl28] 


222.6 ± 16.7tl:l [1231 
— 

371+1?^ ± 25 


280 ± 19 ± 28 ± 34 [124] 

285 ± 19 ± 40 

286 ± 44 ± 41 
323 ± 44 ± 12 ± 34 
194 ± 35 ± 20 ± 14 


— 
— 


Lattice data 








UKQCD [1181 
Fermilab Lattice [1211 
M. Wingate et al. [122] 


201±3±6±9± 13 

— 


229(3)li2 
249±3±7±11±10 
290 ± 20 ± 29 ± 29 ± 6 


— 


QCD Spectral Sum Rules 








S. Narison: [119, 1201 


203 ± 23 


235 ± 24 


1.15 ±0.04 




fB [MeV] 


/s, [MeV] 


fsJfB 


This work 


155+15 

— 12 


239t? 


1.54+°-»g 

— U.Uo 


Experimental data 








BELLE [m] 


99Q+36 +34 






Lattice data 








UKQCD [118] 
M. Wingate et al. [1221 
Lattice world averages 


200 ± 30 [131] 


220(6)tig 
260 ± 7 ± 26 ± 8 ± 5 
230 ±30 [132] 


1.13(1)11 
1.16 ± 0.04 [1311 


QCD Spectral Sum Rules 








S. Narison [119. 1201 


207 ± 21 


240 ± 24 


1.16 ±0.04 



Table 3.1: Pseudoscalar fp decay constants for B and D mesons 
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J 


fo-lfD" 


This work 
UKQCD [1181 


q 1+0.9 
8.6(3)1^ 


D-O_0.4 

8.3(2)1^ 


1 41+0.06 
-'■•^-'--0.05 

1-04(1)12 




/b. [MeV] 


hi [MeV] 


/b* / fsi 


This work 
UKQCD [118] 


35.6111 
28(1)11 


23.0+|-° 
25(1)1^ 


1 c;c;+0.07 
■■-•"3^-0.06 

1.10(2)12 



Table 3.2: fy = my / fv for B* and D* mesons 



D*,\P' \ j;/(o) I B, P 



-i (£(A)m(^^'))* {vo + l)hAAw) 

+ i [^{X) {P')) -v [v^ {w) + v'^. Ha., {w)) 



(3.17) 

where v = P/ms and v' = P'/mD,D* are the four velocities of the initial B and 
final D, D* mesons, w = v • H and e^iyap is the fully antisymmetric tensor with 

In the limit of infinite heavy quark masses mc, rrib — oo HQS reduces the 
six form factors to a unique universal function ^{w) known as the Isgur-Wise 
function [6l|, l62| 



h+{w) = /lAi(u') = hA-iiw) = hv{w) = ^{w) 
h-{w) = hA2{w) = 



(3.19) 
(3.20) 



Vector current conservation in the equal mass case implies the normalization 

e(i) = 1 (3.21) 

Away from the heavy quark limit those relations are modified by QCD corrections 
so that one has 



hj{w) = {aj + Pj{w) + -fj{w) + Oil/ml^)) ^{w) 



(3.22) 



w is related to the four momentum transferred square via 



q — rrig + m^y u* — 2 wniB mo, D* 



(3.18) 
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w 


P+ 


P- 


PV 


PAx 


PA2 


R A 

PA-i 


1 n 

i.U 


9 fi 
Z.D 


— 0.41 


1 1 Q 

li.y 


1 ^ 
— 1.0 


1 1 n 


9 9 


1.1 


-0.3 


-5.4 


8.9 


-3.8 


-10.3 


-0.2 


1.2 


-3.1 


-5.3 


6.1 


-5.9 


-9.8 


-2.5 


1.3 


-5.6 


-5.3 


3.5 


-7.9 


-9.3 


-4.6 


1.4 


-8.0 


-5.2 


1.1 


-9.7 


-8.8 


-6.6 


1.5 


-10.2 


-5.2 


-1.1 


-11.5 


-8.4 


-8.5 


1.59 


-12.1 


-5.1 











Table 3.3: QCD corrections Pj{w) in % as evaluated in Ref. |l35l | 



The aj are constants fixed by the behavior of the form factor in the heavy quark 
hmit 



aAi 

OiA2 



«A3 





ay = 1 



(3.23) 



The different (3j account for perturbative radiative corrections 134i | while the 7^ 
are non pertur bati ve in nature and are proportional to the inverse of the heavy 
quark masses 
restriction. 



135( 1 . At zero recoil {w = 1) Luke's theorem jl36t | imposes the 



7+(l)=7Ai(l) = 



(3.24) 



so that power corrections to /i+(l) and (1) are of order 0{l/m^ j). In Tables [3^31 
and 13. 4| we collect the values for the different f5j and 7^- in the full interval of w 



values allowed in the two decays. These two tables have been taken from Ref. 135l |. 



w 


7+ 


7- 


IV 


7Ai 


IA2 


IA3 


1.0 


0.0 


-4.1 


19.1 


0.0 


-23.1 


-4.1 


1.1 


2.7 


-4.1 


20.7 


2.9 


-21.4 


-0.7 


1.2 


6.2 


-4.1 


23.1 


6.5 


-19.8 


3.4 


1.3 


10.5 


-4.2 


26.3 


10.7 


-18.3 


8.0 


1.4 


15.3 


-4.4 


30.0 


15.4 


-17.0 


13.0 


1.5 


20.6 


-4.5 


34.3 


20.5 


-15.8 


18.5 


1.59 


25.7 


-4.7 











Table 3.4: Power corrections ^j{w) in % as evaluated in Ref. [13: 
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3.4.1 B^DW decay 



Let us start with the B ^ Dlu case. In the center of mass of the B meson and 
taking P' = —q = —\q\k in the z directiorH we will have for the form factors 
hj^{w) and 



h+{w) 
h-{w) 
where E£){—q) 



1 



1 

^/2mB2m£) 



V'i\q\)+^^^^l^iEn{-q)-mD) 

\q\ 



m?j^ + 9^ and l^'^(lg'l) = 0,3) is given by 



D, 



In our model 1^^(|^|) is evaluated as 



^2mB2ED{-q) ^^{D,-\q\k 



B, 



(3.25) 



(3.26) 



B,0) 

/ NR 



(3.27) 



which expression is given in appendix O 



In the case of equal masses nib = "^c vector current conservation demands that 



h+il) = l 







(3.28) 



In this limit we find that /i+(l) = 1 so that our value for /i+(l) complies with 
vector current conservation. On the other hand h-{w) ^ violating vector current 
conservation. 

In the top panel of Fig. 13.11 we show the values of h+{w) and h^{w) for the 
B ^ D transition as obtained from Eqs. (13.251 I3.27P with the use of the ALl 
interquark potential. The values for h-{w) are not reliable. Actual calculation 
shows that they are of the same size as the deviations from zero that one computes 
in the equal mass case. To improve on this what we shall do instead is to use the 
form factor h^{w) and Eq. (13.220 to extract ^{w) (we shall call it £,+ {w) ) and from 
there we can re-evaluate h-{w) with the use of Eq. (|3.22p . The results appear in 
the middle panel of Fig. [XTl where we al so sh ow the lattice results for ^(w) obtained 
by the UKQCD Collaboration in Ref. 1371 ]. We find good agreement with lattice 
data. Finally in the lower panel of Fig. 13.11 we show the different £,+ {w) obtained 
with the use of the different interquark potentials. As we see from the figure all 
(,+ {w) are very much the same in the whole interval for w. 

The slope at the origin of our Isgur-Wise function is given by 



P 



1 dc+H 



dw 



0.35 ± 0.02 



(3.29) 



We use k for the unit vector in the z direction. 
*In this case w is related to \q \ via \q \ = mo 
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Figure 3.1: Top panel: h^{w) and h-{w) for the B ^ D transition as obtained 
from Eqs. (|3.25l[3r27j) using the ALl interquark potential. Middle panel: h^{w) as 
before, ^+{w) obtained from the values of h^{w) using Eq. (I3.22p . h-{w) obtained 
from using Eq. (j3.22p . We also show the UKQCD lattice data by Bowler 



137l |. Lower panel: Different 5,+ {w) obtained with the above procedure for 



et al. 

the different interquark interactions. Lattice data is also shown. 
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small compared to the lattice value of 
data. 



0.81 



+17 
-11 



extracted from a best fit to 



Differential decay width 

Neglectin g lep ton masses the differential decay width for the process B ^ Dlu is 
given by [1381 ] 



(3.30) 



where Gp = 1.16637(1) x 10 ^ GeV ^ [115*] is the Fermi decay constant, Vet is 
the CKM matrix element for the b ^ c weak transition, and Fo[w) is given by 



Fd(w) 



h+{w) 



1 + r 



h-{w) 



(3.31) 



with r = m£, /m b ■ 

In Fig. 13.21 we show our calculation for Fo{w) \Vcb\ obtained with the ALl 
interquark potential and using three different values of \ Vcb\ corresponding to the 
central and extreme values of t he ra nge for \ Vcb\ favored by the Particle Data Group 
(PDG), \Vcb\ = 0.039 ~ 0.044 [ni]. We also show the experimental data for the 



decays B' ^ D^lv and B^ ^ D+lv obtained by the CLEO Collaboration [l39l |. 
a fit to CLEO data using the form factors of Boyd et al. [140], and the ex perim ental 
data for the decay B^ — > D~^ I v obtained by the BELLE Collaboration 1411 ]. Our 
results are larger than experimental data for w > 1.2. Our total integrated width 
will thus be larger than the experimental one for any reasonable value of \Vcb\- 
From our data we extract the slope sd, w = 1 given by 



1 dFoiw) 



Fd{w) dw 



0.38 ± 0.02 



(3.32) 



w=l 



which is small compare d to the values extracted f rom the experimental data: 
pI = 0.76 ± 0.16 ± 0.08 fliiif] and pj^ = 0.69 ± 0.14 jliU obtained from a linear 
fit to the data, or pj^ = 1.30ifc0.27ifc0. 14 [liil and = 1.16±0.25 [141] obtained 
using the form factors of Boyd et al. 140 ]. Thus, only our results close to w = 1 
seem to be reliable. We can use our prediction for -Fd(I) to extract the value of 
\Vcb\ from the experimental determination of the quantity ] V^fc]FD(l). Different 
values of that quantity appear in Table 13.51 

Our result for -Fd(I) is given by (we do not show the theoretical error which 
is or the order of 10~^) 



Fd{1) = 1.04 



which is in good agreement with other calculations -Fz)(l) 
Fd{1) = 1.04 [83i] or Fd{1) = 1.069 ± 0.008 ± 0.002 ± 0.025 fl43' 



(3.33) 

0.98 ± 0.07 [l4^ . 

From our value 
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0,06 r 



0,05- 



> 



0,04^ 



0,03 



0,02 



- - This work using IV^|^I=0.044 

— This work using IV^|^I=0.0415 
. . . . This work using IV^|^I=0.039 
• Exp. (CLEO Coll.) 

■-■ Fit to CLEO data 
O Exp. (BELLE Coll.) 




1.1 1,2 



1.3 
CO 



1,4 



1.5 



Figure 3.2: F£){r,w) \Vch\ obtained with the ALl interquark potential. Sohd line: 
our results using |V"cfe| = 0.0415. Dashed line: our results using \Vcb\ = 0.044. 



Dotted line: our resu 
CLEO Collaboration 
factors of Bo yd et al. 
oration [l4l|] . 



ts u sing \Vcb\ = 0.039. Circles: experimental data by the 



139 1. Dashed-dotted line: fit to CLEO data using the form 



140t |. Diamonds: experimental data the by BELLE Collab- 





\Vcb\FD{l) 


CLEO Collaboration [139] 
BELLE Collaboration [141] 


0.0416 ± 0.0047 ± 0.0037 
0.0411 ± 0.0044 ± 0.0052 



Table 3.5: \Vcb \ Fd{^) values obtained by different experiments. 



for i*D(l) and the experimental values for |14b|FD(l) we can obtain \Vcb\ in the 
range 



\Vcb\ = 0.040 ±0.006 



(3.34) 



This result agrees with a recent determination based on the analysis of the A;, K g I vi 
reaction using the same model, from where \Vcb\ = 0.040ib0.005 was obtained 1441 ]. 



3.4.2 B ^ D*lu decay 

Working again in the center of mass of the B meson and taking P' = —q = — \q\ k 
in the z direction we will have for the form factors hv{w), hji-^{w), hji^iw) and 
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hA^iw) the expressions 



hv{w) = 
h-Aiiw) = 

hA2 (w) = 

hAsiw) = 

with Ed* (- 



V2 



rriB 



y/2 1 

,(*) 



rriD* 
niB 



(*) 



A'o/om En*{-q)A\;/,{\q 



+ 



V2 



m 



4\im . V2 Eo*{-q)A':i,m) 

"r I -*io 



.(*) 



D* 



y/mBmo* 



\q\' 



mo* 



(3.35) 



m 



D* 



and where I) and A^^^^{ 



q\) are given by 



In our model and |) are evaluated as 

= ,j2mB2ED*{-q) j^j^[d*,X -\q\k 
= .j2mB2ED*{-q) ,rniD\\ 



D\ A 
D\ A 



1^1 k 
\q \ k 



1')m(0) 



B, 
B, 



(3.36) 



A 



A, /X \ 



NR 



k 



1')m(o) 



B, O) 

/ NR 

B, o) 

I NR 



(3.37) 



with expressions given in appendix [Dl 

In the top panel of Fig. 13.31 we show our results for the hv{w), Ha^ (w), Ha^ (w) 
and hA:j{w) form factors, obtained with the ALl interquark potential and the use 
of Eq. ()3.35p . In the lower panel of the same figure we show the ratios 



Riiw) 



hv{w) 
hAiiw) 

hA^ jw) + r hA2 jw) 
hAiiw) 



(3.38) 



where now r = mD*/'mB. These ratios are expected to vary very weakly with w. 
We find indeed that this is so in our case being our values of Ri{w) and R2{w) 
within 4% of unity. In Table [3^ we give now our results for -Ri(l) and i?2(l) and 
compare them to different experimental! and theoretical determinations. We find 
discrepancies of the order of 15 ~ 33% for Ri{l) and 13 ~ 46% for i?2(l)- 



^The experimental results by the CLEO and BABAR Collaborations have been obtained with 
the assumption that Ri{w) and R2{w) are constants. 
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1 ' 1 ' 

ALl 






1 . 1 . 1 . 1 



CO 



1 ' 

ALl 




- R2(co) 
Ri(co) 





1 1,2 1,4 

CO 



Figure 3.3: Top panel: hv{w), hAi{w), hA2{w) and hA^iw) form factors obtained 
using Eq. (j3.35p . Lower panel: Ri{w) and R2{w) ratios. In both panels the ALl 
interquark potential has been used. 



One can understand these discrepancies by evaluating the different S,{w) func- 
tions obtained from the form factors with the use of Eq. (I3.22p and the P and 7 
coefficients of Neubert given in Tables 13.31 and 13. 4[ The results appear in the top 
panel of Fig. 13.41 One can infer from the figure that our results for hA-^iw) are 
not reliable. Also we somehow miss the correct normalization for hv{l)- On the 
other hand the values of (w) and Ea ^ (w) are equal within 4% and in reasonable 
agreement with lattice data from Ref. 137l |. 

To improve the nonrelativistic quark model prediction, and similarly to what 
we did in subsection 13.4.11 we will take S^Aiiw) as our model determination of the 
Isgur-Wise function S,{w) and we will reevaluate the form factors with the use of 
Eq. (I3.22p . What we obtain is now depicted in the middle panel of Fig. 13.41 In 
the lower panel we give the different (w) obtained with the different interquark 
potentials. They do not show any significant difference. 
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iii(l) 


R2{1) 


This work 


1.01 ± 0.02 


1.04 ± 0.01 


CLEO [145] 


1.18 ±0.30 ±0.12 


0.71 ±0.22 ±0.07 


BABAR (Preliminary) [146] 


1.328 ± 0.055 ± 0.025 ± 0.025 


0.920 ± 0.044 ± 0.020 ± 0.013 


Caprini et al. [147] 


1.27 


0.80 


Grinstein et al. [1481 


1.25 


0.81 


Close et al. [103] 


1.15 


0.91 



Table 3.6: Ri{l) and i?2(l) 



The slope of the ^yi^(w) function at the origin is given by 



0.55 ±0.02 



(3.39) 



to be compared to the lattice result p"^ = 0.931*^59 137]. In this case we are within 
lattice errors, but one can not be very conclusive due to the large value of the 
latter in this case. 

Finally in Fig. 13.51 we give the ratio ^^{w)/^Aiiw) evaluated with the ALl 
interquark potential. We see the differences between the two Isgur-Wise functions 
are at the level of 3-7%. 



Differential decay width 

Neglecting lepto n masses the differential decay width for the process B 
is given by [149l | 



= ^ {niB - mn^ f m^, y^^^I^ {w + 1)^ 



1 + 



4:W 1 — 2wr + r 



w + 1 {1- r)' 



Fl,{w) (3.40) 



where Fd* {w) is defined as 



Fd*{w) = hA^w)^ 



Hl{w) + Hliw) + Hl{w) 



\ w+l 



Aw l—2wT+r^ 



(l-r)2 



(3.41) 



The Hj{w) are helicity form factors given in terms of the Ri{w) and R2iw) ratios 
as 



w — 1 

Ho{w) = 1 + - [1-R2{w) 

1 — r 



H±iw) 



\/l — 2wr + 
1 — r 



(3.42) 
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Figure 3.4: Top panel: different functions obtained from the hv{w), hA^iw) 
(w) and (w) form factors using Eq. (I3.22|) an d the ALl interquark potential 



Lattice data by K. C. Bowler et al. from Ref. 1371] are also shown. Middle 
panel: form factors obtained from ^Ai{w) with the use of Eq. (I3.22|) . Lower panel: 
Different ^aA''^) obtained with the different interquark potentials. 
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Figure 3.5: Ratio of our two Isgur-Wise functions calculated with the ALl in- 
terquark potential. 



Similarly to Fig. [321 in Fig. l3.6l we show our results for the quantity Fd* (w) \ Vcb\ 
evaluated with the ALl interquark potential and using the values of \Vcb\ corre- 
sponding to the central and extreme values of the range for \Vcb\ favore d by the 
PDG. We also show the experimental data by the CLEO Collaboration jl45l | for 
the B~ — > D*^lv reaction (squares), and for the — > D*^lv reaction (circles) 
together with a best fit, and the experimental data by the BELLE Collabora- 
tion |15Q] for the B^ — > D*~^ I v reaction (diamonds). We find good agreement with 
CLEO data for small w values. Disagreement starts already at around w = \.\ 
where our results start to go above the experimental data. BELLE data are sys- 
tematically below our results. 

Also our slope at the origin 



2 _ 1 dFD*{w) 

Pd* 



Fd* (w) dw 



0.31 ± 0.02 (3.43) 



=1 



issmaller than the value obtained by the BELLE Collaboration p^, = 0.81 it 0.12 
|l50| using a linear fit to their data. All this means that our total width would 
be larger than the experimental one for any reasonable value of Vcb. On the 
other hand experimentalists are able to extract the value of \Vcb\ Fj:)*(l). Different 
experimental results for that quantity appear in Table 13.71 
Our result for i<£)*(l) is given by 



Fd* (1) = /lAi (1) = 0.983 ± 0.001 



(3.44) 
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* 
C 



0,05 



0,04 



0,03 



0,02 



0,01 




_ - This work using IV_.j^l=0.044 

— This work using IV^|^h0.0415 

— This work using IV^j^l=0.039 
• Exp. D** 1 V (CLEO Coll.) 
■ Exp. D*» 1 V (CLEO Coll.) 

■- Fit to CLEO data 
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CO 



Figure 3.6: FD*{r,w) \Vcb\ obtained with the ALl interquark potential. SoHd line: 
our results using \ Vcb\ = 0.0415. Dashed line: our results using \ Vcb\ = 0.044. Dot- 
ted line: our results using {V^-hl = 0.039. Circles and squares: experimental data 
by the CLEO Collaboration i_145:!]. Dashed-dotted line: fit to CLEO Col laboration 
data. Diamonds: experimental data by the BELLE Collaboration 14l|. 





\Vcb\FD^{l) 


(CLEO Coll.) 


[145] 


0.0431 


±0.0013 ±0.0018 


(DELPHI Coll.) 


[151] 


0.0392 


±0.0018 ±0.0023 


(BELLE Coll.) 


[150] 


0.0354 


±0.0019 ±0.0018 


(BABAR Coll.) 


[152] 


0.0355 


±0.0003 ±0.0016 



Table 3.7: \Vcb\ Fd*{1) values obtained by different experiments. 

Comparison with the experimental data for \ Vcb\ Fd* (1) allows us to extract values 
for \Vcb\ in the range 

\Vcb\ = 0.0333 ~ 0.0461 (3.45) 

One can not be more conclusive due to the dispersion in the experimental data 
for \Vcb\ Fd'{1). From DELPHI data alone we would obtain \Vcb\ = 0.040 ± 0.003 
in perfect agreement with our determination using the B ^ D reaction data. We 
should also say that our value for Fc)*(l) is larger than the lattice determination 
F£)*(l) = 0.919l^Qg35 by S. Hashimoto et al. [153] normally used by experimental- 
ists to extract their \Vcb\ values. A new unquenched lattice de term ination of this 
quantity by the Fermilab Lattice Collaboration is in progress 1321 ]. 
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3.5 Strong coupling constants qh* 



Hit 



In this section we will evaluate the strong coupling constants gH*Hn where H 
stands for a. B oi D meson. To this end we shall make use of the partial conser- 
vation of the axial current hypothesis (PCAC) which relates the divergence of the 
axial current to the pion field as 

d^jZ{x)=Uml'^^-{x) (3.46) 



where = 130.7±0.1±0.36 MeV [llSl is the pion decay constant, = 139.57 MeV 
115l | is the pion mass, and (x) is the charged pion field that destroys a vr" 



and creates a tt~^. Using the Lehmann-Symanzik-Zimmermann reduction formula 
one can relate the matrix element of the divergence of the axial current to the 
pion emission amplitude as 

{H, P'\q>^jd^l(0)\H*,XP) = -iU^^A^^l^^iP',P) (3.47) 

q 

where g = P-P' and ^iJ^j^^(P', P) is the pion emission amplitude for the process 
H* Hn given 

A^hIh.{P',P) = -gH'Hn{q') {q-eix){P)) (3.48) 

The matrix element on the left hand side of Eq. (|3.47|) has a pion pole contribution 
that can be easily evaluated to be 

2 

{H, P' I ji^(O) I H* , XP)pion-pole = ~^/7r ^2 _ ^2 '^if*H7r(-^'' 

(3.49) 

SO that we can extract a non-pole contribution 

{H, P'\q^'jil{0)\H*,XP)non-pole = A^iilH^{P',P) 

(3.50) 

which is the one we shall evaluate within the quark model. For the matrix element 
on the left hand side of Eq. (|3.50p we can use a parametrization similar to the one 
used in Eq. dXTTD 

{H, P' {qf^ J'}^^^{0)\H* ,XP)non-pole = 

q''^-i£iX)f,{P) (w + l) HaAw) 

+i (^e(A) {P ) • v'^ {v'^ (w) + Vfj, /lAg (w)) > ^/mJrmF (3.51) 



"Corresponding to the emission of a tt^ 
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with the result that 



gH*H-K{l^) = -^1 + 1) hA^{w) + w ( — ^ /i^2(^) - ^Aziw] 



niH* 
rriH 



(3.52) 



The evaluation of the form factors is done is a similar way as the one described in 
subsection I3.4.2[ The results that we get for = are 



<7D*D.(0) = 22.1 ±0.4 



5i?*B.(0) = 60.5 ±1.1 



(3.53) 



to be compared to the experimental determination gn* £)^(m^) = 17.9±0.3±1. 9 by 
the CLEO Collaboration [l5i]^e lattice results 90*0^(^1) = 18.8±2.3l2:o 
and gB'Biri^) = 47 ± 5 ± 8 [156l |. or a recent determinati on u sing QCDSR for 



157^ . Older QCDSR 



which 90*0^^1) = 14.0 ± 1.5 and gs-'BTriO) = 42.5 ± 2.6 
results give smaller values for both coupling constants. For instance, the calcu- 
lation within QCDSR on the light cone in Ref. ^158i | give go'DTrim'^) = 12.5 ± 1 
and (7B*B7r(0) = 29±3[l|. The latter are small compared to lattice data or the 
experimental determination of gD*D-K{'m'^) by the CLEO Collaboration. 
From our results we obtain the ratio 



R 



gB*Bn{0) Ib* 



1.105 ±0.005 



(3.54) 



in good agreement with HQS that p redic ts a value of 1 with l/mc^b corrections 
appearing in next to leading order [l33l | Our result in Eq. (j3.54p is also in 
agre emen t with the one obtaine d com bining lattice data for /b* and /d* from 
Ref. |ll8l |. for gB* b-k{0) from Re f. |l55l] . and the experimental CLEO Collaboration 
data for g^)* 



from Ref. 1541 ]. In this case one gets 



xp.—Latt. 



1.26 ±0.36 



(3.55) 



where we h ave added errors in quadratures. A calculation using light cone QCDSR 
gives 



1581 ] 



R\ 



QCDSR 



0.92 



(3.56) 



^ Values for b oth coupling constants obtained prior to 1995 within different approaches can 
be found in Ref. [l58{ and references therein. 

*Note the strong coupling constant used in Ref. [13S ] is given in terms of ours as 
{gH'H^U)/{2mH') with H = B,D. 



Chapter 4 



Semileptonic and non-leptonic 
decays of the B~ meson 



4.1 Introduction 



Since its discovery at Fermilab by the CDF Collaboration 159 . 16Cl| | the Be meson 



has drawn a lot of attention. Unlike other heavy mesons it is composed of two 
heavy quarks of different flavor (6, c) and, being below the B-D threshold, it 
can only decay through weak interactions making an ideal system to study weak 
decays of heavy quarks. 



Using HQSS Jenkins et al. |72[ | were able to obtain, in the infinite heavy quark 
mass limit, relations between different form factors for semileptonic Be decays into 
pseudoscalar and vector mesons. Contrary to the heavy-light meson case where 
standard HQS applies, no determination of corrections in inverse powers of the 
heavy quark masses has been worked out in this case. So one can only test any 
model calculation against HQSS predictions in the infinite heavy quark mass limit. 

With both quarks being heavy, a nonrelativistic treatment of the Be meson 
should provide reliable results. Besides a nonrelativistic model will comply with 
the constraints imposed by HQSS as the spin-spin interaction vanish in the infin- 
ity heavy quark mass limit. In this chapter we will study, within the framework 
of a nonrelativistic quark model, exclusive semileptonic and nonleptonic decays 
of the B~ meson driven by a 6 ^ c or c — > d, s transitions at the quark level. 
We will not consider semileptonic processes driven by the quark b ^ u transition. 
Our ex perie nce with this kind of processes, like the analogous i? — > vr semileptonic 



decay 16ll |. shows that the nonrelativistic model without any improvements un- 
derestimates the decay width for two reasons: first at high q'^ transfers one might 
need to include the exchange of a B* meson, and second the model underestimates 
the form factors at low q'^ or high three-momentum transfers. We will concentrate 
thus on semileptonic B~ — > cc and B^ B transitions. As for two-meson non- 
leptonic decay we will only consider channels with a least a cc or 5 final meson. 
In the first case we will include channels with final D mesons for which there is 
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a contribution coming from an effective b ^ d, s transition. As later explained 
this is not the main contribution to the decay amplitude and besides the momen- 
tum transfer in those cases is neither too high nor too low so that the problems 
mentioned above are avoided. 

The observables studied here have been analyzed before i n th e coiitext o f dif- 
ferent models like the relativistic constituent quark mo del 162l. Il63l . Il64l ]. the 



ativistic QCD [172, ll73|, Il74 



1 to the Bethe-Salpeter equation 



17C 



quasi-potential approach to the relativistic quark model 165l . Il66l| , the instanta- 
neous nonrelativistic app roac 
Bethe-Salpete r equation 



167, 168, 161, the 



17111 . the three point sum rules of QCD an d no nrel- 
175( 1 . the QCD relativistic potential model [176l |. the 



relativistic constituent quark model formulated on the light front 177l |. the rela- 
tivistic quark-meson mod el jl78 l or in models that use the Isg ur, Scora, Grin stein 
and Wise wave functions [82, ISSj like the calculations in Refs. HHllsH. We 



will compare our results with those obtained in these latter references whenever 
is possible. Besides, we will perform an exhaustive study for exclusive semilep- 
tonic and nonleptonic B~ decays, paying an special attention to the theoretical 
uncertainties affecting our predictions and providing reliable estimates for all of 
them. 

In the present calculation we shall use physical masses taken from Ref. 1151 ]. 
For the Be meson mass and lifetime we shall use the central values of the recent 
experimental determinations by the CDF Collaboration of m^^ = 6285.7 it 5.3 it 
1.2MeV/c2 [18^ and tb, = (0.463+|]:j^^^ ± 0.036) x IQ-^^g [;l83]. This new mass 
value is very close to the one we obtain with the different quark-quark potentials 
that we use in this memory from where we get = 6291.6^33 MeV. 

We shall also need CKM matrix elements and differ ent m eson decay constants. 
For the former we shall use the ones quoted in Ref. |l62l | that we reproduce in 
Table HTJ All of them are within the ranges quoted by the Particle Data Group 
(PDG) [115]. 



\Vud\ \Vus\ \Vcd\ IK. 



\Vcb\ 



0.975 0.224 0.224 0.974 0.0413 



Table 4.1: Values for Cabibbo-Kobayashi-Maskawa matrix elements used in this 
work. 



For the meson decay constants the values used in this work are compiled in 
Table 14.21 They correspond to central values of experimental measurements or 
lattice determinations. The results for fp and fx* have been obtained by the 
authors in Ref. [162] using r lepton decay data. Our own theoretical calcula- 
tion, obtained with the model described chapter O give fp = 0.189 ~ 0.227 GeV, 
Jk* = 0.180 ~ 0.220 GeV depending on the i nter -quark interaction used, results 
which agree with the determinations in Ref. [l62( ]. We shall nevertheless use the 
latter for our calculations. For /^^ we have been unable to find an experimental 
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termma- 



result or a lattice determinatioij^. The re are at least two theoretical dete 
tions that predict /^^ = 0.484 GeV (mI] and frj, = 0.420 ±0.052 GeV ^8^. Again 
our own calculation gives values in the range frj^ = 0.485 ~ 0.500 GeV depending 
on the inter-quark interaction used. Here we will take frj^ = 0.490 GeV. 



fK-,KO fK*-,K*o 



0.1307 [115] 0.130 [115] 0.210 [162] 0.1598 [115] 0.217 [162] 

fvc fj/^ 

0.490 0.405 [186] 



Td- 



D* 



0.2226 [123] 0.245 [187] 0.294 [188] 0.272 [187 



Table 4.2: Meson decay constants in GeV used in this work. 



4.2 Meson states 



The meson states were introduced earlier in chapter [3l Here we will need the 
ground state wave function for scalar (0"^), pseudoscalar (0~), vector (1~), axial 
vector (l"*"), tensor (2"*") and pseudotensor (2~) mesons. Assuming always the 
lowest possible value for the orbital angular momentum we will have for a meson 
M with scalar, pseudoscalar and vector quantum numbers: 



^ There is a determination by the CLEO Colaboration 'l84\ using factorization aproximation 
and thus model-dependent. 
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1 1 



E 



2' 2' 



0-1,012 



— 5 (^(°") 



si,S2,-m) (1,1,0| m, -m,0) Yi„i(p ) 



1 



;^'^C1,C2 '^£,A),(52,/2)(^) 



1 1 

2'2^ 



2^(P) 



,0 



si,S2,0 ioo(p ) 



si,S2,A loo(p) 



(4.1) 



For axial mesons we need orbital angular momentum L = 1. In this case two 
values of the total quark-antiquark spin Sqq = 0, 1 are possible, giving rise to the 
two states: 



,/2) 



1 



-'CI , C2 



■1) ^uh''~'^(\P\) 



-^5,,,,, (-1) (/>);y^ 
1 1 



(l+,5„-=l) 

2 



(IpI) 



E 



2'2' 



1 1 
2'2' 



si,S2,0 Yix{p ) 



si,S2,A-m) (l,l,l|m, A-m, A)Yim(p ) 

(4.2) 
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Finally for tensor and pseudotensor mesons we have the wave functions: 
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Cl , C2 




i\p\) 




T 2 



1 1 



-, 1 si, S2, A — m 



) 



(l,l,2|m,A-m, A) Y\ni{v) 





Cl , C2 



-,-,1 si,S2,X-m 



(-1)<A 



) 



(2,l,2|m,A-m,A)y2m(p) (4.3) 



All phases have been introduced for later convenience, and radial wave function 
are computed using the potentials described in chapter [2j The use of different 
inter-quark interactions will provide us with a spread in the results that we will 
consider, and quote, as a theoretical error added to the value obtained with the 
ALl potential, that we will use to get our central results. 



4.3 Semileptonic cc decays 

In this section we will consider the semileptonic decay of the B~ meson into 
different cc states with 0~, 1"*", 1~, 2+ and 2~ spin-parity quantum numbers. 
Those decays correspond to a 6 — > c transition at the quark level which is governed 
by the current 



(4.4) 



with a quark field of a definite flavor /. 
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4.3.1 Form factor decomposition of hadronic matrix elements 

The hadronic matrix elements involved in these processes can be parametrized in 
terms of a few form factors as 



cc{l-),\Pcc Jf(0) 



Be, Pb, 

Be, Pb, 



cc(0" 



), Pec 



Be, Pb. 



P^F+{q^) + q^F^{q^ 



cc(l"), AP„ 
-1 



rriB^ + me 



~-fivaP c^(A) V ^ce 



Pre)P 



•a r,P 



Be, Pb. 



i { ("T-Bc - rnee) £*X)^{ Pce ) Ao{q ) 



rriB^ + me 



{P^A^{q^)+q,A^{q^)) 



Cc{2+),\Pe-e 4\0) 



Be, Pb. 



cc(2+), AP,5 J^'^(O) - Jr^(O) 
e^uap ef/)* ( Pee )P5 P° Q " T^{q^) 
^< ^Uf^siPc-c)P'Ti{q') 

+P''P'el^)^,{Pe-e, ) {P, T2{q^) + q, n{q^)) | 

(4.5) 



Be, Pb, 



In the above expressions P = Pbc+Pcc, Q = Pbc~Pcc, being Pb^ and Pec the meson 
four-momenta, uib^ and niec are the meson masses and e(A)^(-P) and £[X)nv{PX 
are the polarization vector and tensor of vector and tensor mesons respectivelyo 
The latter can be evaluated in terms of the former as 



e'^^){P) = Y.(l,l,2\m,\-m,\) ( ^ ^A-™) ( ^ ) 

m 

Besides the meson states in the Lorenz decompositions of Eq. 
malized such that 



M, A' P' I M, A P ) = A i2^f 2 Em{P)5{P' - P ] 



(4.6) 

5]) are nor- 

(4.7) 



Note the factor 2Em of difference with Eq. ()3.6p 

For the 0+, 1+ and 2" cases the form factor decomposition is the same as for 
the 0~, 1^ and 2+ cases respectively, but with — -/^^(O) contributing where </y^(0) 
contributed before and vice versa. 



Note we have taken A to be the third component of the meson spin measured in the meson 
center of mass. 
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The different form factors in Eq. (j4.5p are all relatively real thanks to time- 
reversal invariance. , V, Aq, A^, A- and Ti are dimensionless, whereas 
T2, T3 and T4 have dimension of E^'^. They can be easily evaluated working in 
the center of mass of the Be meson and taking q in the positive z direction, so 
that Pec = —q = —\q\k. 



Be ^rjcl vu Xd)l i^i decays 



Let us start with the B^ decays into pseudoscalar rjc and scalar XcO cc mesons. 
For B~ — > r/c transitions the form factors are given by: 



whereas for Be 

F-{q') 



1 



2mB^ 
1 



V'{\q\) + 



{Fr,S-q) + mB, 



(4i 



XcO 



transitions we have: 



-1 



-1 
2mD 



A\\q\) + ^^{EU-q. 



A\\q\) + 



191 
A\\q\) 



(4.9) 



with Vf'ilql 



and A^(l^l) (^ = 0, 3) calculated in our model as 



Vc, 



k 



Jy "(0) 



B-,0 



2mBjEr,^i-q) 



XcO, 



J 



cb /I 



NR 

(0) 



\q \ k 



Jy 



Be, 0) 

/ NR 



Be,0 



2mBjE. 



XcO\ 



NR 



XcO, -\q\ k 



Jl^^(O) i?-,0)^^ (4.10) 



/ NR 



which expressions are given in appendix lEl 

In Fig. 14.11 we show our results for the -F+ and F_ form factors for the semilep- 
tonic Be Tje, XcO transitions. The minimum ((^ value depends on the actual 
final lepton and it is given, neglecting neutrino masses, by the lepton mass as 
9m,m ~ "''F- "^^^ form factors have been evaluated using the ALl potential. For 
decays into r/c, and for the sake of comparison, we also show the results obtained 
using the BHAD potential. As seen in the figures the differences between the form 
factors evaluated with the two inter-quark interactions are smaller than 10%. 

In Table 14.31 we show and F_ evaluated at g^j^ and q^^^ for a final light 
lepto n (I = e, /i) and compare them to the ones obtained by Ivanov et al. in 
Ref. [163], and, when available, by Ebert et al. in Ref. 165]. For the Be rje 
transition we also show the corresponding values for the Fq form factor defined as 

Fo{q') = F+iq') + , , F.iq') (4.11) 



m 



mt 
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1,5 r 




I ^ L 

4 6 
[ GeV^ ] 



ALl 



Figure 4.1: and F_ form factors for B~ — 77c and B~ — > XcO semileptonic 
decays evaluated with the ALl potential. In the first case, and for comparison, 
we also show with dotted lines the results obtained with the Bhaduri (BHAD) 
potential. 



Our results for the rjc case are in excellent agreement with the ones obtained by 
Ebert et al. Compared to the results by Ivanov et al. we find large discrepancies 
for 



• B: 



J/'^lvu hcli^i, Xcili'i decays 



Let us now see the form factors for the semileptonic B^ decays into vector J/^ 
and axial vector he {Sgg = 0) and Xci {Sqq = 1) cc mesons. For the decay into 
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^min 


"imax 


Be XcO I T^l 




"imax 


F+ 
This work 

[163] 
[165] 


0.49+0-01 
0.61 
0.47 


1.00_o.oi 
1.14 
1.07 


F+ 
This work 
[163] 


0.30+0-01 
0.40 


64+0-01 

'-'•"^-0.01 

0.65 


F_ 
This work 
[163] 


-0.11+°°i 

—0.01 

-0.32 


-0.25+003 
-0.61 


F_ 
This work 
[163] 


-0.62+0-01 

—0.03 

-1.00 


-1.35+0-05 
-1.63 


This work 
[165] 


0.49+0.01 
0.47 


0.91+0.01 
0.92 









Table 4.3: an d F evaluated at q^^^ and qr^^^ compared to the ones obtained 
by Ivanov et al. jl63l ] and Ebert et al. Ref. [165]. Our central values have been 
obtained with the ALl potential. For the r]c channel we also show Fq (see text for 
definition). Here / stands for I = e, fi. 



J/^ the form factors can be evaluated in terms of matrix elements as: 



A+{q') 



A-{q' 



Ao{q' 



\/2 ruB, \q\ 



V^a'=-i(1?1) 



2mD 



\q I "^Bc 



-V2 



\q\mj/m 



^Uid^D 



'imB^ 



|9| "^Bc 



AUm) + 



mB,+ Ej/^,{-q) 



A 



A=Ol 



-V2 



mB,Eji^{-q)+m]j^ 



\q\mj/si, 



iV2- 



1 



-4L-i(k1) 



(4.12) 
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with and |) calculated in our model as 



J/^, A -\q\k 



jcb fj, 



(0) 



B-, 



J/^, A -\q\k 



J 



ch fi 



(0) 



B-, 



2mB^2Ej/^{-q) ( J/^, A 



Jy "(0) 



4' "(0) 



B-, 6) 

/ NR 



B-,6) 

/ NR 



(4.13) 



which expressions are given in appendix lEl 

The form factors corresponding to transitions to the Xci and he axial vector 
mesons are obtained from the expressions in Eg. (14. 121) by just changing 



Vti\q\) 



(4.14) 



and using the appropriate mass for the final meson. Obviously in Eq. ()4.13p J/^ 
has to be replaced by Xci or he- 

In Table [13] we show the results for the different form factors evaluated at q'^^^ 
and (/max *2^se where the final lepton is light {I = e, fi). For the decay into 

J/'f we also show the combination of form factor^: 



A-{q^) (4.15) 



Our results for the Be J/^ decay channel are in agreement with the ones 
obtained by Ebert et al. They also agree reasonably well, with the exception 
of A-, with the ones obtained by Ivanov et al. For the other two cases the 
discrepancies are in general large. 

All the form factors are depicted in Figs. 14.21 and 14.31 



•B^ ^ ^'(3836)/z^i, Xc2li^i, decays 

Finally let us see the form factors for the B~ decays into tensor Xc2 and pseudoten- 
sor ^(3836) mesons. For the decay into Xc2 the form factors can be evaluated 



■^This combination is called by the authors of Ref. [l65l |. 

''Note that while the ^^(SSSG) was still quoted in the particle listings of the former Review of 
Particle Physics [llSl ]. it has been excluded from the more recent one 188]. We shall nevertheless 
keep it in our study to illustrate the results to be expected for a ground state pseudotensor 
particle. 
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B- ^ J/^l-Di 


^min 


(7^ 


B^ ^ hcl vi 


^min 


g2 

^max 


V 

This work 
[163] 
[1651 


-0.61-0.03 
-0.83* 
-0.49 


-1.26+001 
-1.53* 
-1.34 


V 

This work 
[163] 


-0.040 - 0.003 
-0.25 * 


— 0.0 78 - 0.003 

-0.365* 


A+ 
This work 
[163] 
[165] 


0.56+003 
0.54 
0.73 


1.13+001 
0.97 
1.33 


A+ 
This work 
[163] 


-0 8f^+0-"l 

-1.08 


-1.90+0 06 
-1.80 


A_ 
This work 
[163] 


-0.60_o.03 
-0.95 


-1.24_o.oi 
-1.76 


^- 
This work 
[163] 


'-'■^^-0.02 

0.52 


0.36-0.06 
0.89 


This work 
[163] 
[165] 


1.44+0-08 

1.64 
1.47 


9 ro+0.01 

2.50 
2.59 


^0 

This work 
[163] 


0.28+0 02 
0.44 


52+0.02 
0.54 


Ao 
This work 
[1651 


0.45+0.03 
0.40 


0.96-0.01 
1.06 








Be ^ Xci 


^min 


^max 








V 

This work 
[163] 


'-'■•^^-0.02 

1.18* 


1.86-0.12 
1.81* 








A+ 
This work 
[163] 


-0.44-0.03 
-0.39 


-0.78+0-04 
-0.50 








A_ 
This work 
[163] 


96+° °'' 
1.52 


1.97-0.13 
2.36 








This work 
[163] 


-0.50-0.02 
-0.064 


—0.32-0.02 
0.46 









Table 4.4: Aj^, A- and Aq form factors evaluated at g^-^^ a nd qL ^^ compared 
to the ones obtained by Ivanov et al. [163] and Ebert et al. Ref. |l65l |. Our central 
values have been evaluated with the ALl potential. For the J/'^ channel we also 
show Aq (see text for definition). Here I stands for / = e, /i. The asterisk to the 
right of a number mea ns w e have changed its sign to account for the different 
choice of e°^23 j^g£_ 
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Figure 4.2: V (solid line), A^ (dashed line), A_ (dotted line) and (dashed- 
dotted line) form factors for the B~ J/'^ semileptonic decays evaluated with 
the ALl potential. 



in terms of matrix elements as: 



1 f TnP P\ Tn?' 

ue) = ^i^[-^2^^Tx=om)-^i^^{E^A-^)-mB:)A^,^m) 



= ^;;f]|2^TA=+i(k1) (4.16) 
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Figure 4.3: V (solid line), yl+ (dashed line), (dotted line) and Aq (dashed- 
dotted line) form factors for B~ he and B~ Xci semileptonic decays evalu- 
ated with the ALl potential. 

with Vjf^(|g I) and Aj,_^(|g |) calculated in our model as 
= {xc2, X - \q\k\4"'{0)\B-,0) 



Be, O) 



A^xiWl) = { Xc2, A - \q\ k I '^(O) I B-, 



/ NR 



= ^2mBjE^A-<l) ^^{Xc2, X-\q\k\rJ"'{0)\B-,6 



NR 

(4.17) 
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which expressions are given in appendix lEl 

The form factors corresponding to transitions to ^'(3836) are obtained from 
the expressions in Eg . (14.161) by just changing 

v^^m ^ -Kxm (4.18) 

and using the appropriate mass for the final meson. Besides in Eq. ()4.17p Xc2 has 
to be replaced by ^'(3836). 

The results for the different form factors appear in Fig. 14. 4[ In Table 14.51 we 
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[GeV" ] 
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T4 [GeV"'] 



ALl 







2 3 4 

q' [ GeV' ] 



Figure 4.4: Ti (bold solid line), T2 (dashed line), T3 (dotted line) and T4 (thin 
solid line) form factors for B~ Xc2 and B~ — > ^(3836) semileptonic decays 
evaluated with the ALl potential. 



4.3. SEMILEPTONIC B~ 



cc DECAYS 



51 



7^ 

imax 



for the case of a final light 
163] • For the 



show Ti, T2, T3 and T4 evaluated at q^^^ and g„ 
lepton, and compare them to the values obtained by Ivanov et al. 
B~ — > Xc2 transition we find a reasonable agreement between the two calculations. 
For B~ ^'(3836) there is also a reasonable agreement for the absolute values of 
the form factors but we disagree for some of the signs. 





2 

'Jmin 


2 

'?max 


B- *(3836)ri5. 


2 

9min 


2 

Qmax 


Ti 

This work 
[163j 


u-y'-o.oi 

1.22 


1.95-0.06 
1.69 


Ti 

This work 
[163] 


0.29-0.02 
0.052 


0.22-0.01 

0.35 


This work 
[163] 


-0.012_o.ooi 
-0.011 


_0 025+oo»i 
-0.018 


Ta [GeV^-^] 
This work 
[163] 


_0.oio+"oo« 
0.0071 


_0,oi8+o.ooi 
0.0090 


T3 [GeV-^] 
This work 
[163] 


OQ^3+o.ooi 
0.025 


0.030-0.001 
0.040 


Ta [GeT/-2] 
This work 
[163] 


0.025-0.002 
-0.036 


0.052-0.004 
-0.052 


Ta [GeV"-^] 
This work 
[163] 


00^3+0.001 
0.021* 


0.030-0.001 
0.033* 


r4 [GeV-2] 
This work 
[163] 


0.022-0. 001 
-0.026* 


0.045 - 0.003 
-0.038* 



Table 4.5: Values for Ti, T2, T3 ar id T 4 evaluated at g^j^ and 
the ones obtained by Ivanov et al. 1631 ] . Here Z stands for / 
Table! 



jj^a^x compared to 
e, /i. Asterisk as in 



4.3.2 Decay width 

For a i?c at rest the double differential decay width with respect to and x;, being 
xi the cosine of the angle between the final meson momentum and the momentum 
of the final charged lepton measured in the lepton-neutrino center of mass frame 
(CMF), is given bjH 

d'T ^ _3_ X'/\q\ml^,ml^ - mf ^^^^ C^^puP.) 

dq^dxi G'im^ Stt^ '^mB^ q^ o cc i, u 

(4.19) 

where A(a, 6, c) = {a + b — cf' — Aab, mi is the mass of the charged lepton, TLap 
and are the hadron and lepton tensors, and Pbcj Pcs, Pi, Pu are the meson 
and lepton four-momenta. The lepton tensor i^ 

C^^{puPu) = ^(pfp?,+P^iPZ-9''^PrPu^ie''^'"'piaPup) (4.20) 



^We shall neglect neutrino masses in the calculation. 

"The =F signs correspond respectively to decays into l~vi (for B~ decays) and I'^vi (for B^ 
decays) . 
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As for the hadron tensor it is given by 

nap{PB.,Pcc) = /i(A)a(^'B.,Pcc-)/i(A);3(^i?.,^cc-) (4.21) 



with 



hix)a{PB.,Pcc) = {cc, XPM'miB,, Pb^) (4.22) 

The quantity 

n^p{PB.,Pcc) C^^ipuPu) (4.23) 

is a scalar and to evaluate it we can choose Pec along the negative z-axis. This 
implies also that the CMF of the fina l leptons moves in the positive z-direction. 
Furthermore we shall follow Ref. |l63l ] and introduce helicity components for the 
hadron and lepton tensors. For that purpose we rewrite 

nap{PB.,Pcc) C'^iPuPu) = n^'PiPB^^Pcc) 9aa 9pP C'^iPuPu) (4.24) 

and use 189l | 



S'mi' = X] S'rr e{r)/i(9)'^(r)i.(^) ; 5« = 1, 9±i,o = "1 (4.25) 

with £^-j^^{q) = q^/q^ and where the e(j,)(q'), r = ±1, are the polarization vectors 
for an on-shell vector particle with four-momentum q and helicity r. Defining 
helicity components for the hadron and lepton tensors as 

-HrsiPB^^Pcz) = ey)M)'^"'{PB.^Pc^Hs)pi<l) 

CrsiPhPu) = £(r)aiq)^"'^{Pl,Pu)e*,^i3iq) (4.26) 

we have that 

T~(-afs{PBc,Pcc) ^"^{PhPu) = ^ ^ Qrr Qss'HrsiPB,, Pec) Crs{Pl,Pu) 

r=t,±l,0 s=i,±l,0 

(4.27) 

Let us start with the lepton tensor. We can take advantage of the fact that the 
Wigner rotation relating the original frame and the CMF of the final leptons is the 
identity to evaluate the lepton tensor helicity components in this latter reference 
system 

Crs{puPu) = e(^r)a{<l)C-°''^{PuPu)e\,)p{q) = e^r)aiq)^"''^iPi,Pu)els)(3iQ) 

(4.28) 
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were the tilde stands for momenta measured in the final leptons CMF. For the 
purpose of evaluation we can us€0 

= {Eii\pi\),-\pi\ -\pi\xi) 

= i\pi\,\pi\^l-xlO,\pi\xi) (4.29) 

with \pi\ the modulus of the lepton three- momentum measured in the leptons 
CMF. 

The only helicity components that we shall need are the following. 





,Pu) 












<Py) 


= ^Ot{Pl,Pu] 


) = -Axi , 
q^ 


£+1+1 (pz. 


<Py) 


= il'-mf) 


(^4(1 ± xz) - 2(1 




,Pu) 


= {i-mf) 


(^4(1^x0-2(1 


^OoiPl: 


<Pu) 


= 4{q'^-mf 


)(i-.r'7^ 



(4.30) 

As for the hadron tensor it is convenient to introduce helicity amplitudes defined 

as 

V),(Pb.,Pcc-) =4),(g) /ifA)(^i?c,^'cc) , r = i, ±1, (4.31) 
in terms of which 

-HrsiPB^.Pcc) = ^h(^X)r{PB.,Pcc)hl^^^iPB^,P,,) (4.32) 
A 

We now give the expressions for the helicity amplitudes evaluated in the original 
frame. 

• Case 0- ^ 0~,0+. 

2 2 

/it(PB.,Pcc-) = '^^ F+(g^) + v^F_(g^) 



h^{PB.,Pc^ = A' F+iQ 



/i+i(Pb,,Pcc-) = h^i{PB^,Pc-c) = (4.33) 



^Note this is in accordance with the definition of xi and the fact that we have taken Pec in 
the negative z direction. Furthermore there can be no dependence on the ipi azimuthal angle so 
that we can take Pj, and then p^, in the OXZ plane. 
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• Case 0- ^ 



h(\)tiPBc,Pc< 



h^X)+liPB.,Pcc) 
hiX)-l{PB.,Pcc) 
h(X)oiPBc,Pcc) 



i6 



AO 



-A-{q' 



Jl ^2 



rrir. 



y(g^) + {ttib, - nice) ^o(g^ 



-i6 



A+l 



/ A^/2(r,mg ^,mgg; ^^^^^^ ^ _ ^ A^f^2 

2 „2 ^2 



-^^ y{Q ) + ("ific - iT^cc) ^o(r ) 

"T-Bc + "ice 



/ — Q nice 2 

^'^AO ("T-Bc - "ice) ^ -F= Ao{q 

\ 2niccyq 

Jl ^2 ^2 



\{q\ni%^,mic) A+{q^) 



IniccJq^ "iBc + "ic 



(4.34) 



Case 



2-,2^ 



h{X)t{PB,,Pa 



\x)+i{Pb,,Pcc) 



h(X)oiPB,,Pcc) 



.^A(g2,m| m; 
-^^aoa/tt 



4m2g yV 



/1(A)-i(Pb,,^'cc) = ^f^A+l^ 



( ri(g2) + - m^,) T2{q^) + r3(<?2) ) 

1 X^/Hq^,ml^,mlc) 

'^'-'T2 

Tr[q^)-X^'\q\nil^,nilc)T^{q^) 
1 AV2(g2^m| ,m2,-) 



AO 



^/2 2 mcc 

ri(<z') + A^/2(^2^^|^^^2_)j.^(^2^ 
2 Al/2(g2^^| 2_) 



3 



4m2g 

2 2 

9 -"i« 



)ri(g2) + A(g2,m|^,m2,)r2(g2)) 

(4.35) 



We see that the hehcity ampUtudes, and thus the hehcity components of the 
hadron tensor, only depend on q^. The expressions for the latter are collected in 
appendix [F1 
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We can now define the combinations |l63l ] 





— "W+i +1 + W-i -^1 


Hl 


= Woo 


Hp 




Hs 




HSL 


= T~ito 






Hj 


= TTl^j ; J 
2 



J = U,L, S, SL (4.36) 

with U, L, P, S, SL representing respectively unpolarized-transverse, longitudi- 
nal, parity-odd, scalar and scalar-longitudinal interference. 

Finally the double differential decay width is written in terms of the above 
defined combinations as 



\Vb. 



^(1 + xj) Hu + \{1- xf) Hl ± \xiHp 

+|(1 - x'i)Hu + ^x'iHL + \hs + 3x,#5l} (4.37) 

Note that for antiparticle decay Hp has the opposite sign to the case of particle 
decay while all other hadron tensor helicity components combinations defined in 
Eq. (j4.36p do not change (See appendix IE. 21 for details). The sign change of Hp 
compensates the extra sign coming from the lepton tensor. This means that in 
fact the double differential decay width is the same for or B'^ decay. 
Integrating over xi we obtain the differential decay width 

- 8^'^''^' T2<? ^Hu + Hl + Hu + Hl + Hs 

(4.38) 

from where, i ntegr ating over g^, we obtain the total decay width that we write, 
following Ref. 163], as 



r = Tu + Tl + Tu + Tl + Ts (4.39) 
with Tj and Tj partial helicity widths defined as 

2 Gl , ,2 (g2 - m' f Al/2(g2 I ^2 ) 

10 2 2 o 4.40 

and similarly for Tj in terms of Hj. 
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rsL 



Vc e u, 



Xcoe I^e 
XcQT-Ur 

Xcl e'i^e 

Xcl t^'^lJ. 
Xcl T-Ur 

he e~Pe 

he T~ Ur 

Xc2 e~Pe 
Xc2 M'i^M 
Xc2 t~Pt 

■^(3836) e-pe 
'I'(3836) AJ-p^ 

*(3836) T~Vr 



11.5+0-6 

114+0.6 

r, 70+0.10 

90+°-''''^ 
'^■^'-'-0.03 

89+°-"^ 
0.75+0'O2 10-1 
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0.23 



f-0.02 



10- 



, +0.03 
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0.71] 

0.71^ 
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0.13+0-01 10-5 
28+0-02 10-1 

0.17+0-01 



1.55 



+0.14 
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-0.13 



1-51-0.02 
0.80+nn^ 10" 



0.37+0.05 
0.75+0.09 

0.23 



+0.01 



10" 
10" 
10" 



0.32+0-02 10-6 
0.13+0-01 10-1 
0.59+002 
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10.2+0-^ 

1 74+0.07 
'^-0.01 



0.12+0-01 10-5 
0.28+0.03 10-1 
0.39+0.02 



0.43l°;03 10 
0.18+00? 10 



0.35+0.03 10-1 0.28+0.03 10 



-0.01 
0.21+0.01 10-1 

0.57+0.07 10-8 
0.24+0.03 10-3 



0.35+0.03 10-1 
0.46+0.04 10-2 

2.23+0.12 

2 lfi+0-21 
^•10-0. 01 



0.60+0.06 10-2 0.67+0.05 10-1 



0.77+0.08 10-4 
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0.72+0.09 10-6 

0.14+0.01 10-1 
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-0.27- 0.03 10" 



1-0.02 
-0.02 

0.15+0-01 



0.37'^ 



0.15- 



0.01 



1.17 
1.14 
0.43 



0.08 
0.05 
+0.07 
0.05 
+0.01 
0.02 



0.31 
0.62 



0.03 
0.01 
+0.06 
0.02 



10-6 

10-1 0.13-0.01 10-1 -0.18+0-01 10- 



-0 -^5+0-01 

-0 34+0-01 
'-'-■'^-0.02 



44+0.03 iQ-E 

0.10+0.01 
1.58+0-04 



0.14+0.01 10-5 
0.31+0-01 10-1 
0.29+0.01 



0.11+0.01 10-5 0.37+0 05 10-6 

0.23+0-02 10-1 0.75+0-09 10-2 

0.84+0.07 10-1 0.25+0-02 10-1 

0.32+0.03 10-5 0.11+0.01 10-5 
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-8 n qci+0.01 
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-3 n Qo+o.oi 

'^•'^^-0.02 

-3 0.74 



10 
10 

0.06 10 



10 



'2 n fi8+0.04 
0.b8_o 04 

-2 0.15-0.01 10-'' 

-4 0.25-0.02 10-4 



0.60i0.o3 10- 



-0.48+0.05 10-1 
-0.48+0.06 10-1 

-0.69+0.08 10-3 o.54-ao4 10-5 0.65-o.o5 10 



O.lTlHi 10- 
'•281^2 10" 
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Table 4.6: Partial helicity widths in units of 10 GeV for B^. decays. Central 
values have been evaluated with the ALl potential. 

Another quantity of interest is the forward-backward asymmetry of the charged 
lepton measured in the leptons CMF. This asymmetry is defined a^ 

Afb = ^"'^°7^"'<° (4.41) 
and it is given in terms or partial helicity widths as 

3 ±rp + 4r5L 



4 Tu + Tl + Tu + Tl+T. 



(4.42) 



being the same for a negative charged lepton {B~ decay) than for a positive 
charged one {B^ decay), as Fp for antiparticle decay has the opposite sign as 
for particle decay. 

Finally for the decay channel Be — > J/^ Wi with the J/^ decaying into iJi^ iJi^ 
we can evaluate the differential cross section 



The forward direction is determined by the momentum of the final meson that we have chosen 
in the negative 2-direction. 
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dT 



dx^ 
1 + - 

r 



1 - 4m2 /m2/^ [Tu + Tu- 2{Tl + + 



2 - Jl - Am? /m? 



1 



r 



1 + 2m2 /m2^^ 
4 (2 



1 - 4ml/m]^^ 



u 



(4.43) 



where is the cosine of the polar angle for the final lJi~ jJi^ pair, relative to the 
momentum of the decaying J/^*, measured in the n~ ii^ CMF, T is the 
J/^' decay width into the fj.~ ^'^ channel, and Tj/q, is the total J/^' decay width. 
The asymmetry parameter 



a 



2- Jl-Aml/m']^^ 



Tu + Tu]+2[Tl + Tl + Ts 



(4.44) 



governs the muons angular distribution in their CMF. 
Results 

In Table 14.61 we give our results for the partial helicity widths corresponding to 
B~ decays. For B'^ decay the "P" column changes sign while all others remain 
the same. The central values have been evaluated with the ALl potential and the 
theoretical errors quoted reflect the dependence of the results on the inter-quark 
potential. 

In Table we show the asymmetr y pa rameters. Our values for a* compare 



well with the results obtained in Ref. [1631 ]. The same is true for the forward 



backward asymmetry with some exceptions: most notably we get opposite signs 
for Be Xci and B^ 1'(3836). 

In Fig. 14.51 we show the differential decay width dV /dq^ for the decay channels 
rjcl^vi and XcqI~^i for the case where the final lepton is a light one 
or a heavy one I = t. We show the results obtained with the ALl and BHAD 
potentials, finding no significant difference for the r case, while for the light final 
lepton case the differences are around 10% at low q^. 

In Fig. 14.61 we show now the results for vector and tensor mesons. As before 
only for the case where the final lepton is light we see up to 10% differences 
between the calculation with the ALl and the BHAD potentials. 



We show the distribution corresponding to a final electron. The distribution for a final muon 
differs from the former only for around mi. 
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Figure 4.5: Differential decay width for the B~ — > r}cl~i'i and B~ Xcol~^i 
processes obtained with the ALl potential. For comparison, we also show with 
dotted lines the results obtained with the Bhaduri (BHAD) potential. The dis- 
tribution for a final muon is not explicitly shown. It differs appreciably from the 
corresponding to a final electron only for around m? . 
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Figure 4.6: Differential decay width for the J/'^l~i?i, B~ hcl~i>u B~ — 

Xcil~^li B~ Xc2i'~^l ^-iid B~ ^(3836)/^^/ decay channels obtained with the 
ALl potential (solid lines) and the Bhaduri (BHAD) potential (dotted lines) . The 
distribution for a final muon is not explicitly shown. It differs appreciably from 
the corresponding to a final electron only for around mi. 
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Table 4.7: Asymmetry parameters in semileptonic Be — > cc decays. Our central 
values have been evaluat ed w ith the ALl potential. We also show the results 



obtained by Ivanov et al. 163l | 



Finally in Tables 14.81 14.91 we give the total decay widths and corresponding 
branching rat ios fo r the different transitions. The branching ratios evaluated by 
Ivanov et al. |l62l ]. where they have used the new Be mass determination by the 
CDF CoUaboration [l82], are in reasonable agreement with our results. Discrep- 
ancies are larger for the decay widths in Tab le 14.81 where they use the larger mass 
value mBe = 6400 MeV quoted by the PDG [ll5l |. 



4.3.3 Heavy quark spin symmetry 

As mentioned in the introduction one can not apply HQS to systems with two 
heavy quarks due to flavor symmetry breaking by the kinetic energy terms. The 
symmetry that survives for such systems is HQSS amounting to the decoupling 
of the two heavy quark spins. Using HQSS Jenkins et al. [72] obtained relations 
between different form factors for semileptonic Be decays into ground state vector 
and pseudoscalar mesons. Let us check the agreement of our calculations with their 
results. For that purpose let us re-write the general form factor decompositions 
in Eq. (j4.5|) introducing the four vectors v and k such that 



PBe = mBe V ; Pcc = rricc v + k 



(4.45) 
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Table 4.8: Decay widths in units of 10~ GeV for semileptonic B~ — > cc decays. 
Our central values have been evaluated with the ALl potential. Here I stands for 
I = e, fi. 



V is the four-velocity of the initial Be meson whereas A; is a residual momentum. 
In terms of those we have 

P = Pb^+ Pcc = {rriB, + nice) v + k (4.46) 
q = Pb^ - Pcc = {ruB^ - nice) v - k (4.47) 

and we can write for the r]c (cc(0~)) final state case 



j;i'{o) 



Be, Pb. 



Be, Pb. 

{{mB, + mr,J F+{q'^) + (ttib, -m^jF_(g2)) 
+ {F+{q')-F^iq'))k, 



'{q')v, + j:), '{q')kj{4.48) 



where we have introduced the new form factors 
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Table 4.9: Branching ratios in % for semileptonic — > cc decays. Our central 
values have been evaluated with the ALl potential. Here I stands for I = e, 11. 
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y/2mij^2m^^ 
1 



(4.49) 
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Similarly for the J/'^ {cc{l )) final state case we have 



^^^^ e^.ap sr^JPj/^)v^kPV{q') 



-^1 ("T-Bc -mj/^) el^)^{Pj/^,)Ao{q'^) 

■ ( Pj/^ ) ( {{mB^ + mj/^) v, + k,) 

+ (("T-Bc - mj/isf) -kfj) \\ 

- ^2mBjmj/^ e^,ap e(A) ( ^J/* ) ^ id') 



+ ik-el,^{Pj/^)) ^P{q') 



+ 



[k-el^{Pj/^)) ^f\q')k,~^ 



(4.50) 



with 



^^'-\q') = \ {ruB^-mj/^) A^{q^) 



^2 



(5^) = 



(5^) = 



(9^) = 



1 -2mB, 



g2 + ^ ^ A_iq^ 

TUB, + nij/^ 

— ^ — 

^2mBj2mji^ m-B^ + mj/^ 

\, ^ 1 {A^{q') - A.{q^)) (4.51) 

^2mB^2mj/i^ mj/^ ruB^ + mj/^ 

T,^ "* and ^ arc dimcnsionlcss, £2*^ \ ^ have dimensions of E~^, 

rrib S> rric 3> Aqqd with the 



and has dimensions of E^^. 

We can take the infinite heavy quark mass limit 
result that near zero recoil 



= 



= 



(4.52) 
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Figure 4.7: E^"" ^ (solid line) and ^ (dashed line) of the rjc, and T,\ ' 

(solid line), T,^ ^ (dashed line), ^,'2^ ^ (dotted line) and T,^ ^ (dashed-dotted 
line) of the B~ J/^ semileptonic decays evaluated with the ALl potential. 
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This agrees perfectly with the result obtained in Ref. [tI] using HQSd^. 

In Fig. 14.71 we give our results for the above quantities for the semileptonic 
r]c and B~ — > J/^ decays for the actual heavy quark masses. Even though 

we are not in the infinite heavy quark mass limit we find that Tif* ^ and ^ 
dominate over the whole interval. This dominant behavior would be more so 
near the zero-recoil point where A; ~ and thus the contributions from the terms 



-(0- 

-'2 



and 



are even more suppressed. Thus, even for the 



actual heavy quark masses we find that near zero recoil 



J/* 



B-, Pb. 



(4.53) 



Besides, as seen in In Fig. 14. 8| ^ of the B~ — > r/c, and ^ of the B~ — > J/^ 
semileptonic decays are very close to each other over the whole interval. This 
implies that the result obtained in Ref. [73] near zero recoil using HQSS seems 
to be valid, to a very good approximation, outside the infinite heavy quark mass 
limit. 



Ir 

0,9 
0,8 
0,7- 



- £/' V)(B"^J/»P) 



ALl 



4 6 
[GeV^] 



10 



Figure 4. 



(solid line) of the B^ Vc, and ■* (dashed line) of the 



B^ J/^ semileptonic decays evaluated with the ALl potential. 



Note, however, the different global phases and notation used in Ref. [t^I- 



66 



CHAPTER 4. DECAYS OF THE MESON 



4.4 Nonleptonic cc Mp two— meson decays 



In this section we will evaluate decay widths for nonleptonic B~ — > cc Mp two- 
meson decays where Mp is a pseudoscalar or vector meson. These decay modes 
involve a 6 — > c transition at the quark level an d th e y are gov erned, neglecting 
penguin operators, by the effective Hamiltonian 162l . Il65l . Il76l | 

where ci, C2 are scale-dependent Wilson coefficients, and Qf', Q'2 are local four- 
quark operators given by 



H,ff, = ^ {V^ [ci(^) Qf + C2(/u) Qf] + H.c.} 



(4.54) 



cb 



^e(0)7MU-75)*fe(0) 



^dioh'^ii - 75)^«(o) + K*s ^^(0)7^^(1 - 75)^.(0) 



+V:^ ^rf(0)7^(/ - 75)^c(0) + ^',(0)7^(1 - 75)*c(0) 



Q 



cb _ ,T, 



^rf(0)7^(/- 75)^-6(0) 
+^,(0)7^(1 -75)^fe(0) 



v:,^c{o)rii - 75)^n(o) + v:,^Mrii - 75)^c(o) 

Ks^c(0)7^(/ - 75)^n(0) + ^e(0)7^(/ " 75)^c(0) 



(4.55) 



where the different Vjk are CKM matrix elements. 

We shall work in the factorization approximation which amounts to evaluate 
the hadron matrix elements of the effective Hamiltonian as a product of quark- 
current matrix elements: one of these is the matrix element for the B~ transition 
to one of the final mesons, while the other matrix element corresponds to the 
transition from the vacuum to the other final meson. The latter is given by 
the corresponding meson decay constant. This factorization approximation is 
schematically represented in Fig. 14.91 




M 



B" 



M 



Figure 4.9: Diagrammatic representation of B^ two-meson decay in the factor- 
ization approximation. 



When writing the factorization amplitude one has to take into account the 
Fierz reordered contribution so that the relevant coefficients are not ci and C2 but 
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the combinations 

aiifJ-) = ci(/i) + ^ C2{fi) ; a2(/i) = C2(/x) + ci(^) (4.56) 
with A/'c = 3 the number of colors. The energy scale u ap propriate in this case is 



/i ~ mi, and the values for oi and 02 that we use are 1621 ] 



ai = 1.14 ; 02 = -0.20 (4.57) 



. Mp =Tr-, p-, K-, K*- 

This is the simplest case. The decay width is given by 



iTA i2 IT/ |2 2 '^^^^('"'Bc' '"cc; 



2 2 2 



(4.58) 

with mp the mass of the Mp final meson, and Vp = Vud or Vf = Vus depending on 
whether Mp = vr", or Mp = K~ , K*~ . Tia/siPscJ Pec) is the hadron tensor for 
the B~ — > cc transition and 'H"^{Pf) is the hadron tensor for the vacuum Mp 
transition. The latter is 

TC'^iPp) = P'^P^fp Mp=0-case 

n^^iPp) = [p^pP -mlg"^)fl Mp = l-case (4.59) 

with fp being the Mp decay constant. 

Similarly to the semileptonic case, the product 'HajsiPBc^ Pec) 'H°^^{Pf) can 
now be easily written in ter ms of helicity amplitudes for the B~ — > cc transition 



so that the width is given as 162l | 



Mp = case 



,2 ^2 ^2 



Mp = 1 case 



^'f it/ |2 it/ |2 2 •^"^'^^("''Bc' "^cc; "^f; 2 ^2 



p _ "-^ IT/ 1^ IT/ K ^-^ V Dc' CO PI y 
^-1^^'^'=^' "1 "^^^^ 

X (w^UI'^'I^^f) + n^Cf^ml) + <o^--''("^f)) (4.60) 

with the different Jirr evaluated at = rr?p . In Table 14.101 we show the decay 
widths for a general value of the Wilson coefficient a\ , whereas in Table 14.111 we 
give the corresponding branching ratios evaluated with oi = 1.14. Our results 
for a fin al r]c or J/ ip are iri good agreement with t he o nes obtained by Ebert et 



al. [165( 1 ■ El-Hady et al. [170] and Anisimov et al. |l77l |. but t hey are a factor 2 
smaller than the results by Ivanov et al. [162] and Kiselev 1741 ] . Large discrepan- 
cies with Ivanov's results show up for the other transitions. 
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r flO~^^ GeVl 

-L 1 J. W \_l Vy V 1 
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Table 4.10: Decay -widths in units of 10~^ GeV, and for general values of the 

Wilson coefficient ai, for exclusive nonleptonic decays of the B~ meson. Our 
central values have been obtained with the ALl potential. 



. M-=D-,D*-,D-,D*- 



In this subsection we shall evaluate the nonleptonic two-meson B^ ricD , 

r]eD*-, J/^D-, J/^D*- andB- VeDj , rjeD*- , J/^D', J/^D*' decay widths. 
In these cases there are two different contributions in the factorization approxi- 
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Table 4.11: Branching ratios in % for exclusive nonleptonic decays of the 
meson. Our central values have been obtained with the ALl potential. 



mation. Following the same steps that lead to Eq. (j4.60p we shall get 



lOvrm^^ 



\Vcb\^\VF{- 



2mBe 



nn 



(4.61) 



where now Vp = Vcd for Mp = D ,D* and Vp 
The quantity TCTl incorporates all information on the hadron matrix elements and 



Vcs for Mp 



D~ D*~ 

^ S ) ^ s 
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depends on the transition as 1621 



-j^-j^B-^VcD*- ^ 



'1 hf" ) m^,- fo- + a2 hf'' {ml J rur,, fr,. 



E 

r=+l,-l,0 



+a2 hfr)7'^* ("^J/*) "^J/* /j/* 



(4.62) 



and similarly for D^, -Dg". Note that the helicity amplitudes corresponding to 
B~ D~ , D*~ , Dj , D*~ have been evaluated from the matrix elements for the 
effective current operators ^'d,s(0)7^(/ — 75)^'fe(0) in Eq. (j4.55p . While in practice 



this is a 6 ^ d, s transition, the momentum transfer (m^2 or rrij^^) is neither 
too high, so that one has to include a B* resonance, nor too low, so as to have 
too high three-momentum transferal. Besides the contribution is weighed by the 
much smaller 02 Wilson coefficient. In Table 14.121 we give the decay widths for 
general values of the Wilson coefficients ai and 02, and in Table W7l3\ we show the 
bra nchi ng ratios. We ar e in reasonable agreement with the results by Ivanov et 
al. [162l |. El-Hady et al. [170l | and Kiselev [174] . For decays with a final _Dr, 
the agreement i s als o reasonable with the results by Colangelo et al. |176t | and 



Anisimov et al. 177 ] 



4.5 Semileptonic B 



, B*^ . B^l, B*^ decays 



B*\ ^° 



B*^ decays. 



In this section we shall study the semileptonic B^ 
With obvious changes the calculations are done as before, with the only novel 
thing that now it is the antiquark that suffers the transition (we have c ^ d, s ), 
and thus we have to take into account the changes in the form factors according 
to the results in appendix IE.2i 



4.5.1 Form factors 

In Fig. 14.101 we show the form factors for the above transitions evaluated with 
the ALl potential. For the B and B^ cases we also show the results obtained 
with the BHAD potential. Although they are less visible in the figures, the larger 
differences, up to 25%, occur for the -F_ form factor. 



^^Note the different phases used in Ref. 
^^Our experience with the B — > vr decay 



162|. 

1611 ]. where we have a similar b u quark transition, 



shows that the naive nonrelativistic quark model gives reliable results for g ~ 9 GeV . 
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r [10-15 GeV] 


B-^rjeD*- 
B-^J/^D- 
B- J/^ D*- 
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(-0.195_o.oo8ai -0.066l°-°??a2)2 
+(-0.390_o.oi8ai - 0.209l°-°J^a2)2 
+(0.447_o.oi6«i + 0.167t[J-°i?a2)' 


B~ rjc 

B-^veDr 

B-^J/^DJ 
B-^J/^D*- 


(2.54+'J-^'5ai + 1.93+'J-iUa2)' 
(-l.84_o.o4ai - 1.17t[!-?2a2)2 
(-1.85l°-0^ai - 1.23_o.06a2)' 
(-l.Ol_o.04a1 - 0.60l0-[!2a2)2 
+(-2.00_o.o6ai - 1.7lt[!-?|a2)' 
+(2.17_o.o8ai + 1.42l0-?2a2)2 



Table 4.12: Decay widths in units of 10^ GeV, and for general values of the Wil- 
son coefficients ai and 02, for exclusive nonleptonic decays of the B~ meson. Our 
central values have been obtained with the ALl potential. For vector-vector final 
state we show the three different contributions corresponding to r = +1, —1, 
(see Eg. droU ). 



In Table 14.141 we show F^, F_ and Fq (defined as in Eq. ()4.1ip changing 
the mass of the final meson) of the B~ B , B^ transitions evaluated at Qmin 
and 0^3. and compare them with the results by Ivanov et al. [l6jj and Ebert et 



al. 166l |. Notice that, to favor comparison, we have changed the signs of the form 



factors by Ebert et al. (they evaluate B'^ decay) in accordance with the results 
in appendix IE.2[ The agreement with the results by Ebert et al. is good. We 
also agree with Ivanov et al. for but get very different results for F-. As 
fermion masses are very small the disagreement in the F_ form factor will have a 
negligible effect on the decay width. 

In Table HT5] we show V , Aq and ^40 (defined as in Eq. ()4.15p changing 

the mass of the final meson) of the B~ — > 'B*^ , B*f eval uated at g^-^ and 



and compare them with the results by Ivanov et al. [I64i | and Ebert et al. |166l ]. 
With some exceptions the agreement with Ebert's results is bad in this case. We 
are also in clear disagreement with Ivanov's results. 



4.5.2 Decay width 

In Tables [3. 161 14.171 we give respectively our results for the partial helicity widths 
and forward- backward asymmetries. 

In Fig. 14.111 we show the differential decay width for the B~ B l~i>i, 
B~ Bg l~ui, B~ ^ l~i'i and B~ — > l~i'i transitions (l = e, /i). In Ta- 
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Table 4.13: Branching ratios in % for exclusive nonleptonic decays of the 
meson. Our central values have been obtained with the ALl potential. 



B-^B°l-Di 


^miii 




ymax 




o2. 




ymax 


F+ 
This work 

[166] 


-0.39t° 
-0.58 
-0.39 


03 

03 


-0 70+0-°°'' 

^- '"-0.02 

-0.96 
-0.96 


F+ 
This work 
\16A] 
[166] 


-0.58^° 
-0.61 
-0.50 


01 

02 


-0.92 
-0.99 


F_ 
This work 

\m 


-O.lltg 
2.14 


02 
03 


-0 09+° °^ 

^•"^-0.06 

2.98 


F^ 
This work 
[164] 


-0.08t|^ 
1.83 


01 
02 


-0 n5+°-°i 

"■^3-0. 03 

2.35 


Fo 

This work 
[166] 


-0.39t|^ 
-0.39 


03 
03 


-0 71 +"-"5 

'J- ' ^-0.02 

-0.80 


Fo 

This work 
[166] 


-0.58_o 
-0.50 


02 


-0.86^rol 
-0.86 



Table 4.14: F+, F_ and Fn ev aluated at q'^-^^ a nd q ^^^ compared to the ones 



obtained by Ivanov et al. 164i | and Ebert et al. 166l |. Our central values have 
been obtained with the ALl potential. Here I stands for I = e, fi. 



bles 14.181 14.191 we give the total decay widths and branching ratios and compare 
them with determinations by other gro ups. Our results are in bet ter agreement 
wit h th e ones obtaine d by Ebert et al. 166l |. Colangelo et al. 176], Anisimov et 



al. [177| and Lu et al. [181 1. 



4.5.3 Heavy quark spin symmetry 

In Fig. 14.121 we give Tif' ^ and ^ of the B~ B^ and B~ — > i?^ transitions, 

and \ 11^2 \ ^ ^iid t!~^ ^ of the B~ — > B*'^ and B~ B*^ transitions. 
We can take the infinite heavy quark mass limit on our analytic expressions 
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ALl 

0,2 0,4 0,6 0,8 

[GeV^] 

Figure 4.10: -F+ and form factors (solid lines) for B~ — > B^ and B~ —>■ 
semileptonic decay, and V (solid line), (dashed line), A„ (dotted line) and 
^0 (dashed-dotted line) form factors for B~ —>■ B and B~ — > B^ semileptonic 
decays evaluated with the ALl potential. For the first two cases, and for compari- 
son, we also show with dotted lines the results obtained with the Bhaduri (BHAD) 
potential. 




with the result that near zero recoil 



SI 



-(1-) 



.(0-) 
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(4.63) 



When compared to the results of HQSS by Jenkins et al. [T^] we see differences. 
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Table 4.15: V, j4„, a nd vl n evaluated at q ^j^ a nd g^^x compared to the 
ones obtained by Ivanov et al. 1641 ] and Ebert et al. |166t |. Our central values have 
been obtained with the ALl potential. Here / stands for / = e, /i. Asterisk as in 
Table [ 
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Table 4.16: Partial helicity widths in units of 10 GeV for decay. Central 
values have been evaluated with the ALl potential. 
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Table 4.17: Forward-backward asymmetry. Our central values have been evaluated 
with the ALl potential. We would obtain the same results for 5+ B^ decays. 
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Table 4.18: Decay widths in units of 10 GeV. Our central values have been 
evaluated with the ALl potential. 
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Table 4.19: Branching ratios in % . Our central values have been evaluated with 
the ALl potential. 



In Ref. [71 we finc^!! S 



(1-) 
2 



-(0-) ; 



instead. This is wrong as there is a misprint 



Note the different notation and global phases used. 
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Figure 4.11: Differential decay width for the for the B^ I ui and B^ — s- 

P/, B~ — > B* l~9i and B~ — > I'vi (l = e, fi) transitions. Solid line: 



results for a final e evaluated with the ALl potential; dotted line: results for a 
final e evaluated with the Bhaduri (BHAD) potential; dashed line: results for a 
final fj, evaluated with the ALl potential; dashed-dotted line: results for a final ^ 
evaluated with the Bhaduri (BHAD) potential. 



in Ref. [72] that has not been noted before: the sign of the term inv^ in the 



last expression of Eqs. (2.9) and (2.10) in Ref. [T^] should be a minus [l90|. Also 



from Ref. [t^] one would expect = T,^ ^ contradicting our result in 

Eq. ()4.63p were we find ^ = \ Our result is a clear prediction of the 

quark model and comes from the extra signs that appear due to the fact that it is 
the antiquark that decays (See appendix IE.2p . This difference between quark and 



^*One would have to look at Eq. (2.10) in Ref. [72 ]. even though it refers to decay into 
because that is the reaction where you have antiquark decay in their case. 
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Figure 4.12: 
-ftO 



(0-) 



(solid line) and 



(0-) 



B^ B'^g transitions, and Y,\^ ' (solid line), S2 

line) and 113^ ^ (dashed dotted line) of the B~ — 
evaluated with the ALl potential. 



> B^ and 

^^^'^ (dashed line), (dotted 



(dashed line) of the B^. 



B*^ and B: 



B*^ transitions 



antiquark decay was not properly reflected in their published work 19Cl| |. 



How far are we from the infinite heavy quark mass limit? In Fig. I4.13l we show 



of the semileptonic B^ 



^ and B^ 
^0 



B^g transitions, and ^ of the 

— *u — *o 

semileptonic B~ — > B and B~ — > B^ transitions. The differences between the 
corresponding T,^^ ^ and T,^ ^ are at the level of 10%. The differences are much 
more significant for ^2^^ \ ^ and ^ that we show in Fig. 14.141 In each 
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Figure 4.13: Tif' ^ of the semileptonic B' (solid line) and B^ 

(dotted line) transitions, and Ti\ of the semileptonic B~ B (dashed line) 
and B~ B^ (dashed-dotted line) transitions evaluated with the ALl potential. 



case the three curves shown would be the same in the infinite heavy quark mass 
limit. Clearly in this case corrections on the inverse of the heavy quark masses 
seem to be important. 



4.6 Nonleptonic ^ B Mp two-meson decays 



In this section we will evaluate decay widths for nonleptonic B^ B Mp two- 
meson decays where Mp is a pseudoscalar or vector meson with no b quark content, 
and, at this point, B represents a meson with a b quark. These decay modes involve 
ac— >dorc— >s transition at the quark level an d they ar e governed, neglecting 



penguin operators, by the effective Hamiltonian 162l . Il66l | 



H, 



eff. 



^ {Vcd [ci(/i) Qf + C2(m) Qi'] + Vcs [Cl(/i) Qf + C2(/X) Qf] + H.C.] 



(4.64) 



4.6. NONLEPTONIC B- 



BMf TWO-MESON decays 
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Figure 4.14: —T,^ ^ (solid line) of the semileptonic B^ — B^ and B^ — > tran- 
sitions, and ■* (dashed line), ■* (dotted line) of the semileptonic B~ B*^ 
and B~ — > B^ transitions evaluated with the ALl potential. 



where ci, C2 are scale-dependent Wilson coefficients, and Q'f, Q'i^ Qf*, are 
local four-quark operators given by 



Qf 
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(4.65) 
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We shall work again in the factorization approximation taking into account 
the Fierz reordered contribution so that the relevant coefficients are not ci and C2 
but the combinations 



1 1 

al{^J') = Clin) + — C2 ifJ,) ; a2 (m) = C2 (/i) + 7^ (//) (4.66) 



The energy scale app ropriate in this case is /i ~ rric and the values for ai and 
a2 that we use are 162l | 



ai = 1.20 ; a2 = -0.317 



(4.67) 



Mf = IT-, p", K-, K* 



In this case B denotes one of the B^ , B*^ , B^, B*f . The decay widths are 



Mp = case 

^ = l6^"^^^" l"^" '^^ "^^^^"^^ ("^^^ 

Mp = 1^ case 

r = -^— ^V^rf \Vf\ nipfp 

167rm^ 2mB^ 



X ( Wjur (m^) + (m^) + Wof (ml 



(4.68) 



and similarly for iJ^, i?*" with \Vcd\ —>■ {Vcsl', B*^ —>■ S^, B*f . Vp = Vud or 

yp = depending on whether Mp = 7r~ , p~ or M^? = K*^ , fp is the 

decay constant of the Mp meson, and the different Tirr have been evaluated at 

= m'jp. In Table [^^20] we show the decay widths for a general value of the Wilson 

coefficient ai, and the corresponding branching ratios evaluated with oi = 1.20. 

— *o 

The transition B^ K*~ is not allowed with the new Br mass value from 



Ref. [182[ |. Our branching ratios for a final B^ or B^ are in very good agreement 

^0 -Fr*0 



with the results by Ivanov et al. jl62l |. while for a fin al B or B we are in very 
good agreement with the results by Ebert et al. 1661 ] (with the exception of the 



Be —>■ B*^7r decay). 



4.6. NONLEPTONIC B- 



B Mp TWO-MESON DECAYS 
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Here the generic name B stands for a B or a B*^ meson. The different decay 
widths are given by 

Mp = 0~ case 

r _^IT/ |2|T/|2 2 ^'^'("^^c'"^|-,B»-^"^f) 2 72 ^B^->B-,B^-, 2^ 
^=1^"^"^" l^^l ("^^^ 

Mi? = 1~ case 

(4.69) 

where Vp = Vcd or = Vcs depending on whether Mp = 11^, or Mp = K^, K*^ , 
fp = fp for Mp = K^, K*^ whereas Jp = fp/y/2 for Mp = 1:^, p^, with fp the 
Mi? meson decay constant, and the different T-irr evaluated at = nip. The 
latter have been obtained from the matrix elements for the effective current oper- 
ator ^'2(0)7^(1 — 75)^^(0). The decay widths, for a general value of the Wilson 
coefficient 02 , and the corresponding branching ratios are shown in Table 14.211 
With the exception of the B^ B*^i: ^ ca se, our results are in a global good 
agreement with the ones by Ebert et al. [l6i |. 
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0.29to.ol 0.29 0.21 0.41 0.17 1.06 0.21 
0.013-0.001 0.011 0.0030 0.10 0.0043 


B- ^B*%- 
B- ^Bf p- 
B~ B*f K- 


2.34t°;;t;^ 2.1 1.58 5.08 1.23 6.5 3.50 1.95 
13.4lgj 11 10.8 14.8 16.8 20.2 10.8 12.1 
0.13lo;ol 0.13 0.11 0.29 0.13 0.37 0.16 



Table 4.20: Decay widths in units of 10~ GeV, and for general values of the 
Wilson coefficient ai, and branching ratios in % for exclusive nonleptonic decays 
of the B~ meson. Our central values have been obtained with the ALl potential. 
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r [10-15 Gcv; 





This work 


B- ^ B- tt" 
B- B- p° 
B- B-^ 
B- B- K*° 


^■0^-0. 07 "2 
^■'^-0^10 "2 

35.3^ ^ g 

13.1 7 ^2 


Be B*- ttO 
Be ^ B*- pO 
Be B*- 
Be B*- K*" 


U-'J"J-0.05 "2 
9 O/I+0.27 2 

16.9l^-^a2 

103t«3a2 



B.R. in % 





This work [162] [1661 UM UM UM [llll 


B- -> B- 
B- ^ B- K" 
B- B- K*° 


0.0038^[;;;'oo6 0-007 0.004 0.011 0.004 0.037 0.007 
O.OOSOlg-jJgg^ 0.0071 0.005 0.020 0.010 0.034 0.009 
0.25i[!-^^ 0.38 0.24 0.66 0.27 1.98 0.17 
0.093t!^;Sj;^ 0.11 0.09 0.47 0.32 0.43 0.095 


Be B*- vr" 
Be B*- p° 
Be B*- 
Be B*- 


o.oo25i^:;j;;;]| 0.0020 0.001 0.010 0.003 0.033 0.004 

0.020ti]S]o3 0.011 0.024 0.041 0.031 0.09 0.031 
0.12l°°2 0.088 0.11 0.50 0.16 1.60 0.061 
0.73l^'?[; 0.32 0.84 0.97 1.70 1.67 0.57 



Table 4.21: Decay widths in units of 10~ GeV, and for general values of the 
Wilson coefficient 02, and branching ratios in % for exclusive nonleptonic decays 
of the B~ meson. Our central values have been obtained with the ALl potential. 



Chapter 5 



Doubly heavy baryons 
spectroscopy and static 
properties 



5.1 Introduction 



The subject of doubly heavy baryons has been attracting attention for a long 
time. Magnetic mom ents of doubly charmed baryons were evaluated back in the 
70's by Lichtenberg [19l[ | within a nonrelativistic approach. The infinite heavy 
quark mass limit was already used in the 90's to relate the s pect rum of doubly 
heavy baryons to the one of m esons with a single heavy quark |l92l | , or to analyze 
their ser nilep tonic decay 193| |. A factor of two error in the hyperfine splittings 
of Ref. [192] have been recent ly n oticed by the potential nonrelativistic QCD 



(pNRQCD) calculation of Ref. [194^ 



On the experimental side the SELEX Collaboration claimed evidence for the 
baryon, in the A+iir~7r+ and pD^K~ decay modes, with a mass of M„+ = 
3519 ± 1 MeV/c^ |l95l ]. Those results were challenged by a theoretical analy- 



sis [l96l | which claimed the observed events by the SELEX Collaboration could be 
explained without the involvement of doubly charmed baryons. Other experimen- 
tal coUaborations like FOCUS [igj], BABAR and BELLE have found 
no evidence for doubly charmed baryons so far. At present theH^ has only a one 
star status and it is not listed in the particle summary table 188l ]. 

As discussed before, for hadrons with two heavy quarks one can apply HQSS [7: 
which amounts to the decoupling of the heavy quark spins in the infinite heavy 
quark mass limit. In that limit one can consider the total spin of the two heavy 
quark subsystem (Sh) to be well defined. In this work we shall assume this is 
a good approximation for the actual heavy quark masses. This approximation, 
which is the only one related to the infinite heavy quark mass limit that we shall 
use, will certainly simplify the solution of the baryon three-quark problem. Re- 
cently the authors of Ref. [194l | have developed and effective theory (pNRQCD) 
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Table 5.1: Quantum numbers of doubly heavy baryons analyzed in this study. S, 
are strangeness and the spin parity of the baryon, / is the isospin, and S*^ is 
the spin parity of the heavy degrees of freedom. / denotes a light u or d quark . 



suitable to describe baryons with two and three heavy quarks. 

Solving the three-body problem is not an easy task and here we shall do it by 
means of a variational approach. The approach, with obvious changes, was a lready 
applied with good results in the study of baryons with one heavy quark 200( 1 . This 
method, that leads to simple and manageable wave functions, is made possible by 
the simplifications introduced in the problem by the fact that we can consider Sh 
to be well defined. 

Our simple variational calculation reproduces the results for static proper- 
ties obtained in Ref. [tH] by solving more involved Faddeev type equations. Our 
method has the advantage that we provide explicit and manageable wave func- 
tions that can be used to evaluate further observables. Static properties like 
masses and magnetic moments of doubly heavy baryons have also been studied 
in other models. Masses have been calculated i n th e relativistic quark model as- 
suming a light quar k he avy diquark structure 20ll |. the potential approach and 
sum rules of QCD j202l ]. the nonperturbative QCD approach 203], the Bethe- 
Salpeter equation applied to the light quark heavy d iqua rk [ 2o3 |, the nonrelativis- 
tic quark model wi th harmo nic oscillator potential 205l | or with the use of QCD 
derived potentials [2] 



2071], the relativistic quasi-potential quark model [2081 ] . 



with the use of the Feynman-Hellman theorem and semi-empirical mass formulas 



5.2. THREE BODY PROBLEM 
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Figure 5.1: Definition of different coordinates used through this work. 



within the framework of a nonrelativistic constituent quark model [209, l210(], in 
effective field theories (g^, 194], or in latti c e nonrela tivistic QCD [211i] . There 
are also lattice QCD determinations 212 . 2131 . 214]. S imil arly, magnetic mo- 
ments have been evaluated in a nonrelativistic approach 1911 ]. in the relativistic 
three-quark model 215|, th e relativistic quark model usin g different forms of the 
relativistic kine matics f216^, in the skyrmion model 217 ]. in the Dirac equation 
formalism 2181 1. or using the MIT bag model |219l ]. 

In Table [5T] we summarize the quantum numbers of the doubly heavy baryons 
considered in this study. 



5.2 Three Body Problem 
5.2.1 Intrinsic Hamiltonian 

In the Laboratory (LAB) frame (see Fig. 15. 1|) . the Hamiltonian [H) of the three 
quark (/ii, /12, <?, where /ii, /12 = c, 6 and q = l{u, d), s) system reads: 



H= E \^^-^A+y>^^>^2 + Vh,, + Vn,, (5.1) 
j=/ii,ft,2,g 



where ruh^ , m/^j , ruq are the quark masses and the quark-quark interaction terms 
Vjk depend on the quark spin-flavor quantum numbers and the quark coordinates 
{xhi,Xh2, Xq for the hi, /12, q quarks respectively). To separate the Center of Mass 
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(CM) free motion, we go to the hght quark frame {R, 7^1,7^2), 

^ _ ruh^Xh^ + ruh^Xhi + mgXg 

'"1 — 2^/11 Xq 

r2 = Xh2 - Xq (5.2) 

where R and n, r2 are the CM position in the LAB frame and the relative positions 

of the hi, /i2 heavy quarks with respect to the light quark q. The Hamiltonian 
now reads 

H = + (5.3) 
2M 



i=i,2 

^sp ^ _ NO ^ Vh^girj, spin), 3 = 1,2 (5.4) 



2^j 

where M = ruh^ + mn^ + mq, ^ij = {l/irih^ + l/ruq) ^ and Vj= i9/%, J = 1, 2. 
The intrinsic Hamiltonian i/™* describes the dynamics of the baryon. Apart from 
the sum of the quark masses M, it consists of the sum of two single particle 
Hamiltonian H^^ , each of them describing the dynamics of a heavy-light quark 
system, plus the heavy-heavy interaction term, including the Hughes-Eckart term 

(Vi • V2)- We will use a variational approach to solve it. 

5.3 Variational Wave Functions 

For the interactions presented in chapter [2l we have that both the total spin and 
the internal orbital angular momentum given as 

5" = {oh^ + + <?g) /2 

L = li + I2, with Ij = — i fj X Vj, j = 1,2 (5.5) 

commute with the intrinsic Hamiltonian and are thus well defined. We are inter- 
ested in the ground state of baryons with total angular momentum J = 1/2, 3/2 
so that we can assume the orbital angular momentum of the baryons to be L = 0. 
This implies that the spatial wave function can only depend on the relative dis- 
tances ri, r2 and ri2 = |n— ^2!- Furthermore when the heavy quark mass is infinity 
(m/i 00), the total spin of the heavy degrees of freedom, 5heavy = {^hi + '?h2) /2) 
commutes with the intrinsic Hamiltonian, since the spin-spin terms in the poten- 
tials vanish in this limit. We can then assume the spin of the heavy degrees of 
freedom to be well defined. 
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With these simplifications we have used the following intrinsic wave functions 
in our variational approactfl 

• Ehj^h2, ^hih2-type haryons: 



2' 2 



Ms^Ms.Mj \ \hih2;lMs^) ® \q;\Ms,)^l[^^{ri,r2,ri2) 



1 



(5.6) 



where Mj is the third component of the baryon total angular momentum 
while \ h1h2; Sh, Msf^) and \q; \Msq) represent spin states of the /11/12 subsys- 
tem and the light quark respectively. For hi = /12 we need ^^'^^ {ri , r2,ri2) = 

^/iihi (^2) '^i) ^'12) to guarantee a complete symmetry of the wave function un- 
der the exchange of the two heavy quarks. 

• '^h^h2^ Kih2-tyP'^ baryons: 
3 



^ ( 2' 2 

Ms, Ms, ^ 



Ms, ,Ms^,Mj\ I /ii /i2 ; IM5 J C3 I (?; ^Mcj, ) ^^kT (^1 ' ^2 > ^2 ) 



(5.7) 



Similarly to the case before for hi = h2 we need ^'^^^^ (?'i,''2,n2) = 

*hihr('^2,n,n2). 

• -h,h2^ ^'h,h2-^yP^ baryons: 

\^'h,h2^^'h,h2-^J=\^Mj) = \hih2;m)®W,\Mj)^i;^^{ri,r2,ri2) (5.8) 

In this case hi 7^ /i2 and we do not need the orbital part to have a definite 
symmetry under the exchange of the two quarks. 

The spatial wave functions ^'(ri, r2, ri2) in the above expressions will be de- 
termined by the variational principle: 6{B\H'^'^^\B) = 0. For simplicity, we shall 
assume a Jastrow-type functional fornH: 



*fih2(n,r2,ri2) = NF^{ri2)(t^hAn)(t>h2,{r2) (5.9) 

^We omit the antisymmetric color wave function and the plane wave for the center of mass 

motion which are common to all cases. 

^ A similar form lead to very good results in the case of baryons with a single heavy quark [2O0I ] . 
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where is a constant, which is determined from normahzation 

r+oo r+oo r+1 

= Svr^ / dnrl / drarl/ dfi \^f,^^,^{n,r2,ri2)\ (5.10) 
Jo Jo J-i 

with /X being the cosine of the angle between the vectors ri and r2 (ri2 = ( ''i + 
rl — 2rir2/i)^/^). 

The functions (phiq and 4>h2g will be taken as the 5"— wave ground states ^j{rj) 
of the single particle Hamiltonians fZj^ of Eq. (j5.4p modified at large distances 

(phjqirj) = {l + ajrj)ipj{rj), j = 1,2 (5.11) 

The heavy-heavy correlation function will be given by a linear combination 
of gaussian^ 

i^^(ri2) = «i = l (5-12) 

The value of one of the aj parameters can be absorbed into the normalization 
constant A^, so that we fix ai = 1. The values that we get for all parameters are 
given in appendix [HI 



5.3.1 Infinite heavy quark mass limit of variational wave functions 

In the infinite heavy quark mass limit the wave function of the system should look 
like the one for a "meson" composed of a light quark and a heavy diquark. The 
two heavy quarks bind into a 3 color source diquark that to the light degrees of 
freedom appears to be pointlike |l93l ]. In our model the pointlike nature of the 
heavy diquark comes about through the one-gluon exchange Coulomb potential 
which binds the two heavy quarks into a distanc^ 

rhih2 oc ; fJ-hih2 = : (5-13) 

that tends to zero if both quark masses go to infinity. 

For our wave functions we can define the probability distribution Phih2 for the 
two heavy quarks to be found at a distance r/^^^j 



Ph^h2 Kha) = J d-'n J d-'r2 d{n2-rh^h2) \^Zh2i''2,r2,n2)\ (5.14) 

In Fig. 15. 2^ we show the Phih2 probability distributions for the Ecc, H^c, '^bb and 
Qcc, ^6c) ^bb baryons evaluated for the ALl potential. We see how the maximum 



•^Note that should vanish at large distances because of the confinement potential. The 
confinement potential is also responsible for the non-vanishing values of the parameters Uj , j = 
1,2 in Eq. (EHJ. 

''The relation can only be approximate due to confinement and the interaction with the light 
quark. 
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Figure 5.2: Phih2 probabilities (see Eq. (l5.14p ) for the Hcc, '^bc, ^bb and 
f^cc) ^bcy ^bb baryons evaluated with the ALl potential. 



moves, as expected, to lower rh^h^ values as the quark masses increase. 

On the other hand the relative wave function for the light quark with respect 
to the heavy quark subsystem should tend to the one of a light quark relative to 
a pointlike diquark. This limit is not evident in the coordinates we work as one 
has first to separate the diquark internal orbital wave function from the total one. 
To see that this limiting situation is in fact reached within our variational ansatz 
let us introduce the new set of coordinates 

^ _ rnh^Xhi + ruh^Xhi + rngXg 
rrihi + nih^ + nig 
ri2 = Xh-, - Xh2 
^ _ ruh^Xhi +rnh2Xh2 ^ _ rrih^n + mh2r2 



in terms of which the Hamiltonian now reads 

H = -X^ + ^-t 
2M 

rint 



(5.15) 



W'^' = M + Hh,h2+Hgh,h2 (5.16) 



where 



-^2 
Vl2 



Hhih2 = -7, VVh^h2{n2, spin) 

If 1 l\-^2 TO,, ^ 

Hqh,h2 = -o L^ ,^ ^7;r)^i + ^ I ^ ^12, spin) 

+yh2q{rq ri2, spin) (5.17) 

mh^+ mh2 
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^cc 








^bb 


\V\^ 0.974 


0.975 0.991 


0.959 


0.966 


0.984 



Table 5.2: Absolute value square of the V projection coefficient defined in 
Eq. ([CTD 

with Vi2= d/df^^ and Vg= d/dft^. Defining now 

Hqh,h2 = -\ ( ^ ]^ + Vg +Vh^q{rq, Spin) + Vh^q{fq, Spiu) (5.18) 
^ \ '"'/ii ~r '"'/i2 '"'9/ 

one would have 

= i??ii/i2 + -ff°hi/i2 + {Hqh^h2 - Hqh^h^) (5.19) 

Hfij^fi^ is the Hamiltonian for the relative movement of the two heavy quarks while 
^qhih-z Hamiltonian for the relative movement of the light quark with respect 

to a pointlike heavy diquark where the two heavy quarks are located in their 
center of mass. Both movements are coupled through the term {Hqh^h2 ~ ^qhih2)- 
If the heavy quark masses get arbitrarily large the average ri2 value tends to 
zero and thus one can neglect the effect of Hqh^h^ — -f^g/j^/j^. In that limiting 
situation the light and heavy quark degrees of freedom decouple completely and 
the internal Hamiltonian reduces, as it should, to the sum of the Hamiltonians 
Hhih2? corresponding to the relative movement of the two heavy quarks, and 
-^g/ii/i2' corresponding to the relative movement of the light quark with respect to 
the pointlike heavy diquark subsystem. 

As to the wave function it should reduce in that limit to the product ^hih2{''"i2)- 
i/i2 ('"'?)' being <I'hih2(n2), *^ghi/i2(^g) ^^e ground state wave functions for Hh^h^, 
^qhih2 respectively. To check that we reach that limiting situation we have evalu- 
ated the projection V of our variational wave functions onto ^hih2 (^"12) • ^qfuJi2 ('''?) 
evaluated with the actual heavy quark masses. This projection is given bjo 

V = I dVi I d'T2 {^U,{n,r2,n2y ^h,h2{n2) ^lh,h2i^,) (5.20) 

and the values for \P\^ that we obtain for the Hcc, '^bc-, '^bb and Qcc, ^bc^ ^bb 
baryons using the ALl potential are given in Table 15.21 We see how ["Pp in- 
creases with increasing heavy quark masses indicating that for very high heavy 
quark masses the total wave function tends to the one of a light quark relative 
to a pointlike diquark times the diquark internal wave function. In summary, our 
variational wave functions respect the infinite heavy quark mass limit. 

We note by passing that the product ^hih2{''"i2) • ^qhih2(^'i^ corrected by a 
correlation function in the variable ri2 • rq would have been a good variational 



*Note Tq can be expressed in terms of ri,r2,ri2 and that d^rid^r2 = d^ri2d'^rq 
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orbital wave function. For instance, a calculation using the ALl potential and an 
orbital wave function given just by that product ^hih2{^i2) • ^^g/ii/ijC'^'?) gives for 
the expectation value of ff™* for the Hcc, S^c, '^tb baryons 

/^intu^ii ^ 3640 MeV; = 6943 MeV; = 10198 MeV (5.21) 

which are respectively 28 MeV, 24 MeV, and 1 MeV larger, and therefore worse, 
than our best values in Table 15.31 The correlation function would clearly be 
needed, being its role more important for the "cc" and "6c" systems. The im- 
provement in going from a "cc" to a "6c system" is not as good as one would 
naively expect, the reason being that, considering fi2 to be a small quantity, 
Hqhih2 ~ -f^5hi/i2 ^^^^ order in fi2 for the "6c" system while, for symmetry rea- 
sons, it is necessarily of second order for the "cc" and "66" ones. In any case one 
sees how the expectation values obtained with this simple ansatz get closer to our 
best values in Table 15.31 as the heavy quark masses increase. 

5.4 Static Properties 
5.4.1 Masses 

The mass of the baryon is simply given by the expectation value of the in- 
trinsic Hamiltonian. Our results (VAR) appear in Table 15.31 where we com- 
pare them with the ones obtained in Ref. [751] with the use of the same inter- 
quark interactions but within a Faddeev approach (FAD). For that purpose we 
have eliminated from the latter a small three-body force contribution of the type 



^123 = constant /m/i^m/ij'Tig that was also included in the evaluation of Ref. [751]. 
We will show their full results in the following tables. Whenever comparison is 
possible we find an excellent agreement between the two calculations. In some 
cases the variational masses are even lower than the Faddeev ones. Besides we 
give predictions for states that were not considered in the study of Ref. [t^ . 

In Tables 15.41 and 15.51 we compare our results with other theoretical calcula- 
tion^. Our central values correspond to the results obtained with the ALl poten- 
tial, while the errors quoted take into account the variation when using different 



Note in Refs. [2091 . |210| | the Hi,c, H^^. and Qbc, i^'bc baryons are defined such that the total 
spin of the hght q quark and the heavy c quark are well defined, being for S.bc, f^bc and 1 for 
Hjij., f2j^. They are thus linear combinations of our states. The different spin functions are related 
by 

1 \ 1/2 , / 1 \ '''^ a/3 / 1 ^ 

!gc;0)®|6;i)j = i (^|6c; 0) ® Ig; i) j - ^ 1) ® jg; i) 

1 \'/^ a/3 / 1 X'/^ 1 / 1 ^ 



gc;l)®16;-)J = - 1- i^\bc; 0) \q; -) j _ _ |^|6c; 1) ® |g; -) j (5.22) 

In order to extract their predictions for the Hbc, Hj^. and ^Ibc, ^ic baryons with total spin of the 
two heavy quarks well defined, we have assumed that the above relations, but with coefficients 
square, are also valid for the masses. Note this may be incorrect as we are neglecting a possible 
interference contribution. 
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ATI 


A T O 

ALZ 


A t;> 1 
Ar i 


A TDO 

Arz 


T)TT A T^ 


^cc 


VAR 
FAD [75] 


3612 
3609 


3619 
3616 


3629 
3625 


3630 
3628 


3639 
3633 




VAR 


3706 


3715 


3722 


3729 


3722 


'^bb 


VAR 

FAD [75] 


10197 
10194 


10180 
10175 


10207 
10204 


10179 
10176 


10202 
10197 


^bb 


VAR 


10236 


10219 


10245 


10219 


10235 




VAR 
FAD [75] 


6919 
6916 


6912 
6913 


6933 
6928 


6917 
6907 


6936 
6934 


^bc 


VAR 


6948 


6942 


6957 


6944 


6965 


^bc 


VAR 


6986 


6981 


7000 


6987 


6993 




ATI 


A T O 


A TD 1 

Ar i 


A T>0 


T3TT A T^ 




VAR 
FAD [75] 


3702 
3711 


3718 
3718 


3711 
3710 


3710 
3709 


3743 
3741 




VAR 


3783 


3802 


3800 


3802 


3805 




VAR 
FAD [75] 


10260 
10267 


10249 
10246 


10259 
10258 


10226 
10224 


10274 
10271 


0* 
^^bb 


VAR 


10297 


10287 


10301 


10269 


10302 


^bc 


VAR 
FAD [75] 


6986 
7003 


6986 
6996 


6990 
6996 


6969 
6971 


7013 
7023 


^'bc 


VAR 


7009 


7010 


7011 


6994 


7033 




VAR 


7046 


7047 


7055 


7037 


7057 



Table 5.3: Doubly heavy H and baryons masses in MeV. VAR stands for the 
results of our variational calculation. FAD stands for the results obtained in 
Ref. ^] using the same interquark interactions but within a Faddeev approach. 
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potentials. The same presentation is used for the results obtained in Ref. [75l | 
for which we now show their full values including the contribution of the three- 
body force. All calculations give similar results that vary withi n a few per cent. 
From the experimental point of view the SELEX Collaboration 195f | has recently 
measured the value of This experimental value is 100 MeV smaller that 

our result. On account of what has been said in the introduction to this chapter 
, one sho uld take this experimental value with due caution. Note also that in 
Ref. |l95l | the systematic error is not given. There ar e also different lattice deter- 
minations for baryons with two equal heavy quarks 2121 . 213l . 214]. Our results 
agree within errors with the lattice data for baryons with two c quarks, while they 
are roughly 100 MeV below lattice results for doubly 6-quark baryons. The best 
over a,ll agreern ent with lattice data available so far is achieved in the calculation of 
Ref. [209l . l210l | where they use the Feynman-Hellmann theorem and semiempirical 
mass formulas in the framework of a nonrelativistic quark model but without the 
use of an explicit Hamiltonian. With i n fu ll dyn amical calculations, ours and the 
relativistic calculations in Refs. 201 . 204, 20a] have the best overall agreement 
with lattice data. 
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6820 


6950 


6965 ± 901" 




6840 ± 236 


—be 


69481b^ 




6963 


6850 




7070 






6850 


7000 


7065 ± 90+ 






—be 


69861" 




6980 


6900 




7100 






6900 


7020 


7060 ± 90 







This work Exp. [195] Latt. [212] Latt. [213] Latt. [214] 
3612+1'' 3519 -j- I 3605 ± 23 3549 ± 95 

3706+^^ 3685 ± 23 3641 ± 97 

10197li? 10314 ± 47 

10236l?7 10333 ± 55 



Table 5.4: First panel: doubly heavy H masses in MeV as obtained in different 
models. Our central values, and the ones of Ref. 75|, have been evaluated with 
the ALl potential. Second panel: we compare our resul ts w ith the experimental 
value for M^^^ measured by the SELEX Collaboration [195] (Note the cautions 
that appear on this expe rime r ital mass in the introduction to this chapter), and 
lattice results from Refs. (ill, I213I . \214\ . Entries with ^ should be taken with due 
caution (see text). 



Ther e are also independent determi natio ns of mass splitti ngs i n lattice QCD 2121 . 

2ll,|2li], nonrelativisitic lattice QCD and pNRQCD [l9i]. In Table [52] we 
compare those results to the ones obtained in the present calculation and in other 
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This work [75] [201] [202] [203] [204] [205] [207] f208] [209,210] [65] [211] 

TlZ 3702+"" 3710+f 3778 3590 3860 3760 3749 3590 3760 3740 ± 70 3710 3698 ± 65 

3783+^^ 3872 3690 3900 3826 3690 3890 3820 ± 80 3760 

flib 10260t34 102 67143 10359 10180 10340 10340 10180 10320 10370 ± 100 

nib 102971^8 10389 10200 10380 10200 10360 10400 ± 100 

ribc 69861 17 7OO3I32 7088 6910 7130 7050 6910 7050 7045 ± 90^ 6954 ± 225 

ni^ 70091 15 7116 6930 7110 6930 7090 7105 ± 90^ 

n^c 7046I9" 7130 6990 7130 6990 7110 7120 ±90 



This woA Latt. [212] Latt. [213] Latt. [214] 

TlZ 3702+*" 3733 ±9l38 3663 ± 97 

3783+^^ 3801 ±9l34 3734 ± 98 

Qtb IO26OI34 10365 ±40+12 lo** 

IO297I28 10383 ± 39+|lJ^ 



Table 5.5: Same as Table [53] for doubly heavy baryons. 



models. 



Our central results e valuated with the ALl potential are larger tha n th e ones 
obtained in lattice QCD [ilj . [2li . l2li ] and lattice nonrelativistic QCD [ml. The 
agreement is better when we use the BHAD potential of Ref. 73|] for which we get 
the l owes t results. Similar results are obtained by the rel ativistic calculation of 
Ref. [201 1 . whereas for the rela tivistic calculations in Refs. [ioi, [ioi, I2OS] and the 
nonrelativistic one in Ref. 207l | the agreement with lattice QC D and non relativistic 
lattice QCD data worsens. As for the calculations in Ref. j209l . bid ] we do not 
quote their results due to the large theoretical errors involved. 







This woA [201] 


[202J [204] 


[207] [208] 


[19J 


[21^ 


Latt. [212] 


Latt. [213] Latt. [214 


Ms- 

cc 




94t?i 


107 


130 


120 


132 


150 


120 ± 40 


70± 13 




80 ±11 87 ±19 


"bb 
~bc 


-Ms,, 


391^ 


35 


40 


40 


40 


50 


34 ±4 


20 ±7 


20 ±6 




-Ms,'. 


67l?o 


47 


80 


90 


80 


70 




43± 11 






~bc 


-Ms,^ 


29+^ 


30 


30 


60 


30 


50 




9± 7 










8lli^ 


94 


100 


140 


100 


130 




63 ±9 




68 ± 7 67 ±16 




-Mn,, 


371^ 


30 


20 


40 


20 


40 




19 ±5 


20 ±5 




-Mn,^ 


60l?6 


42 


80 


80 


80 


60 




39 ±8 








- Afo.c 


231^ 


28 


20 


60 


20 


40 




9±6 







Table 5.6: Ivlass splittings in MeV for doubly heavy H and Vl baryons. Qur central 
values have been obtained with the ALl potential. 
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5.4.2 Charge densities and radii 

The baryon charge density at the point P (coordinate vector r in the CM frame, 
see Fig. 15. ip is given by: 

{ehiS^{r- yh^) + et^S^ir- yu^) + eq5^{r- yq)] 

= Pe{n\h,+Pe{r)\h,+Pe{r)\q (5-23) 

where eh^,h2,q quark charges in proton charge units e, and from Fig. 15.11 

we havfllfyX^ = Vq + n, yn^ = l/g + ^2 and yq = - {miin + mh2r2) /M. Since our 
L = wave functions only depend on scalars {ri,r2 and ri2) the charge density is 
sphericahy symmetric {Pe{r) = (|^ |))- 
The charge form factor is defined as usual 

:F^{q) = J dVei^>f(r) (5.24) 

and it only depends on \q\. Its value at g = gives the baryon charge in units of 
the proton charge. 

The charge mean square radii are defined 

/r+00 
d\ r^p^ir) =Att j dr r^p^{r) (5.25) 

In Figs. 15.31 15.41 and 15.51 we show the charge form factors for the different doubly 
heavy baryons under study including the two different charge states for the doubly 
heavy S ones. We show the calculations with both the ALl potential and the 
BHAD potential. The differences between the two calculations are minor in most 
cases. 

In Table [5^71 we show the charge mean square radii. With the exceptions of the 
Hfjj, and fij^^, we find good agreement with the results obtained in Ref. |75|] within 
a Faddeev calculation. The possible presence of a 5/^ = contribution in the wave 
functions of Ref. [75[] could be the possible explanation for this disc repancy. We 



also compare with the results obtained, for a few states, in Ref. [216l | with the use 



a relativistic quark model in the instant form. The agreement is bad in this case. 
5.4.3 Magnetic moments 

The orbital part of the magnetic moment is defined in terms of the velocities v of 
the quarks, with respect to the position of the CM, and it reads 

P'^ = y dVidVs (^'f^;,2(^l'^2,n2))* |^^(y/.i X mftjW;,J^ 

+ ^^(y^2 ^ rnh^Vh2)z + ^iVq X rUqVy^)^^ ^'^^^.^(ri, ra, ri2)(5.26) 



^There exists the obvious restriction mhiVhi + mh-iVh^ + mqijq = 0. 
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Figure 5.3: Charge form factor of the Ecc, S^f, and H*^, baryons evaluated 
with the ALl (sohd hne) and BHAD (dashed hne) potentials. We show the two 
possible charge states. 



with ruh^Vh^ = -iVi, mn^Vhi = -iV2 and mgVg = i \ Vi + 

Since our orbital wave function has L = 0, the orbital magnetic moment van- 
ishes. The magnetic moment of the baryon is then entirely given by the spin 



Note that the classical kinetic energy has a term on Vh-^ ■ Vh^ and then the operator nihiVhi 
is not proportional to — i Vy^^ , but it is rather given by 

^h,Vh, = ■ i-iVy^^) ~ ^ ■ i-tVy^J = (-*Vi) (5.27) 

Similarly mh2«h2 = V2)- 
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q [fm ^] q [fm ^] 

Figure 5.4: Charge form factor of the Qca ^bb and Q^^, il^j, baryons evaluated 
with the ALl (sohd hne) and BHAD (dashed line) potentials. 



contribution. 

{B; J, Mj = J\ ^{Bh,), + ^{Bh,)z + ^{B,)z |S; J, Mj = J) (5.28) 
Those matrix elements are trivially evaluated with the results 

^hlh2i "hlh2 ^ o I 7^ 1 



3 2mh2 2 2m^ 



-hih2, ^Ihih, 2m,, + 2m,, + 2m, 

^/xi/x2' ^w/i2 ^ (5.29) 



In Table 15.81 we give our numerical results. Our central values, as the ones 
obtained in Ref. [751] within a Faddeev approach, have been evaluated with the use 
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Figure 5.5: Charge form factor of the Ef,c, H^^, H'^^ and Clhc, ^Ic^'bc baryons eval- 
uated with the ALl (solid line) and BHAD (dashed line) potentials. For the S 
baryons we show the two possible charge states. 
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Table 5.7: Charge mean square radii in fm^ for doubly heavy H and Q baryons. 
Our central values, and the ones of Ref. |75|], have been evaluated with the ALl 
potential. 



of the ALl potential. When compared to the values obtained in Ref. \73\ we find 
very good agreement with just a few exceptions (H^^, H^, ^^J- The discrepancy 
for the latter baryons may come from a possible non negligible Sh = contribution 
to their wave functions in the calculation of Ref. [tH]. In our case we have fixed 
5/1 = 1 which we think is a very good approximation since in the limit of infinite 
heavy quark masses the spin of the heavy quark degrees of freedom is well defined . 

We a lso compare our results to the ones obtained in Refs. 191, bisl . bid . 217 . 

using different approaches. The differences between different calculations 
are in some cases large. Being L = a good approximation the magnetic moments 
are essentially determined by the spin contribution of the quarks. With mi, ^ 
iTT'u, TTT'di iTT's, the Contribution from the b quarks is negligible compared to the one 
of the light quark. This is also true to a lesser extent for the c quark. 



102 



CHAPTER 5. DOUBLY HEAVY BARYONS 1 





This work 


[75] 


[191] 


[215] 


[2161 [217] [218] 


[219] 


' — 'CC 

C++ 

"cc 


U. /00_o.030 
-U.2U«_0 Ogg 


U-'o4_o 029 
U.ZUD^o 086 


0.806 
-0.124 


0.72 
0.13 


0.72 0.89 ~ 0.98 0.778 ~ 
-0.10 -0.47 -0.172 ~ 


0.790 
-0.154 


0.86 
0.17 


-.*+ 

' — cc 

■=*++ 

'cc 


n 01 1 +0.052 

U.O±±_Q J^3Q 

^•O'-0.15 




-0.186 
2.60 




-1.17 0.98 

3.16 ~ 3.18 




0.20 
2.54 


"i)6 
"i)6 


-0 742+0-0'** 

'^^-0.091 
^ ^ —0.021 


-n 742+0-044 

'^^-0.092 

0.251+"°:^? 

0.021 




-0.53 
0.18 


—0.726 ~ 
0.226 ~ 


—0.705 
0.236 


0.61 
0.14 


"i)6 
"i)6 


1 07+0.27 
^■^'-0.13 
1 1 1+0.06 
^•J^J^-0.14 












1.37 
-0.95 


-7O 

"&c 
■3+ 
"be 


0.518«n 
-0.475l°;r8 


o.ossini 

-0.1981°:°^^ 




0.42 

-0.12 








-7/0 

"be 
■3'+ 
' — ^be 


-0.993t°:?f, 






-0.76 
1.52 


-0.385 ~ 
1.50 ~ 


-0.366 
1.54 




-7*0 

"be 
-■* + 

"be 


-0.712t°;«^ 

Q 97+0.27 

-0.14 












-0.39 
2.04 




This work 


[75] 




[215] [2161 


[2171 [218] 


[219j 






0.635+ti;i|^ 


0.635+ill!i^, 


0.688 


0.67 0.72 


0.59 ~ 0.64 0.657 ~ 0.663 


0.84 






, nq+0 009 

U.iJ»_0 017 




0.167 




-0.20 ~ 0.03 


0.39 




f^6"b 


Q ,r,, +0.007 






0.04 


0.105 ~ 0.108 


0.084 




"bV 












-1.28 






0.368tHi? 


0.009«;«i 




0.45 








Q/0 

"6c 


U.Oy:Z_Q Q24 






-0.61 


-0.130 0.125 




Q*Q 


U.ZOi_o 021 










-0.22 





Table 5.8: Magnetic moments, in nuclear magnetons {e/2mp, with nip the proton 
mass ), o f doubly heavy H and O baryons. Our central values, and the ones of 
Ref. [73], have been evaluated with the ALl potential. 



Chapter 6 



Doubly heavy baryons 
semileptonic decay 



6.1 Introduction 

In this chapter we should use the manageable wave functions that we got in chapter 
[5] to study semileptonic decays of doubly, J = 1/2, baryons. All the transitions 
studied here involve a 6 — > c transition at the quark level. We shall evaluate form 
factors, decay widths and angular asymmetry parameters. Previous calculations 
of semilepto nic decay widt hs have been done in differe nt re lativistic quark model 
approaches [iii], or with the use of HQET [iT^. 



6.2 Decay width and angular asymmetries 



The differential decay width for a -6(1/2+) - 
2 „2 d^P' 



AT = S\Vcb\'mB'G'F- 



B'{l/2^) transition reads 



i2TT)W^,ip') (2^)32^p,(fc) (27r)32S^'(fe') 

X {27t)''6\p -p'-k- k')C^^{k, k')na^{p,p') (6.1) 

where \ Vcb\ is the modulus of the corresponding CKM matrix element, m^' is the 
mass of the final baryon, p, p' , k and k' are the four-momenta of the initial baryon, 
final baryon, final anti-neutrino and final lepton respectively, and C and TC are the 
lepton and hadron tensors. 

The lepton tensor is given as 

C''''{k,k') = k'f'k'' + k"'k'' - g'"'k-k' -ief"''"^k'^k/3 

As for the hadron tensor we have 



(6.2) 



n^AP,P) = ^^(i?',r'p-"|^'^(0)7^(/-75)^'(0) 



B,r p 



B',r' p' #(0)7,(/- 75)^-^(0) B,rp) (6.3) 
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with \B, r p) {\B' , r' p')) representing the initial (final) baryon with three-momen- 
tum p {p') and spin third component r (r'). The baryon states are here normalized 
such that (r p\ r' p') = {2tt)'^ {E(p)/m) 5„' 6^{p—p'). The hadron matrix elements 
can be parametrized in terms of six form factors as 

B', r' p' I ^(0)7^(1 - 75)^''(0)| r p ; 

= n^'(/){7M (^i(^) - l^Gr{w)) + {F2{w) - -i^G2{w)) 

+v'^ (Fsiw) - j^Gsiw)) }u?{p) (6.4) 

where u^'^' are dimensionless Dirac spinors, normalized as uu = 1, and = 
Pn/iTT-B {v'^ = p'^/msi) is the four velocity of the initial B (final B') baryon. The 
form factors are functions of the velocity transfer w = v ■ v' or equivalently of the 
four momentum transfer {q = p — p') square = m\ + m^/ — 2mBmB'W. In the 
decay w ranges from w = 1, corresponding to zero recoil of the final baryon, to a 
maximum value given hy w = tUmax = {it^'b + "7^/)/(2?7iBmB/) which depends on 
the transition. 

Neglecting lepton masses, we have for the diff'erential decay rates from trans- 
vers ely ( Tr) and longitudinally (Tl) polarized VF's (the total width is F = F^;, + 
Ft) H 



dTr 
dw 

dTL 
dw 



Gl\Vcb\' 
— n — rriR 

Gl\Vcb\ ' 
247r3 



-m 



>y'w^-lq^{iw - l)\Fiiw)\'' + {w + l)\Gi{w)\^] 
,Vw^-l{{w - l)\T^{w)\'^ + {w + l)\T^{w)\^^ 



(niB ± mB')F^'^{w) + {1 ±w) I mB'F2'^{w) + mBFp^{w) 



Fj" ^Fj{w), Ff^Gjiw), j 



1,2,3 



(6.5) 



One can also evaluate the polar angle distribution 223l |: 

d^r 



du; d cos 6 



3 f^^^dljA (i^2a'cose + a"cos^9} (6.6) 
8 \ dw dw y L J 



where 9 is the angle between k' and p' measured in the off-shell W rest frame, 
and a' and a" are asymmetry parameters given by 



a 



a 



Gl\V^\^ 3 qHw'-l)F,{w)G,{w) 



67r3 drT/du; + 2drL/di(; 

dTr/dw -2dTL/dw 
dTx/dw + 2dTL/dw 



(6.7) 
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These asymmetry parameters are functions of the velocity transfer w and on av- 
eraging over w the numerators and denominators are integrated separately and 
thus we have 



(«') = - 



67r3 Ft 1 + IRlit 



{a") 
Rl/t 



1-2R 



L/T 



1 + 2Rl/t 
Ft 



(6.8) 



6.3 Form factors 



To obtain the form factors we have to evaluate the matrix elements 



B ,r p 



*^(0)7;.(I-75)*'(0) 



B,r p 



which in our model are given by 

->/ ^/ 



'Enip) Eb'{p') 



ruB 



m B' 



B',r' p' 



*'=(0)7^(7-75)*'(0) 



B,r p 



NR 



(6.9) 



(6.10) 



where the suffix "NR" denotes our nonrelativistic states and the factors E/m 
take into account the different normalization. We shall work in the initial baryon 
rest frame so that p = 0, p' = —q, and take q in the positive z direction. Further- 
more we shall use the spectator approximation. Having all this in mind we have 
in momentum space 



lEs'i-q) 



= 2 



B',r' (0)7^(1 -75)*''(0) 



B,rO 



NR 



EB'{-q) 



rriB' 



EE 



1 1 
2' 2 



, Sh 



Sl,S2 - Si,S2 Sh, 



1 1 

2' 2 



S2,r- S2,r 



X 



r' — r + S2,r — S2,r' 



Q> \,Shr' -r + si, S2 -si,r' -r + (^S'f,, ^, 
/"jS j3 [ ^b' rnh2+mq ^ ^ nih^ ^.\* ^b \ 

J d\id\2 [^chM^^ W^^'^' W^^V ^^'»2(«1'92) 

u';>-r+si{qi - q)ltj,{I -lb)u\^{qi) (6.11) 



Eh{Qi) V Ec{qi - q] 

where ^f^^ (^1,^*2) (^f/lj^^i' ^2)) is the Fourier transform of the coordinate space 

wave function ^^^^^ (''i' ''Sj ''12) (^^/^2('^i5 '''27 ^^12)) with ^1, 5*2 being the conjugate 
momenta to the space variables ri, r2- The factor of two comes from the fact 
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that: i) for 5c— baryon decays, the charm quark resulting from the b ^ c transition 
could be either the particle 1 or the particle 2 in the final cc baryon, while ii) for 
66— baryon decays, there exist two equal contributions resulting for the decay of 
each of the two bottom quarks of the initial baryon. 

The actual calculation is done in coordinate space where we have 



'Es'i-q) 



TUB' 



B , r 



B,r O) 

I NR 



' EB'{-q) 
niB' 



Si,S2 - Si,S2 ] { Sh, ^, ^ 



S2,r - S2,r 



r' -r + si,S2-si,r' -r + S2^ (^S'f,,^,^ 



r' — r + S2,r — S2,r' 



X / (iVi(iV2*fL(n,r2,ri2) e 



o A T ' ^ ^ c> 



nib 



nir. 



Eb{l ) Y Ec{l - q) 

(6.12) 



where I = —i Vi represents an internal momentum which is much smaller than 
the heavy quark masses m^, m-c. On the other hand |^| can be larg^. Thus, to 
evaluat e the above expression we have made use of an expansion in /, introduced 
in Ref. 144i |. where second order terms in / are neglected, while all orders in |g | are 
kept. For instance Ec{l — q) is approximated by Ec{l — q) ~ Ec{q) x {l — l-q / El{q)) 
with Ec{q) = y'rnf+if^. 

The three vector and three axial form factors can be extracted from the set of 
equation^ 



B',1/2 -q ^'(0)71^^(0) B,-l/2d 



' EB'{-q ) + ms' 



5', 1/2 -q 



2mB' EB'{-q)+mB 
r (0)73*''(0) B, 1/2 O) 

Fi{\q\) 



'EB'[-q) + rriB^ 



+ 



E3{\q\ 



2mB' ' \EB'{—q)+mB' ruB' 

B', 1/2 - f r(0)7„*'(0) B, 1/2 O'; 



2m B' 



niB' 



one has \q \ — {m% — m%,) /2mB which is ~ ms/'i for the transitions under study. 
^Remember g* is in the z direction. Notice also that, for p — Q, w is just a function of \q\. 
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for the vector form factors and 

.b 



B',1/2 



B',1/2 



B',1/2 -q 



r(0)7i75^''(0) 5,-1/2 



'EB'{-q) + niB' 



i-Gii\m 



' 2mB' 

^'(0)7375^''(0) 5,1/2 



lEB'{-q) + iriB' 



2mB' 

F(0)7o75*''(0) B, 1/2 



^1(191) + 



\q?Gm) 



lEB'{-q) + ruB' 



2m B' 



EB'{-q) + ruB' 



mB'{EB'{-q) + ruB') 



^EEtllGM)] (6.14) 



ruB' 



for the axial ones. All the left hand side terms can be evaluated using Eg . (16.121) 
with the approximation mentioned above. 

For each transition there are only two different coordinate space integrals from 
which all different matrix elements can be evaluated. Those integrals are 



■ q-ri 



'^cUn,r2,ri2) 



■ q-ri 



^•9 ^f/^2(n,r2,rl2)(6.15) 



In appendix lH.ll we relate the form factors to the integrals Z^'^{\q\) and IC^'^{\q\) 
for the different Sh, S'^ combinations, while in appendix IH.2I we give the actual 
expressions we use to evaluate those integrals. 



6.3.1 Current conservation 

In the limit = rric and for B' = B (and thus Sh = S'f^) vector current conser- 
vation provides a relation among the vector F2 and -F3 form factors, namely 



F2iw) = Fsiw) 



(6.16) 



On the other hand the matrix element of the zeroth component of the vector 
current evaluated at = 1 just counts the number of heavy quarks so that we 
should have 

Fi(l)+F2(l)+F3(l) = 2 (6.17) 
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In this limiting situation the integrals I^^{\q\) and }C^^{\q\) are related bjH 

ic^^m) = ^^^^^^i^^m (6.18) 

Besides one has that J^^(O) = 1. 

Using now the relations in Eq. ()H.4p in appendix IH.ll we see that our model 
satisfies the constraint in Eq. (l6.17p exactly. On the other hand we violate current 
conservation . For instance, and again using the relations in Eq. ()H.4p . we obtain 
for w = 1 

F2(1) = F3(1) + 2(1-^) (6.19) 

which shows that current conservation is violated by a term proportional to 
the binding energy of the baryon divided by the sum of the masses of its con- 
stituents. This violation disappears in the infinite heavy quark mass limit. This 
deficiency is shared by the relativistic calculation of Ref. 220l | and it is avoided 



in two other |l79l . H by the neglect of binding effect!0. Improvements on vec- 



tor current conservation would require at minimum the introduction of two-body 
currents [Z^, going thus beyond the spectator approximation, that we have not 
considered in this analysis. 

Note also that, for transitions that do not conserve the spin of the heavy quark 
subsystem Sh (i.e. E'^^ E'^Jui) we have in the nib = f^c limit and at zero recoil 
that 

Fi(l)+F2(l) + F3(l) =0 (6.20) 
due to the orthogonality of the initial and final baryon wave-functions. 



6.4 Results 

In Figs. 16.11 16.21 we show the form factors for the different transitions evaluated 
with the ALl potential. Variations when using a different potential are at the 
level a few per cent at most. The results for doubly heavy H decays are almost 
identical to the corresponding ones for doubly heavy decays. The fact that we 
have two heavy quarks and that the light one acts as a spectator makes the results 
almost independent of the light quark mass. 

In Figs. [6^ 16.41 we show now our results for the differential dTx/dw, dVi/dw 
and dT/dw decay widths evaluated with the ALl and BHAD potentials. The dif- 
ferences between the results obtained with the two inter-quark interactions could 

■^One just has to i ntegr ate by parts in the JC'^^{\q\) expression. 

■*!£ we look at Ref. [220l ] and consider for instance transitions where the initial and final baryons 
have total heavy quark spin equal to 0, we see that vector current conservation in the equal mass 
case wo uld r equire the h- (ui) form factor to vanish. This is not accomplished within that model. 
In Ref. J22^ the currents are constructed at the diquark level and the vector part is conserved 
thanks to the neglect of binding effects with the light quark. Finally, the calculation in Ref. [l79| ] 
avoids this problem by using the infinite heavy quark mass limit and thus cancelling binding 
effects. 
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reach 30% for some transitions and for some regions of w. As a consequence of 
the apparent SU(3) symmetry in the form factors we also find that the results for 
doubly heavy H and 17 decays are very close to each other. This apparent SU(3) 
symmetry goes over to the integrated decay widths and asymmetry parameters. 

In Table 16.11 we give our results for the semileptonic decay width (transverse 
Tt, longitudinal Tl and total T) for the different processes under study. Our 
central values have been evaluated using the ALl potential while the errors show 
the variations when changing the interaction. The biggest variations appear for 
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Figure 6.2: Same as Fig. 16.11 for doubly heavy Q{J = 1/2) baryons decays. 



the BHAD potential for which one obtains results which are larger by 7 ~ 12%. 
In Table 16.21 we compare our results to the ones calculated in different models. 
For that purpose we need a value for \Vcb\ for which we ta ke IV ^^I = 0.0413. Our 



results are in reasonable agreement with the ones in Ref. 220f | where they use a 
relativistic quark model evaluate d in the quark-diquark approximation, and with 
the r(Hbc — > Sec) value of Ref. |l79l ] obtained using HQET. The value for the 
latter width but now evaluated in the relativistic three-quark model calculation 
of Ref. 22l| is much smaller than in any other calculation. On the other hand in 
Ref. [2221 ]. where they use the Bethe-Salpeter equation applied to a quark-diquark 
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Figure 6.3: dT/dw, dVi/dw and dTT/dw semileptonic decay widths in units 
of \Vcb\^ • 10~^^ GeV, for doubly H(J = 1/2) baryons decays. Solid line, long- 
dashed line and short-dashed line: dV/dw, dTi/dw and dTx/dw evaluated with 
the ALl potential; dotted line, long-dashed dotted line and short-dashed dotted 
line: dT/dw, dT^/dw and dTx/dw evaluated with the BHAD potential. 



system, they obtain much larger results for all transitions. 

In Table 16.31 we compile our results for the average angular asymmetries a' 
and a" , as well as the Rl/t ratio, introduced in Eq. (l6.8p . The central values have 
been obtained with the ALl potential. Being all quantities ratios the variation 
when changing the inter-quark interaction are in most cases small. 
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;ure 6.4: Same as Fig. 16.31 for doubly heavy Q{J = 1/2) baryons decays. 
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Table 6.1: Semileptonic decay widths in units of iKbP • 10"^^ GeV. Ft and F^ 
stand for the transverse and longitudinal contributions to the width F. The central 
values have been obtained with the ALl potential. I stands for a light charged 
lepton, I = e, ji. 
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Table 6.2: Semileptonic decay widths in units of 10 GeV. We have used a value 
1 14;, I = 0.0413. / stands for a light charged lepton, I = e, fx. 
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Table 6.3: Averaged values of the asymmetry parameters a' and a" evaluated as 
indicated in Eq. ()6.8p . We also show the ratio Rl/t = ^l/^t- The central values 
have been obtained with the ALl potential. / stands for a light charged lepton, 
I = e, fi. 



Chapter 7 



Strong pionic decay of heavy 
baryons 



7.1 Introduction 

In this chapter we shah evaluate strong widths for the Sc — > Ac vr, S* ^ Ac vr and 
H* ^ He TT decays. 

Last decade has seen a great progress on charmed-baryon physics and now 
the ground state baryons with a c quark, with the exception of the 0,*, are well 
estabhshed [usj], and we hav e experira ental informati on on the stro ng o ne-pion 



decay widths for the ^i, 113, E* HH, and H* |222|, H. With 



very httle kinetic energy available in the final state these reactions should be well 
described in a nonrelativistic approach. Although they hav e been aii alyzed before 
in the framework of the constituent quark model (CQM) 230 . l23ll |. no attempt 
was made there to evaluate the full matrix elements. While there have been 
dynamical calculations in other models (see references below), to our knowledge, 
ours is the first fully dynamical calculation within a nonrelativistic approach. 
In our calculation we will use the HQS~constrained wave functions evaluated in 
Ref. [20ft] using the same different interquark interactions that have been used 
in this work, and whose goodness have already been te sted in the study of the 
semileptonic A^ — > Ac and H;, — > He decays in Ref. 80, Il44i |. Again, the use of 
different quark-quark potentials will allow us to obtain theoretical uncertainties 
on the widths due to the quark-quark interaction. The pion emission amplitude 
will be obtained in a spectator model (one-quark pion emission) with the use of 
partial conservation of axial current hypothesis (PCAC) as was done, in the meson 
sector, in the previous study of the strong B*Btt and D*D7r couplings in section 



These reactions, and similar ones, have also been addresse d in heavy hadron 
chiral perturbation theory (HHCPT) |23l|, H H H H, in QCDSR |23^ . 



2371] and within relativistic quark models like the l ight-front quark model (LFQM) 23 
and the relativistic three-quark model (RTQM) SIm^]. 
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7.2 One-pion Strong decay width 

The width for the reaction B ^ B' tt is given by 
1 f d^PB' d^Pn 



2mB J {27ry^2EB'{PB') {2TTr2E^{P^) ' 



X 

2Jb 

where Pb = {mB;0), Pb' = {Eb>{Pb')-,Pb') and P^- = (£7^(^71), -Ptt) are the four- 
momenta of the particles involved. Jb is the spin of the B baryon and s and s' are 
the third component of the spin of the B and B' baryons in their respective center 
of mass systems. Finally J^q^)^{Pb , Pb') is the pion emission amplitude, that will 
be obtained through the non-pole part of the matrix element of the divergence of 
the axial current similarly to the case for mesons in section 13. 5[ 
The quantity 

^j^EEKB;i(^^'^^')f (7.2) 

has to be a Lorentz invariant. Denoting that invariant as \}A BB'nl'^ we can perform 
the integrations to give 

r = IMsB'.f (7.3) 

with \q\ the modulus of the three- momentum of the B' baryon or the vr meson 
given by 

^ X^^ {ml, 171%,, ml) 

2mB ■ 
Throughout the calculation we shall use physical masses taken from Ref. [ll5l ] . 



7.3 Description of baryon states 



We use a similar description for baryons as the one used in chapters \5\ [H but in 
this case we shall work in momentum space. The Fock space representation of a 
baryon B with three- momentum P and spin projection s in the baryon center of 
mass is 



B,sP 



NR 



{Qi,Q2. 



(2vr)3 ^2Ef,{pi)2Ef,{p2)2Ej,ip3 



ai pi 



M 



'^P + Qi 



"2 P2 



"3 P3 



^3 p 
M 



Qi 



M 
-Q2 



(7.5) 
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2 





1+ 


ssc 
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-2 


3 + 
2 





1+ 


ssc 



Table 7.1: Summary of the quantum numbers of ground state baryons containing 
a single heavy c quark. I, and S]^ are the isospin, and the spin-parity of the 
light degrees of freedom and S, are the strangeness and the spin-parity of the 
baryon. In the last column I denotes a light quark of flavor u or d. 

ai, 02 and 03 represent the quantum numbers of spin (s), flavor (f) and color 
(c), a = {s,f,c), of the three quarks, while {Ef^{pi), pi), {Ef^{p2), P2) and 
{Ef^^Ps), P3) are their respective four-momenta, is the mass of the quark 
with flavor /, and M = nif-^ + "^/2 + "^/a- choose the third quark to be the c 
quark while the first two will be the light ones. 
The normalization of the quark states is 

{a'p'\ap) = {2Trf 2E{p) 5\p' - p) (7.6) 

Besides, Vji^oaiOs ( Qij Q2 ) is the momentum space wave function for the internal 
motion, being Qi and Q2 the momenta conjugate to the positions ri and r2 of the 
two light quarks with respect to the heavy one (See Fig. 17. ip . This wave function 
is antisymmetric under the simultaneous exchange ai < — > Ci2,Qi < — > Q2; being 
also antisymmetric under an overall exchange of the color degrees of freedom. It 
is normalized such that 

jd^Qjd^Q2 (CiU(^i'^2))*V'SL«3(^i'^2) = 5.',. (7.7) 

01,02,03 

and, thus, the normalization of our baryon states is 

[b,s'P'\B,sP = 5s', s (2vr)3 d^P' - P) (7.8) 

For the particular case of ground state Ac, T,c, S*, He and H* we can assume 
the orbital angular momentum to be zero. We will also take advantage of HQS 
and assume the light-degrees of freedom quantum numbers are well defined (See 
Table 17. ip . In that case we hav^ 

^We only give the wave function for the baryons involved in tt^ emission. Wave functions for 
other isospin states of the same baryons are easily constructed. 
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'5, = 1/ 



(7.9) 



Here £ciC2C3 is the fully antisymmetric tensor on color indices being Ec^c^czI ^f^^- 
the antisymmetric color wave function and the (/>j' Qi, Q2 ); with 5; the total 

spin of the light degrees of freedom, are the Fourier transform of the corresponding 
coordinate space wave functions (i;^^^ jgl^i ''^2)5 see below). Their dependence on 
momenta is through \Q\\, \Q'2\ and Q\ ■ Q2 alone, and they are symmetric under 
the simultaneous exchange /i < — > f2,Qi < — > Q2- Their normalization is given 
by 



^Uj,SQi.Q2) 



(7.10) 



7.4 Coordinate space wave functions 



Wave functions for baryons with a heavy quark are quite similar to the ones 
presented in chapter [5l but there are som e di fferences. In this section we briefly 
explain how they are constructed in Ref. 200l |. 



7.4. COORDINATE SPACE WAVE FUNCTIONS 
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7.4.1 Intrinsic Hamiltonian 



Working in the heavy quark frame {R,fi,f2) (See Fig. 17. ip . where R, fi and r2 
are the center of mass position in the laboratory frame and the relative positions 
of the light qi and q2 quarks with respect to the heavy one h, the center of mass 
motion can be separated from the intrinsic Hamiltonian, 



ft , TT 

int 



2M 

Vi • V2 



Hint = M+^Hf + Vq^^q^ (fi - fa, spin) 



nih 

j=l,2 



V? 

i?f = -l7^ + yh,g,{rj,spin), i = l,2 (7.11) 



with M = ruq^ + ruq^ + ruh, fJ.j = {l/rug^ + l/rrih) ^ and Vj = d/d^^. , j = 1, 2. 
Apart from M, Hint consists of the sum of two single particle Hamiltonians {H^^), 
describing the dynamics of the light quarks in the mean field created by the heavy 
quark, plus the light-light interaction term that includes the Hughes-Eckart term 
(Vi-Vs). 

For the inter-quark interaction V, the potentials presented in chapter [2] were 
used. 



122 



CHAPTER 7. STRONG PIONIC DECAY OF HEAVY BARYONS 



7.4.2 Variational coordinate— space wave function and its Fourier 
transform 

To s olve the three-body problem the following variational ansazt was used in 



Ref. [200f | for the spatial part of the wave function 



^l,,,hi^^^^2) = N^h,An)4>h,,{r2)F'''{rr2) (7.12) 

with rj = |r^|, ri2 = \'ri — f2\- is a normalization constant, (j)hq is the s-wave 
ground state solution ( ijjhq. ) of the single particle Hamiltonian ( H'^^ ) corrected 
at large distances in the form 

{rj) = (1 + aj Tj) iphq, (rj), j = 1,2 (7.13) 

and finally F^^q^{ri2) is a Jastrow correlation function in the relative distance of 
the two light quarks and was constructed as a sum of gaussians 

4 

F^^{ri2) = J2 «i e"^^' ('^12+'^^)', ai = 1 (7.14) 
i=i 

The variational parameters for the different states obtained with the different 
potentials can be found in the appendix of Ref. 200]. The iphqj functions can be 
easily obtained through a Numerov algorithm. 

In momentum space we will have 

^S,,,,,(Qi,Q2) = 7^3/ [ d\2 e^^^^"^^ e'"^--'^- N^hqAn)cl>hq,{r2)F''(.ru) 



(27r)3 

J d'P iV0,,,(|P|) ^Hq,i\Ql +Q2- P\) F^KIQl - ^1) 

(7.15) 



being 4>hq^ , (j)hq2 and F^^ the Fourier transforms of (phq^ , 4'hq2 and F^' respectively. 
They depend solely on the modulus of the momentum. Clearly 

<t>qUhiRQl^RQ2) = <f>qUHiQl,Q2) (7.16) 

with it! a rotation. For the purpose of evaluation we can take Qi in the z direction 
and Q2 in the OXZ plane. 

7.5 Results 

7.5.1 Sc ^ AcTT decay 

Let us do the A+ 7r+ case. The matrix element of the divergence q^J'^^{0) 

of the axial current determines the vr"'" emission amplitude as 



' non-pole 



^ffs++A+^+ s'iP' ) 75 u^++ siP) (7-17) 
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9i: + + A+n + 


r(Er - K^^) 


r(E: ^ A^Tv^') 


r(EL'^A+^-) 




[MeV] 


[MeV] 


[MeV] 


This work 


21.73 ±0.32 ±0.08 


2.41 ±0.07 ±0.02 


2.79 ±0.08 ±0.02 


2.37 ±0.07 ±0.02 


Experiment 




2.3 ±0.2 ±0.3 [224] 
2.05toil ± 0-38 [2261 


< 4.6 (CL=90%) [225] 


2.5 ±0.2 ±0.3 [224] 
l-55to.tl ± 0.38 [226] 


Theory 










CQM 




1.31 ±0.04 [230] 
2.0251^:1^8^ [231J 


1.31 ±0.04 [230] 


1.31 ± 0.04 [230] 
l.939tl;llt [231J 


HHCPT 


22, 29.3 [232] 


2.47, 4.38 [232] 


2.85, 5.06 [232] 


2.45, 4.35 [232] 






2.5 [233] 


3.2 [233] 


2.4 [233] 
1.94 ±0.57 [234] 


LFQM 




1.64 [238] 


1.70 [238] 


1.57 [238] 


RTQM 




2.85 ±0.19 [239. 240] 


3.63 ±0.27 [239. 240] 


2.65 ±0.19 [239, 240] 



Table 7.2: Coupling constant S's++A+7r+ ^'^d total widths r(i;^+ A+vr^), 
r(S+ A+vr*^) and r(S[! — > A+tt") (See text for details). Experimental data and 
different theoretical calculations are also shown. 



where the coupling constant gj,++j^+^^ , in analogy to the pion coupling to nucleons 

and nucleon resonances, has been chosen to be dimensionless, and Uj.++ g{P), 

Uj^+ gi{P') are Dirac spinors normalized to twice the energy {u^u = 2E). The 
width is given by 

r(E++ ^ A+vr+) = J|L- 4,,^,^, ( (M^.. - M^,f - m^) (7.18) 

with \ q\ the modulus of the final baryon or pion three-momentum. From Eq. (j7.17p . 
taking P = 0, P' = — ]g ] A; in the z direction, s = s' = 1/2, and taking into account 
the different normalization of our nonrelativistic states, we have 



5s++A+7r^ 



/tt ]( 



X ((M,.. - E,^m ^ + 1,1 3) (7.19) 



with 



^v+;A',^ = iVi?(^- 1/2 4^(0) l/2 0)_ , (7.20) 



/ NR, non-pole 



1/2 1/2 

The ^ are easily evaluated using one-body current operators and their 

expressions can be found in appendix HI 

In Table 17.21 we present the results for g^++^+_^j^ and the widths r(S++ 
A+^+), r(S+ ^ A+^°) and r(SO ^ A+vf-). To get the values for r(S+ ^ A+vr^) 
and r(S2 ^ A+7r~) we use 5's++A+7r+ ^^"^ rnake the appropriate mass changes in 
the rest of factors in Eq. ()7.18| . 



124 



CHAPTER 7. STRONG PIONIC DECAY OF HEAVY BARYONS 



Our results show two types of errors. The second one results from the Monte 
Carlo evaluation of the integrals needed to obtain the 5'2++A+7r+ coupling constant. 
The first one is the usual theoretical uncertainty coming from the use of different 
potentials. The results are in very good agr eement with the experimental data by 
the CLEO Collaboration in Refs. liii,^. The value for r(S++ A+7r+) also 
agrees with the experimental data by the FOCUS Collaboration in Ref. [22g]. The 
agreement with FOCUS data is not good for the r(E|? — > A^vr") case, although 
our result is still within experimental errors. Our results show variations as large 
as ~17% between different charge configurations. This is due to the little kinetic 
energy available in the final state that makes the widths very sensitive to the 
precise masses of the hadrons involved. In this respect there is a new precise 
determination of the mass by the BABAR collaboration = 2286.46 it 0.14 
MeV/c^ [2411], which is ro ughly 1.5 MeV/c^ above the value quoted by the Particle 



Data Group in Ref. 115l |. With this new value our calculated widths would get 
reduced by 9%. This reduction comes from phase space factors while the coupling 
fl's++A+7r+ changes only at the level of 0.1%. 



As for the other theoretical determinations, the CQM calculation in Ref. j230l | 
uses exact s < — > c symmetry to relate the Sc — > AcVr decay to the non-charmed 
S* — > Avr analogue decay. Their results are sm aller than the experimental data 
obtained by the CLEO Collaboration. In Ref. 23ll | a unique coupling constant 
is fixed in order to reproduce all experimental information on S* Acvr widths. 
That coupling is latter used to predict the Sc — > AcVr widths. This coupling suffers 
from large uncertainties and thus the theoretical e rrors on the predicted widths are 
also very large. In the HHCPT calculation of Ref. |232| a simple CQM argument is 
used in order to obtain the unknown coupling constant in the HHCPT Lagrangian. 
Furthermore the authors allow for a renormalization of the axial coupling for 
light quarks. The largest of the two values quoted corresponds to the case where 
that coupling is unrenormalized and then is given by = 1. The smaller number 
quoted corresponds to the use of a renormalized value of g^ = 0.75. The case 
^nd = ]^ is ii^Q Qjig that compares with our calculation. Their results for the 
widths almost double ours and are not in agreement with experiment. Their 
simple determination of the coupling consta nt ca n not be correct. The values 
obtained in the HHCPT calculation of Ref. H are closer to our results and 



experimental data. There the authors determine the needed coupling constants 
from the analysis of analogue decays involving non-charm ed b aryons. In Ref. [234 ], 
also within the HHCPT approach, and similarly to Ref. j23ll ], the authors fix the 
unknown coupling in the Lagrangian using the experimental data for the S* — > AcTT 
decays. From there the y predict the r(S|? A^vr") obtaining a value very close 



to the one in Ref. 2311 ]. and that suffers als o from large u ncertainties. The two 



relativistic quark model calculations of Refs. 238 . 2391 . 240 ] give results that differ 
by almost a fac tor of tw o. Our results are closer to the ones obtained within the 
RTQM of Ref. 
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r(E*++ ^ Aj7r+) 










[MeV] 


[MeV] 


[MeV] 


This work 


36.20 ±0.75 ±0.13 


17.52 ±0.74 ±0.12 


17.31 ± 0.73 ± 0.12 


16.90 ±0.71 ±0.12 


Experiment 




14.1+i;^ ± 1.4 [227] 


< 17 (CL=90%) [225] 


16.6^i;'^ ± 1.4 [227] 


Theory 










QCDSR 


13.8 ^ 24.2 [236] 










32.5 ±2.1 ±6.9 [237] 








CQM 




20 [230] 


20 [230] 


20 [230] 


HHCPT 




25 [233] 


25 [2331 


25 [233] 


LFQM 




12.84 [238] 




12.40 [238] 


RTQM 




21.99 ±0.87 [239. 2401 




21.21 ±0.81 [239, 240] 



Table 7.3: Coupling constant 9s*++A+7r+ ™^ total widths r(I]*++ A+7r+), 
r(E*+ A^TT^) and r(S*'^ A^vr^). Experimental data and different the- 
oretical calculations are also shown. Note that in order to compare with our 
definition oi_2r,* ++ a'^ it+ have multiplied the coupling constants evaluated in 
Refs. |23i,[23jby 2M^+/f^. 



7.5.2 S* ^ AcTT decay 

Let us analyze the case with a tt~^ in the final state, — > A+tt^. Similarly to 
the Sc decay before we now have 



.5s*++A+,^^ 



f Jl"(0) sP 



non~pole 



(7.21) 



where we have introduced the dimensionless coupling constant g^* 



(P) is a Rarita-Schwinger spinor normalized to twice the energy {u'^^Ui, 



and 



-2E). The width is given by 



^ ^ ' ^ 247rM2 



4*++ A 



^((M^*+++M^+)2-m2) (7.22) 



Taking again P = 0, P' = —\q \ k in the z direction, and s = s' = 1/2 we obtain 
from Eq. (l73Tl) 



fe++A+7r- 



V3 2M^+v^^s^-^+^^Ac+(9: 



with 



=^kK+, 1/2 4"^(0) 1/20) (7.24) 

Ac, \ / Af_R, non-pole 
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The expressions for ^^{^lY^^. (/i = 0,3) can be found in appendix HI 
Results for g^,++^+^+ and the total widths r(S*++ A+vr+), r(S*+ A+tt' 
and r(S*~ A^vr^) appear in Table [731 Our value for the latter two are 
obtained with the use of 92*++A+7r+ ^^'^ with the appropriate mass changes in 
the rest of factors in Eq. ()7.22p . Our central value for r(S* — > A^vr"'") is 
above the central val ue of the latest experimental determination by the CLEO 



Collaboration in Ref. [2271 ] . For some of the potentials used, API and AP2, the 
results obtained are within experimental errors. The central value for T(Ti*~^ 
A^vr*^) is slightly above the upp er experimental bound determined also by the 
CLEO Collaboration in Ref. [225l], but again, we obtain results which are below the 
experimental bound using the API and AP2 potentials. As for r(S*~ A^tt^) 
we get a nice agreement with experiment. Our results for the different charge 
configurations differ by 4% at most . With the new value for given by the 



BABAR Collaboration in Ref. 2411] they would get reduced by 3%. Our results 
are globally in better agreement with experiment than the ones obtained by other 
theoreti cal c alculations d with perhaps the exception of the QCDSR calculation 



of Ref. [237 ]. 



7.5.3 S* ScTT decay 

Once more we analyze the case with a tt"*" in the final state, H*+ ^ H°7r+. What 
we obtain is 

' JjY ^ I non-pole 



A]::V^,^4P,p') = -{E^,, s'p'iq^JiimiK^, sP^ 

q.uso AP')u'i,+ (P) (7.25) 



where again we have introduced a dimensionless couphng o^*+-^o +• The width is 
given as 

T(E*+ ^ H07r+) = — ( (M-,+ + M^of - ml) (7.26) 

Taking now P = 0, P' = —\q \ k in the z direction, and s = s' = 1/2, o„*+„o , is 
evaluated from Eq. (l7.25p to be 



/tV2 \^\^2M^,+ {E^oiq) + M=o) 

X ((M..+ - E^o{q)) ^ + \q\ A'J,\%1 3) (7.27) 



^We did not show those cases were data on E* Acvr widths were used to fit parameters of 
the models. 
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9^.+^0,+ r(H;+-.H»^+) r(H;+-,H+^») 

^ [MeV] [MeV] 


r{H;"^S+,r-) r{3*"^S^7r") 
[MeVI [MeVI 


This work 


-28.83 ± 0.50 ± 0.10 1.84 ± 0.06 ± 0.01 1.34 ± 0.04 ± 0.01 


2.07 ± 0.07 ± 0.01 0.956 ± 0.030 ± 0.007 


Theory 
LFQM 
RTQM 


1.12 [2381 0.69 [2381 
1.78 ± 0.33 p39, 2401 1.26 ± 0.17 [539, 2401 


1.16 12381 0.72 [2381 
2.11 ± 0.297539, 2401 1.01 ± 0.157539, 2401 




r(H;+ ^H0,r++H+7r0) 
[MoVI 


r(H;o^H+.-+HO.O) 

[MeVI 


This work 


3.18 ± 0.10 ± 0.01 


3.03 ± 0.10 ± 0.01 


Experiment 


< 3.1 (CL=90%) [2281 


< 5.5 (CL=90%) [229J 


Theory 
CQM 

HHCPT 


< 2.3 ± 0.1 [2301 
1.191- 3.97172311 
2.44 ± 0.85 [2341 


< 2.3 ± 0.1 [2301 
1.230 - 4.074 12311 
2.51 ± 0.88 [2341 


LFQM 
RTQM 


1.81 [238^^ 
3.04 ± 0.501239. 2401 


1.88 [238| 
3.12 ± 0.337239. 2401 



Table 7.4: Values for the coupling g„*+„Q^_^^ and decay widths r(H*+ — > E'^tt'^), 
r(H*+ ^ H+ttO), r(H*o ^ H+vr") and r(H*° ^ HOyr^). Experimental upper 
bounds for the total H*+ and 'E*^ widths, and different theoretical calculations 
are also shown. 



with 



=^^(hO, 1/2 -\q\k\ JZ{0)\E: + , 1/2 o) (7.28) 

— c "ci A* \ I NR, non-pole 



which expressions can be found in appendix HI 

iiin will, iir, 1 I ^ 7^/-*+ . '='+^0 



Results for the coupling (/„*+^o^+, the widths r(H*+ — > H^7r+), r(H*+ EJtt' 
r(H*o ^ H+TT") and r(Hfo ^ E^^tt^), and the total widths r(H*+ ^ E^^tt+ + 
H+vrO) and r(H*o ^ E^^tt^ + E+n') appear in Table El Our values for r(H*+ ^ 
H^vr'^), r(H*° — > H+vr^) and T{E*^ ^c'^^) obtained with the use of 
, /(M_+ + M=o), and with the appropriate mass changes in the rest of 
factors in Eq. (17:26]) . For r(H*+ ^ H+7r°) and r(H*° H°7r°) an extra 1/2 
isospin factor should be included. Our results for r(H*+ — > E^tt^ + Etn^ ) ar e 



slightly above the experimental bound obtained by the CLEO Collaboration [228l |. 
As for the S* decay case above, the API and AP2 potentials gives results closer 
to experiment. For T{E*^ '^t'^~ + H^tt^ ) ou r result is well below the CLEO 



Collaboration experimental bound in Ref. 229f |. Isospin breaking due to mass 
effects is clearly seen when comparing the predictions for r(H*+ E^7:~) and 
r(H*+ — > Efn^). One finds a factor 1.4 difference when a factor of two would be 
expected from isospin symmetry. Again, the fact that there is little phase space 
available makes the results very sensitive to the actual mass values. Compared 
to other theoretical calcul ations ou r results agree nicely with the ones obtained 
within the RTQM of Refs. [23^, [m^], while they are larger than most other deter- 
minations. 



Chapter 8 

Conclusions 



This work has been devoted to the study of different properties of hadrons with 
one and two heavy quarks c and/or b. All calculations have been done in the 
framework of a nonrelativistic constituent quark model. In order to check the 
sensitivity of our results to the inter-quark interaction we have used five different 
quark-quark potentials. The spread in results gives an estimation of theoretical 
uncertainties and so has been quoted. Most observables studied change only at the 
level of a few per cent when changing the interaction potential. Another source 
of theoretical uncertainty is the use of nonrelativistic kinematics in the evaluation 
of the orbital wave functions and the construction of our nonrelativistic states, 
although we think that a good part of these relativistic effects are contained in 
an effective way in the parameters of the quark-quark potentials, which were 
adjusted in the original works to reproduce the experimental spectra of mesons. 
Mesons wave functions have been evaluated using a Numerov algorithm, while 
for baryons the three body problem has been solved using a variational approach 
with ansatz wave functions that include simplifications that arise from HQS and 
HQSS. To compute hadron decays we have relied on the impulse approximation, 
which induces small axial and vector current conservation violations due to binding 
effects. These binding effects, though they affect some of the weak form factors 
analyzed, have little influence on the decay widths. 

• In chapter [3] we have evaluated leptonic decay constants of D, D* , B and B* 
mesons, and form factors and decay widths for the semileptonic B ^ DlV and 
B D*lV. 

Our analysis of leptonic decay constants shows that in a nonrelativistic cal- 
culation the equality /yMy = fpMp is satisfied within 2%. This equality is 
expected in HQS in the limit where the heavy quark masses go to infinity. The 
nonrelativistic result suggests that the HQS infinite mass limit sets in already at 
the rric scale, something that is not supported by lattice data for D mesons where 
one finds deviations as large as 20% from the above equality. 

One also finds problems in the semileptonic B ^ D and B — > D* decays. We 
have seen how the /i„ (w) form factor of the B ^ D decay and the /lyij (^) form 
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factor of the B D* decay are not reliably calculated in the nonrelativistic quark 
model where one finds large deviations from the HQET relations of Eq. (I3.22p . We 
have tried to remedy this failure by evaluating the Isgur-Wise function from a 
form factor whose calculation in the quark model we trust: h^{w) for the B ^ D 
decay, and {w) for the B ^ D* decay. Those Isgur-Wise function were later 
used together with HQET constraints in Eq. (j3.22p to recalculate all other form 
factors. The two Isgur-Wise functions thus determined show an overall reasonable 
agreement with lattice data but in both cases the slope seems to be too small. 
This deficiency multiply its effects when one goes to larger w values and as a 
consequence the quantities F]:){w)\Vcb\ and Fj:)*{w)\Vcb\ go above experimental 
data for w values larger than 1.2, and for any reasonable value of \Vcb\- A failure 
of some kind is expected in a nonrelativistic calculation as w increases. For w = 1.2 
the three momenta of the final meson amounts to 66% of its mass and relativistic 
corrections in the wave function could start to be important. On the other hand 
we believe our results are sound at zero recoil {w = 1). That enables us to obtain 
\Vcb\ = 0.040 lb 0.006 from the B ^ D decay. This result is in perfect agreement 
with the value \Vcb\ = 0.040 ± 0.005 previously obtained using the same model in 



the analysis of the A;, — > Ac semileptonic decay in Ref [1441 ] , and in good agreement 
with a recent determination \Vcb\ = 0.0414 ± 0.0012 ± 0.0021 ± 0.0018 by the 
DELPHI Collaborati on |l51i1 . and with the average value \ Vcb\ = 0.0416 ± 0.0006 
favored by the PDG |l88l |. The experimental situation concerning the B — D* 
reaction is not so clear as different experiments give values for .Fn* (w) \ Vcb\ which 
are hardly compatible. From DELPHI Collaboration data [15ll | we would get 
\Vcb\ = 0.040 lb 0.003 in agreement with the above result. 

Finally we have made use of PC AC to evaluate the strong coupling constants 
gB*B-K and gD*DiT- Our results are larger than experimental data or the results 
provided by lattice and QCDSR calculations. In this case the final meson is 
nearly at rest and one would expect a nonrelativistic calculation to perform better. 
The main difference with the other observables analyzed is that here the two 
active quarks are the light ones. On the other hand the value for the ratio R = 
30*^(0) /g* ™ '^-^ agrees with the HQS prediction and with the one 
evaluated using a combination of lattice and experimental data, being also close 
to a QCDSR determination. 

Results of this chapter were published in Ref. |242l | 

• In chapter Owe have made a comprehensive and exhaustive study of exclusive 
semileptonic and nonleptonic two-meson decays of the Be meson. We have left 
out semileptonic processes involving a 6 — > u transition at the quark level to avoid 
known deficiencies both at high and low q'^ transfers (see section 14. ip . For simi- 
lar reasons we have only considered two-meson nonleptonic decay channels that 
include a, cc or B meson in the final state. Our model respects HQSS constraints 
in the infinite heavy quark mass limit but hints at sizable corrections away from 
that limit for some form factors. Unfortunately such corrections have not been 



worked out in perturbative QCD as they have been for heavy-light mesons [135l |. 
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Our results for the observables analyzed are in a general good agreement 
(whenever comparison is possible) with the results obtained within the quasi- 



16,1 Il6fil| 



potential approach to the relativistic quark model of Ebert et al. 

The branching ratios for the leptonic B~ — > cc and B~ — > B decays are also 
in reasonable agr e eme nt with the relativistic constituent quark model results of 



Ivanov et al. 162l . Il63l | . 



For the nonleptonic B^ 



ricMp 



and B^ J/^Mp two-meson decay 



channels with Mp = tt , p , K , K* , we find also reasonable agreement (b etter 
for the J/^ channel) with the Bethe-Salpeter calculation by El-Hady et al. |l7Cll ] 



and the light front calculation by Anisimov et al. |l77l|. w hile our results are 



a factor of two smaller than the ones by Ivanov et al. |l62l |. and Chang et al 



1671 . Il68l . Il69l |. the latter obtained within the nonrelativistic approach to the 



Bethe-Salpeter equation. For the two-meson decay channels with XcO) Xci, he, Xc2 



or ^'(3836) as the final cc meson and Mp = tt^ , , K~ , K* ~ ou r results are 
generally a factor of two smaller than the ones of Ivanov et al. |l62 |. whereas for 
some channels (xcO tt" , XcO P~ j hc'^~ ■> he p~ , XdLKT and P~ ) we find very good 



agreement with the results by Chang et al. [167 . IGSl . [l69.] . The disagreement 
with Ivanov et al. extends to the two-meson decay channels with final cc a nd D 
mesons. There we find good agreement with the results by El-Hady et al. 17Cll | 
and the ones obtained within the sum rules of QCD and nonrelativistic QCD by 
Kiselev fl74]. 

As for the two-meson nonleptonic decay channels B^ B Mp with Mp = 
TT^ , p~ , K~ , K*~ we are in a reasonable good agreement with the results by 
Anisimov et al. [l77|. F or th e case of a final B^ or B*f the agreement with the 



results by Ivanov et al. 1621 ] is very good. For the B^ 



B- 



Mf = vr°, /)°, K"", K 



'*0 



Mp c ase with 
and apart from the results by Ebert et al. [16a], we find 



no global agreement with other calculations, neither do they agree with each other. 

From the above comparison one sees that there are different models producing 
sometimes very different results for the same observables. Accurate experimental 
data is needed to shed light into this issue. 



Results in this chapter were published in Ref. 2431 ] 



• In chapter [S] we have computed wave functions, spectra and other static 
properties of doubly heavy baryons. In order to build our wave functions we have 
made use of the constraints imposed by the infinite heavy quark mass limit. In 
this limit the spin-spin interactions vanish and the total spin of the two heavy 
quarks is well defined. With this approximation we have used a simple variational 
approach, with Jastrow type orbital wave functions, to solve the involved three- 
body problem. We have also checked that our wave functions have the correct 
behavior in the infinite heavy quark mass limit. 

Among the static properties, our results for the masses are in very good agree- 
ment with previous results obtained with the same inter-quark interactions but 
within a more complicated Faddeev approach [75|. In some cases we even get 
lower, and thus better, masses. We have calculated charge densities (charge form 
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factors) finding that the corresponding mean square radii are again in good agree- 



ment with the Faddeev calculation of Ref . 75|] . We have also evaluated magnetic 
moments. Being the total orbital angular momentum of the baryon L = 0, the 
magnetic moments come from the spin contributions alone. With the exception 



of H^^, and we agree perfectly with the Faddeev calculation in Ref. 75l |. 
For the magnetic moments of H^^, and Jl^^, the discrepancies between the two 
calculation are very large. The origin might be attributed to the possible presence 
of a non-negligible = component in the wave functions of Ref. [75l] . In our 
case we have Sh = ^ which we think is a good approximation based on the infinite 
heavy quark mass limit. This assertion seems to b e corroborated by the results 
obtained in the relativistic calculation of Ref. [215^, at least for the and 
cases. 



Results in this chapter were published in [244]. 



• In chapter [6l and using the wave functions obtained in chapter [5l we have 
evaluated form factors, decay widths and angular asymmetry parameters for semilep- 
tonic decays involving the ground state of doubly heavy H and baryons. We 
have a small vector current violation by an amount given by the binding energy 
over the mass of the baryon that disappears in the infinite heavy quark mass limit. 

Our results for th e tot al decay widths are in reasonable agreement with the 



ones obtained in Ref. [220l | using a relativistic quark model in the quark-diq uark 
approximation, while they are much smaller than the ones obtained in Ref. H 
by means of the Bethe-Salpeter equation applied to a quark-diquark system. 

For the weak form factors the results exhibit an apparent SU(3) symmetry 
when going from S to baryons. This is due to the fact that we have two heavy 
quarks and the light one acts as a spectator in the weak transition. This apparent 
symmetry appears also in the decay widths and asymmetry parameters. On the 
other hand SU(3) violating effects are clearly visible in some static quantities like 
charge form factors and radii and magnetic moments, that depend strongly on the 
light quark charge and/or mass. 



Results in this chapter were published in 244]. 



• Finally in chapter [7] we have evaluated the widths for the strong decays 
S* — > Acvr, Sc — > A ^TT and H* HcVr. We use HQS constrained wave functions 
taken from Ref [2o3] and that were obtained solving the nonrelativistic three-body 



problem in a similar way as the one used in chapter O 

The pion emission amplitude has been obtained with the use of PCAC from 
the analysis of weak current matrix elements. Our results are rather stable against 
the quark-quark interaction, with variations in the decay widths at the level of 
6-8%. We find an overall good agreement with experiment for the three reactions. 
This agreement is, in most cases, better than the one obtained by other models. 



Results in this chapter were published in Ref. 2451 ] 



For all single and double heavy hadrons studied in this thesis, we have shown 
that NRQM calculations respect the HQS and HQSS constraints obtained in the 
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infinite heavy quark mass limit. On the other hand the deviations from that hmit 
for the actual heavy quark masses are not always in agreement with HQET predic- 
tions, and we had to include corrections at some instances to improve the NRQM 
results. On the other hand we find NRQM results for many observables agree 
with experiment, lattice results and quark model calculations where relativity is 
explicitly included. 



Appendix A 



S(^X){P) Polarization vectors 



Different sets of polarization vectors used in this work: 
P = 



Spin or helicity bases < 



r <+i)(^) = (0'-73'-7I'0) 

<-i)(^) = (0'71'-71'0) 
I 4)(^) =(0,0,0,1) 



P = \P\k 



Spin or helicity bases < 



'(+1) 
^(-1) 



V2' V2' 

73' ~72' 



-\P\k 



^ <o)(^) =(^,0,0,^) 



f <^,,(P) = (0,-^,-^,0) 



Spin base < 



el,,{P) = {0,^^,-^^,0) 



'(+1) 
'(-1) 



Helicity base < 



-(+1) 



1 ' 0) 
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Matrix elements for the 
leptonic decay of pseudo scalars 
and vector mesons 



The matrix element needed for the evaluation of the pseudoscalar decay constant is given 
by 



P,0 



NR 



Vs2,h{p)l0l5Us^jA-p) 

U-S2J2{P)70US,JA-P) 



d\p\ 

I' 



4Ef,{p)Ef,{p) 



\P\ 



{Ef,{p)+mf,){Ef,{p)+mf,) 
where we have used the fact that Vgj-{p) = 75 u-sj{p)- 



(B.l) 
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Similarly for the vector meson case 



(27r)V2i^/i(p)2i?/.(p) 

^S2,/2(P)73HS1,/1(-P) 



U-S2,h{p) 7375 (-P) 



-73 



1 + 



^|2 



(B.2) 



Appendix C 



Expression for the V^{\q 
matrix element 



In this appendix we give the expressions for the y^(|(f |) matrix elements introduced in 
chapter [3] for the B ^ D decay. They are just a particular case of the ones presented in 
appendix |E] 



V^\\q\) = ^2mB2ED{-q) ( D, -\q\ k 



(JyTiO) 



B, 



NR 



™/2 



2E,{\q\k+p)2Eb{p) 



■■UsiA-W\k - p)l'' u,,^^h{~p) (C.l) 



where /2 represents a light u oi d quark. 
Going a little further we have 



Ec{\q\k + p) + mc] [ Eb{p) + nib 



AE,{\q\k+p) Eb{p) 



( 



\p?+p.\q\ 



^ [Ec{\q\k + p) + nicj [Eb{p) + mbJ ^ 
In the case of equal masses rri}, — rric and for \q \ ~ {w ~ 1) we will obtain 



(C.2) 



(C.3) 
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Similarly 

VH\q\) = 



[Ec{\q\k + p) + m}j (^Eb{p) + mb 
\ 4E,{\q\k+p) Eb{p) 

„ / Pz , Pz + \q\ 



Eb{p)+mb Ec{\q\k + p) + rric 



(C.4) 



Appendix D 



Expressions for the ^ and 
\ matrix elements 



In this appendix we present expressions for the V^^ and matrix elements introduced 
in chapter [3] for the semileptonic B ^ D* decay. They are just a particular case of the 
ones presented in appendix IeI 



\q\k 



B, 



NR 



TO/2 



1 



2E,{\q\k+p)2Eb{p) 



siA-W\ k~p)jf,Us'^^i,{-p) 



- ^2mB2ED4~'l) nb{d*^^ -kl^ 
= ^2mB2Eo^i-q-) J d'pY. (^Zc^,s..fJ 



(^a)m(O) 

TO/, 



B, 0^ 



NR 



1 



2E,{\q\k+p) 2Eb{p) 



\q\k+p)] T.€Uis.J.)iP) 

(D.l) 
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)(| ""^^ \q\k+p\)] 



Ee{\q\k + p) + m, 



Eb{p) +mbj 



'iEe{\q\k+p) Eb{p) 

Pz Pz + \q\ 



V2 



Eb{p) + mb Ee{\q\k+p) 



1 (liD')n "^f2 



sj2mB2ED'{-q) / d 



\q\k+p\)] 



(D.2) 



A 



1 + 



^pI-\p\^ -Pz\q\ 



y {Ec{\q\k+p)+mcj lEb{p) + mbj J 



(D.3) 



rUc + m/2 



\ ^E,{\q\k+p) Ebip) 

/ Pz , Pz + \q\ 



+ 



\Eb{p)+mb Ec{\q\k+p) 



(D.4) 



^^--^^-(-^"^ / i^-i^^+^^i ) 



Eci\q\k+p) +mc] [Eb{p)+inb 



( 



AE^{\q\k+p) Eb{p) 



2pl-\p\^+p,\q\ 



y iEe{\q\k+p) + mc\lEb{p)+mb\J 



(D.5) 



Appendix E 



Expressions for the 

vnm, v^^^^{\q 

AI'{\q\),A';,Mq\) 



and 



A 

matrix elements 



4^ 



q 



E.l Transitions involving quarks 



Here we give general expressions valid for transitions between a pseudoscalar meson Mj 
at rest with quark content Qfjlf^ a-nd a final Mp meson with total angular momentum 
and parity = 0~, 0"*", 1~, 1"*", 2~, 2''", three-momentum — |(f |fc and quark content Ifjlf^- 
In the transition it is the quark that changes flavor. The phases of the wave functions are 
the ones chosen in Eqs. (|4.1l4.2l4.3p . We generally have 



V''(|g I) - A^^m) - ^2mi2EF{-q) ( Mpir). A - \q\ k 



NR. 



1 



2^ 2^ V{S'^J[),(S2J2)^PI) ^(suh),(s2j2)yP 



Mi{Q-), ) 



^PEj72Ei, 



S', Sl,S2 



NR. 
q\k) 



-\q\k~p) -l5)Us,.fA 



™/2 



\q\k-p) 



(E.l) 



where V^'{A^') represent any of the (A^), V^'^^ or F^^^^ (^t(A))' where 

and Ef^ are shorthand notations for ^/{ (- |g|fc-p) and |g|fc-p) 

respectively. Defining also -By = By and By = By -l-TOy we arrive at the following 

final expressions: 
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• Case J"" = 0- 



( 



(0-)) 



™/2 



my; + 



l<flfc 



1 + 



(0-)) 



( ■^Si 



iJ Li - T). \ 



(E.2) 



• Case r = 0+ 



(0-)) 



P-,:r^\S 



^^■,^^1 [ P-{ m^!+mf, W\k-P) ^P-{- rn,Zl,^ \q\k - p) 



(o~)) 



Ef[Efi 



P. 1 



EfiEfi 



EfiEf^ 
( ^\q\-Pz) P- 



+( 



m/; + TO/2 



m/j + TO/2 



ruf^ + ruf^ 



\q\k-p 



(E.3) 



Case J"" = 1- 



I? I 



+ 



'"/(+™/2 



\q\-Pz 



fi 



E 



fi 



(E.4) 
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and similarly 



\q\-Pz 



P ^W\k 

'-^\q\-pA 



+ ■ 



%/2 



TO/i + 



\Q\k 



(0-)) 



Ef^Ef^ 



P ^W\k 



I 



^/(+'"/2 



Ef^Ef^ 

7 



(E.5) 



Case J'^ = 1+ 



f(o-)) 



^^W\k 
ruf^ + ruf, 

ails ul^_^T•^^ 



Ef[Ef^ 



^2mi2EF{-q) / d^p 



47r|p| V^/^/: 



;(Mf(i+,s„-=i)) 



(0-)) 



TO/^ + TO/2 



\ ^Ef^Ef^ Ef,Efi 
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mf> + ruf^ 



^fl^n 2 



1 1 
+ 



A=-l 



(k1) 



P ^lalfc 



Ef^Efi ( Py+Pz+PzW\ mf,+mj^ Pl+Pt-PzW\^^ 



E 



E 



(l+,S',,-=0)3 



A=0 



(l9|) 



(0-)) 



/i 



™f2 



Ef^Efi I rnji +nif^ 

'A — E 



-Pz 



+ ■ 



(l + ,S„-=l)),|^|A* 7(M,(0-)) 



(IpD) « 



m/2 g 



(E.6) 



and similarly 



5,5=0)1 



m) = 



sj2mi2EF{-q) / d'^'p 



47r|p I 



Ef^Efi 



^(Mf(1+,S,, 



=0)) 



(0-)) 



™/2 



mf> + ruf^ 



\q\k 
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.V3 



i—sj2mi2EF{-q) 



-\v\ Vf^'f- 



Att\p\ 



7 (MHO-)) 

/li /2 



ruf^ + ruf^ 



\q\k 



+ 



mf^-irif^ pl\q\ 



(E.7) 



• Case = 2- 



p — ^lalfc 



2 



Ef[Efr p,{pl+pl)\q\ 



\ ^Ef^Ef, Ef,Efi 



= i^^2mi2EF{-q) / d^p 



4Tr\p\ 



(0-)) 



P ^W\k 



Efi^Efi 



[pI-pD 



lo-l -V + \a\\ 



Et 



J 



-PzPv 



Ef^ Ef 



:^V2mr2EFi-q) I d^p ^ (C/f C/f" 



P ^W\k 



Ef^Efi 

\\^eJ7eJ; 



Pz{pI+pI) 



(E.8) 
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and similarly 



(0-)) 



\q\k 









2 2\ 




- Py) 


\ 



-VzVxm 



(E.9) 



Case r = 2+ 



= i ■ 



.^/3 



|g|fc 



/ Pi-Pi-PzW\ raf,+L,, Pi-Pi+PzW\ mf,'+m,, 



(E.IO) 



and similarly 

<i.°(lfl) = 



(0-)) 



m/; + ruf, 



-\q\k 



EfiEfi (p'+pI- ^p' - ^P^WlT^^jritj; pI+pI- ^p' + iTJ^fp^ 



E 



l(Mi{0-)) 

^fl ■ /2 



Ef 



\q\k 



E 

^4^ 



+ 



H 



+™/2 



Ef[Efi 



E.2. TRANSITIONS INVOLVING ANTIQUARKS 



149 







1 ( 




47r|p ' 


4Ej 


Efi 


{-{ 



do-)) 



1 ff 

^ (- m,,+m,, l9 - P) ■ ( m,"+m,, |g 1^ " P) 

Efi^Efi 



V 



Ef,Ef-i 



2p.( |g| - P.) • {—^ |g| -Vz) 



+ {pI+pI) (-Pz + 



ruf^ + mf^ 



2(m// +m/J ^ 



(E.ll) 



E.2 Transitions involving antiquarks 



When it is the antiquark that suffers the decay the expressions are modified as de- 
scribed below for a general transition between a pseudoscalar meson Mj at rest with 
quark content Qfjlf^ a-nd a final Mp meson with total angular momentum and parity 
J'^ = 0^, 0+, 1^, 1+, 2~, 2+, three-momentum — |^|fc and quark content qfiQf^- We have 



MiiO-), 



NR 



(_1)S2-S2 



S'^ Sl,S2 



\q\k+p) 7''(7 - 75) 



\q\k+p) 
(E.12) 



where V^'{A^') represent any of the V^' {A^'), V^^^ (A^^^) or 1/^^^^ (^t(a))' ^^d Ef^, Ej^ 
are shorthand notation for Ef^{—— — T^^lqlk + p), Ef^{ ^ "^im , 1^1^ + We can use 
now that 

VsAp) = {-ir'^^-^Cul0) ; = -(-l)(V2)-« „T^(p^) ct (e.13) 

where C is a matrix given in the Fermi-Dirac representation that we use by 

C = iiV (E.14) 

and that satisfies 



(E.15) 
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Using the above information and making the change of variable p — > —p we can rewrite 



V»(l<ll)-^"(lfl) = -/2vi,2Er(-g) 

•' S' S1.S2 ^ ^ ll ' 1-2 



^ ^4 (- ^ ""iL - P ) - 75) Us2j2 { _ \q\k - p) 



(E.16) 

By comparison with the corresponding expressions involving quarks we find that, apart 
from the changes in the masses involved, there is an extra minus sign for the vector part, 
and, due to Clebsch-Gordan re-arrangements and the fact that Yijn{—p) = (—1)' Yim{p), a 
global sign given are the orbital and spin angular 

momenta of the initial (final) meson. 

In any case this implies a change of sign in the relative phase between vector and axial 
contributions, which in its term produces a sign change in the tensor helicity components 
combination Hp due to the fact that H+i+i goes into TC^i^i and vice versa. All other 
tensor helicity components combinations defined in Eq. (|4.36p keep their signs. 

A simple way of anticipating the above result is the following: the current for qf^ 
decay into g/^ is 

^/,(0)7''(/-75)*/^(0) (E.17) 

But for antiquarks the fields that play the similar role as the ^ fields play for quarks are 
the charge conjugate ones In terms of the latter the above current is written as 

(0)7^(7- 75)*/^ (0) =W;,(0)7n-^-75)*5,(0) (E.18) 



Now this is similar to the current for quark decay but with an extra minus sign in the 
vector part. Whatever other changes might come from reorderings in the wave functions 
we will have an extra relative sign between the vector and axial part. 



Appendix F 

Expressions for the helicity 
components of the hadron 
tensor 



In this appendix we give the expressions for the non-zero hehcity components Ti^s of 
the hadron tensor, as defined in Eq. (|4.32p . corresponding to a cc transition. The 

different cases correspond to the ones discussed in the main text. 

• Case 0" ^ 0-,0+ 



ml, - mt. 



■Hto{PB^,Pcc) = notiPB^,Pcc) = X'^^Q^rnl^^rnl^) 



2 2 



noo{PB^,Pcc) 

• Case 0" 

HttiPB^^Pcc) = 



q- 

X{q^,m%,mls) 



(F.l) 



nt {Pb^ ,Pcc) = not {Pb^ , Pec) 



({niB. nic) {Aoiq') A+{q')) A^{q')) 

\ ms, + nice J 



(mB^-mcc) (Aoiq^) - A+{q^)) 



ms, + rricc 



A-{q^ 



{mB, - mcc) 



^X-q'-mZs , A(g2,m|^,m2g) A+{q^) 



H+1+i{Pb^,Pcc) 
H-1-i{Pb^,Pcc) 



■Aoiq' 



+ nice 
niBr^ + nice 



V{q^) + (ms, -TOcc) Ao(g^) 



Viq^) + iniB,~nicc) Ao{q^) 
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APPENDIX F. HELICITY COMPONENTS 




) 



(F.2) 



For a B~ B transition (with B representing any of the B — B , B , B^ , B*^), 
where it is the c antiquark that decays, we have to change the mass of the final meson in 
the expressions above and take into account the changes in the form factors that derive 
from the discussions in appendix IE. 21 



. Case 0- ^ 2-,2+. 

H»(Pb.,Pcc-) = ^'^';"f'r'^'^ {T,{q') + {ml^-ml,)T,{q') + q'n{q'))' 

Ut oiPB^, Pec) = no t {Pb^ , Pec) 

= ^"^'^hf'"^ W) + {ml^~ml,)T,{q^)+q^T,{q^)) 

X ( {m\, -q'- ml,) T^{q^) + X{q\ml^,ml,) T^iq^) ) 
H+i+i(Pb.,Pcc-) = ^^'"'"V""''^ (T,{q')~X'/'{q',ml^,ml,}T,{q')y 

H-i-i(Pb.,Pcc-) = ^^"''g'tgr''^ {T,{q') + \'/'{q',ml^,ml,)Uq')y 
noo{PB.,Pe.) = ^^''I'^t'T''^ {{ml^-q'~ml,)T,{q') + X{q',ml^,ml,)T,{q'))' 



(F.3) 



Appendix G 



Variational wave function 
parameters for doubly heavy 
baryons 

In Tables IG.ll (doubly heavy S baryons) and IC2I (doubly heavy il baryons) we give the 
variational parameters of the orbital wave functions evaluated with the different quark- 
quark interactions analyzed in this work. 
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bi 


di 


Ci2 


b2 


d2 




kf 


d[i 


(24 


f)4 


d^ 


^fi-] 




^cc 


ATT 


1.78 


0.37 


0.38 


1.04 


0.28 


0.57 


1.24 


0.66 


1.36 


1.56 


0.08 


0.44 


0.44 




AL2 


2.48 


0.09 


1.25 


1.50 


0.01 


0.34 


1.55 


1.46 


2.20 


2.21 


1.82 


0.45 


0.45 




API 


2.14 


0.14 


0.83 


1.45 


-0.01 


0.67 


1.95 


1.77 


2.69 


2.98 


1.34 


0.40 


0.40 




A P9 


1.86 


0.11 


0.48 


1.25 


0.14 


0.62 


1.62 


1.17 


1.92 


2.20 


1.26 


0.41 


0.41 




BD 


1.71 


0.18 


0.37 


1.23 


0.16 


0.85 


1.59 


0.96 


1.46 


2.07 


1.22 


0.44 


0.44 


"* 

'CC 


ATI 

ALl 


1.26 


0.40 


-0.03 


0.82 


-0.30 


0.55 


0.89 


0.54 


0.76 


1.66 


0.15 


0.83 


0.83 




A T O 


1.18 


0.53 


-0.02 


0.65 


-0.36 


0.58 


0.85 


0.58 


1.29 


1.71 


0.12 


0.83 


0.83 




API 


1.93 


0.17 


0.66 


1.37 


0.00 


0.67 


2.29 


1.56 


1.66 


2.99 


1.28 


0.76 


0.76 






1.18 


0.46 


-0.02 


0.64 


-0.38 


0.43 


0.79 


0.65 


1.59 


1.73 


0.12 


0.75 


0.75 




BD 


1.23 


0.31 


-0.01 


0.72 


-0.22 


0.53 


0.97 


0.63 


1.22 


1.79 


0.17 


0.73 


0.73 




ATT 

ALl 


2.37 


0.23 


1.26 


1.31 


0.85 


1.47 


1.93 


1.04 


1.17 


2.09 


2.13 


0.39 


0.39 




AL2 


2.35 


1.81 


1.54 


1.57 


1.41 


2.03 


2.11 


0.27 


0.69 


1.85 


2.20 


0.39 


0.39 




API 


2.66 


1.52 


1.57 


1.48 


1.42 


1.99 


1.94 


0.37 


0.85 


1.77 


2.31 


0.35 


0.35 




A P9 


2.22 


1.93 


1.63 


1.53 


1.59 


2.01 


2.23 


0.19 


0.69 


1.82 


2.08 


0.36 


0.36 




BD 


2.74 


1.35 


1.42 


1.42 


1.20 


2.04 


1.90 


0.43 


0.85 


1.80 


2.54 


0.40 


0.40 


" * 


ATI 


3.17 


0.15 


1.04 


1.89 


0.24 


0.14 


1.47 


1.79 


2.25 


2.84 


1.47 


0.54 


0.54 


A T O 


2.82 


0.11 


0.56 


1.56 


0.39 


0.57 


2.14 


1.47 


1.47 


2.80 


1.89 


0.52 


0.52 




API 


2.69 


0.18 


0.88 


1.73 


0.33 


0.75 


2.33 


1.80 


1.90 


3.14 


1.76 


0.48 


0.48 




AP2 


2.78 


0.08 


0.41 


1.60 


0.33 


0.71 


2.05 


1.45 


1.56 


2.87 


1.91 


0.49 


0.49 




BD 


3.40 


0.13 


0.19 


2.74 


-0.15 


-0.07 


2.13 


2.04 


2.64 


3.85 


3.02 


0.51 


0.51 


"6c 


ALl 


2.24 


0.21 


0.39 


1.73 


0.01 


-0.01 


2.18 


-0.78 


1.62 


0.77 


1.84 


0.01 


1.51 


AL2 


2.40 


0.12 


0.28 


1.67 


-0.06 


-0.01 


2.15 


-0.82 


1.29 


0.87 


1.49 


0.01 


1.53 




API 


4.81 


-0.01 


2.78 


2.03 


0.05 


-0.15 


3.45 


0.52 


2.31 


0.91 


1.09 


0.01 


1.37 




AP2 


2.85 


0.05 


0.98 


1.91 


-0.03 


-0.15 


2.66 


-0.19 


0.68 


0.99 


1.55 


0.02 


1.34 




BD 


2.60 


0.14 


0.58 


1.62 


0.05 


-0.04 


2.47 


-0.66 


1.65 


0.94 


1.28 


-0.03 


1.54 


S6c 


ALl 


2.82 


0.15 


0.46 


2.37 


-0.15 


-0.09 


1.06 


-0.37 


1.31 


0.70 


1.35 


-0.02 


1.09 




AL2 


2.97 


0.08 


0.57 


2.14 


-0.11 


-0.05 


1.15 


-0.40 


0.92 


0.72 


1.39 


-0.06 


1.20 




API 


2.91 


0.10 


0.56 


2.16 


-0.10 


-0.08 


1.13 


-0.50 


1.21 


0.72 


1.25 


-0.06 


1.05 




AP2 


4.75 


-0.02 


1.98 


2.74 


0.06 


0.38 


1.47 


-0.16 


3.03 


2.13 


-0.04 


-0.07 


1.09 




BD 


2.64 


0.21 


0.48 


2.38 


-0.08 


-0.12 


1.02 


-0.57 


1.25 


0.56 


1.74 


-0.07 


1.20 


"* 
"6c 


ALl 


1.85 


0.21 


-0.19 


2.20 


0.10 


-0.02 


1.26 


-0.66 


1.27 


0.70 


1.71 


0.07 


1.72 


AL2 


1.93 


0.13 


-0.14 


2.25 


-0.09 


-0.05 


1.10 


-0.56 


0.92 


0.62 


1.67 


0.06 


1.74 




API 


1.93 


0.17 


0.06 


2.54 


0.17 


-0.10 


0.97 


-0.57 


1.07 


0.55 


1.76 


0.06 


1.56 




AP2 


2.16 


0.07 


-0.12 


1.91 


-0.08 


-0.09 


1.23 


-0.78 


0.84 


0.66 


1.14 


0.07 


1.48 




BD 


2.79 


0.12 


0.55 


2.14 


-0.05 


-0.03 


1.65 


-0.83 


2.08 


0.79 


1.39 


-0.01 


1.70 



Table G.l: Variational parameters of the three body orbital wave function for 
doubly heavy E baryons. a's are dimensionless, d's have dimensions of fm and b's 
and a's have dimensions of fm~^. Whenever the two heavy quark are different the 
suffix hi stands for a b quark and /i2 stands for a c quark. 
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di 


(22 


62 


d2 


0,3 


63 




(24 


64 


C/4 








ALl 


1.60 


0.29 


0.30 


1.29 


1.20 


1.70 


1.27 


0.25 


0.94 


0.85 


1.33 


0.53 


0.53 




AL2 


1.21 


0.17 


0.25 


1.35 


1.22 


2.04 


1.66 


0.22 


0.69 


0.98 


1.63 


0.48 


0.48 




API 




nil 




1 70 


1 26 


1 79 
J_. 1 z 


1 .ou 


12 


U.Uo 


u.ou 










A P9 


1.63 


-0.24 


0.62 


1.75 


2.01 


2.86 


2.16 


0.10 


0.16 


1.35 


2.56 


0.49 


0.49 






1.34 


0.09 


0.37 


1.46 


1.30 


1.91 


1.69 


0.28 


0.73 


1.17 


1.55 


0.47 


0.47 




ALl 


1.29 


0.16 


-0.01 


0.90 


-0.49 


0.96 


0.99 


0.75 


0.67 


2.05 


0.13 


0.95 


0.95 




AL2 


1.41 


0.10 


0.01 


0.72 


-0.47 


0.99 


1.00 


0.77 


0.71 


2.00 


0.15 


0.89 


0.89 




A PI 




n 1 n 


U. -LU 






L.OO 


1 .uu 


u.ou 


-U.UO 


9 (\R 


41 


n QA 
u.y^ 


n QA 




A P9 


2.32 


0.04 


0.22 


1.97 


-0.29 


1.33 


1.07 


0.47 


-0.65 


2.31 


1.70 


0.92 


0.92 






1.66 


0.15 


0.16 


1.41 


-0.27 


1.14 


1.26 


0.80 


0.30 


2.21 


0.71 


0.73 


0.73 




ALl 


1.28 


0.91 


-0.07 


0.71 


1.05 


1.86 


2.32 


0.26 


1.38 


0.99 


1.27 


0.51 


0.51 




AL2 


1.20 


0.82 


-0.04 


0.65 


1.28 


1.78 


2.56 


0.15 


1.22 


1.03 


1.41 


0.45 


0.45 




A PI 


1 1 

L. lo 


n on 
u. yu 


n no 


n 79 


1 24 


1 70 


2 48 


n 9n 


1 A'^ 


1 12 


1 12 


44 


44 




A P9 


2.00 


0.03 


0.22 


0.58 


1.70 


2.81 


3.22 


0.06 


-0.63 


0.42 


3.32 


0.47 


0.47 






1.47 


0.58 


-0.09 


0.79 


1.16 


1.68 


2.48 


0.25 


1.25 


1.00 


1.52 


0.46 


0.46 


^Ib 


ALl 


1.14 


1.01 


0.01 


1.03 


0.44 


3.05 


2.03 


0.36 


1.20 


2.48 


1.37 


0.68 


0.68 


AL2 


0.97 


1.30 


0.08 


0.98 


0.75 


2.79 


2.36 


0.19 


1.21 


2.40 


1.38 


0.63 


0.63 




API 


1.05 


1.16 


0.11 


1.29 


0.56 


3.25 


2.16 


0.29 


0.86 


2.54 


1.31 


0.63 


0.63 




AP2 


0.94 


1.30 


0.33 


1.11 


0.76 


3.39 


2.55 


0.11 


0.79 


2.45 


1.75 


0.65 


0.65 




BD 


1.52 


1.12 


0.19 


1.15 


0.74 


2.91 


2.05 


0.37 


1.29 


2.44 


1.41 


0.57 


0.57 




ALl 


2.20 


0.32 


1.00 


1.92 


0.05 


-0.09 


2.89 


-0.24 


1.79 


0.75 


1.87 


-0.14 


2.42 


AL2 


2.21 


0.22 


0.76 


2.03 


-0.07 


-0.16 


2.94 


-0.30 


1.56 


0.78 


1.92 


-0.18 


2.34 




API 


5.34 


0.04 


2.24 


2.83 


0.01 


-1.34 


4.37 


0.29 


0.64 


2.14 


-0.28 


-0.17 


2.35 




AP2 


3.33 


0.01 


1.71 


2.14 


-0.04 


-0.36 


3.00 


-0.20 


1.23 


0.71 


1.88 


-0.17 


2.34 




BD 


3.06 


0.14 


0.83 


1.94 


0.18 


0.44 


2.16 


-0.11 


1.67 


0.75 


1.78 


-0.24 


2.26 




ALl 


2.15 


0.14 


-0.04 


1.22 


-0.32 


-0.02 


2.33 


-1.04 


1.10 


0.84 


0.98 


-0.20 


1.84 




AL2 


2.03 


0.15 


0.03 


1.39 


-0.36 


-0.01 


2.37 


-0.81 


1.22 


0.75 


1.78 


-0.25 


1.84 




API 


1.89 


0.18 


-0.02 


1.49 


-0.40 


0.04 


2.10 


-0.61 


1.01 


0.81 


1.90 


-0.24 


1.78 




AP2 


2.34 


0.02 


0.51 


4.50 


-0.01 


0.53 


1.82 


-0.08 


3.95 


2.69 


5.47 


-0.23 


1.78 




BD 


2.17 


0.18 


0.15 


1.50 


-0.06 


-0.01 


2.46 


-0.89 


1.41 


0.75 


1.97 


-0.27 


1.85 




ALl 


2.59 


0.10 


1.01 


1.90 


-0.01 


-0.27 


2.88 


-0.13 


0.87 


0.51 


2.88 


-0.06 


2.54 


AL2 


2.69 


0.11 


1.15 


1.65 


0.10 


0.08 


2.80 


-0.07 


1.28 


0.64 


2.90 


-0.13 


2.51 




API 


2.59 


0.14 


0.96 


1.88 


-0.01 


-0.08 


2.82 


-0.23 


1.15 


0.52 


3.09 


-0.10 


2.59 




AP2 


3.80 


-0.01 


1.75 


1.81 


0.03 


0.69 


2.31 


0.03 


2.48 


0.56 


2.94 


-0.11 


2.60 




BD 


2.67 


0.12 


1.25 


1.88 


0.06 


-0.21 


2.97 


-0.11 


1.15 


0.52 


3.16 


-0.20 


2.35 



Table G.2: Same as Table lUH for doubly heavy Q baryons. 



Appendix H 



Expressions for the ^1,^2,^3 
and Gi^G2^G^ form factors 



In this appendix we study the vector Fi, F2, -Fa and axial Gi, G2, G3 form factors, defined 
in chapter [H] to study semileptonic decays of doubly heavy baryons. 

H.l Form factors in terms of the T^'^(\q\) and JC^'^{\q\) 
integrals 

In this section we relate the vector Fi, F2, F3 and axial Gi, G2, G3 form factors that we 
evaluate in the center of mass of the decaying baryon, to the integralsX-^ '^(I'f l)i '^(I'f I) 
defined in Eq. (j6.15|) . To simplify the expressions it is convenient to introduce 



Bfi^n EB'{-q) + mB' / 2i?c(g) , 



rrir 



n EB'{-q) + mB' / 2i?c(g) „ .|^|. . „ „ 

y 2EB'[~q) y Ec[q) + rric 

(H.l) 

• Cases Sh = I, S',^ = or Sh = 0, S;^ = 1: 

hm 2 / /C^'^(|g|) / me : ^ 



= 



EB'{-q) + mB' Vi\Ec{q)+mc 2 \E^{q) mt 

Fii\q\) , hi\q\) 

EB'{-q)+mB' mB' 



F^i\q\) + F2i\q'\) + ^^^^F,im = 
rriB' 

(H.2) 

where Ec{q) ~ \/ m\ + (p- . From the above expressions we have that F2 — F^. 
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mB'{EB'{-q) +mB') 



-gl(|g|) + G2(|g|) + ^^g3(|g|) ^ 

EB'{-q) +mB' 



+ 



V3 \^£;c(g)+mc 2 V^^c(«) rnb 

Case = 1, 5; = 1 



EB'{-q) + mB' 3 \^£;c(g)+mc 2 \E^{q) 

FA\q\) _ 



EB'{-q) + mB' niB' 

2 / I^'^iWl) /C^'^(|g|) (ml 
E,{q)+me 2 V^^2(9) 



F,m+F2m + ^^^Fs{\q\) 

TUB' 



(H.3) 



(H.4) 



mB'{EB'{-q) + mB') I 

3 I 1^ ^ 2 (£;e(9) + me) V£^?(9 ) ^6 



-Gi(|g|) + G2(|<r|) + ^^G3(|9|) 



EB'{-q) +mB' 
4 / 



T^^(|g|) /C^^(|g|) / 



3 I Ee{q)+mc 2 \E^{q) nib 



(H.5) 
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H.2 I^'^{\q\) and /C^'^(|g |) integrals 

To evaluate ^iWl) and K,^ we use a partial wave expansion of the orbital wave 

functions: 

oo 

*f,,(ri,r2,ri2) = 5]/f(ri,r2)Pz(M) 

oo 

*^Uri,r2,ri2) =^/,^'(ri,r2)PKM) (H.6) 



where /i is the cosine of the angle made by fi and f2 and Pi (/i) is a Legendre polynomial 
of rank /. The radial functions fP{ri,r2), fP (?'i,f2), are evaluated as 



21 + 1 



= dA*PKM)*f..(r-i,r2,ri2) (H.7) 



In terms of those we have 



; ^7:;^ Jo M Jo M 

x/f (ri,r2)/,^(ri,r2) (H.8) 

with j' s being spherical Bessel functions. 
And 



B'Bi 

E E E E E("i)^''^'"^'^/ V(2i + 1)(2;" + mv" + 1) 



/C^^(|g|) 

16^2 



X (/,/',;"! 0,0,0) (?",r,i| 0,0,0) (/',L,r| 0,0,0) w{i"i'iL-ii"') 

X / rfrir^j,»( "'lZ"^ |g>i) / dr2r22jr(^k>2)/f (ri,r2)J7^/f(ri,r2) (H.9) 
Jo M Jo M 

with W(l"l'lL;l,l"') being a Racah coefficient and fl^ the differential operatoiQ 



For the actual evaluation we restrict the I, I' values to /, = 0, • • • , 6. 



Note that the Racah and Clebsch-Gordan coefficients restrict L to the two possible values 
L = l±l. 



Appendix I 



Expressions for the ^^^/^^ 



matrix elements 



The values for the ^^'b^^'^ evaluated using one-body current operators and their 
expressions are given by 



41/2,1/2 



fc°(Qi - "^^^ m k, Q2 + =5^ w\ k) 



M 



At 



M 



A+ 



S ( 2' 2'^ 



1 



2i;d(Qi - |g I k)2Eu{Qi) 



si, -si,0 

{Qi - \q\ k) iQi) 



(1.1) 



With = m„ + md + nic- For = 0, 3 we get the final expressions 



,1/2,1/2 



M 



At 



(EaiQi - \q\ k) + md) {Eu{Qi) + to„) 
\ 2Ei{Q^-\q\k)2Eu{Qi) 

Ql , Ql- W\ 



Eu{Q-i) + mu Ed{Qi-\q\k)+mdJ 
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(1.2) 
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41/2,1/2 _ 

^E++A+,3 ~ 



At 



m 



A+ 



[EdiQi - \q\ k) + md) (k(<3i) + m„) 



X 1 + 



2Sd(Qi - |g I fc)2K(Qi) 
1 



(2{Qlf-\Q,\^-Ql\q\) 



(1.3) 



For we just have 



ll/2,l/2 



= -V2 A, 



1/2,1/2 



(1.4) 



Finally 



4I/2-I/2 _ 



^ I d^gi d^g^ 0f;=^(Qi, Q2) Ul\Z'M - '^^=^ W\ ^' ^2 + ^ W\ k) 



5:(i,i,i.i,-.i,o) (1,1,0 



1 



^j2Ea{Qi-\q\k)2Eu{Q^) 



si,-si,0 
I'd Si {Qi - \q\ k) iQi) 



(1.5) 



with M=o = + mg + mg. For /x = 0, 3 we get the final expressions 



41/2,1/2 _ 



9| fc, Q2 + |g| A:) 

M =0 



(£^d(gi - I? I fc) + TOd) (Eu{Qi) + m„) 
^ 2^d(Qi - \q\ k)2E^{Qi) 

Ql , Ql- W\ 



Eu{Qi)+mu Ed,{Qi - \q\k) + md, 



(1.6) 
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?.3 

V2 



2£;d(Qi - |g I k)2Eu{Qi) 



X 1 + 



7 z z ^7 z V (mif - \Qi? - Ql\<i\) 

{Ed{Qi-\q\k)+md){E^{Qi)+mu)^ \ 



(1.7) 
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